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20.  ABSTRACT  (contd) 

The  matched  filters  were  evaluated  using  complex  terrains  such  as  woodlands, 
roadsides,  and  villages  in  which  the  principal  targets  (i.e.,  tanks,  trucks  and 
armored  personnel  carriers)  were  located.  For  each  terrain  and  target,  the 
signal-to-clutter  ratio  (s/c)  was  measured.  These  ranged  from  5.5  to  11  dB  for 
the  target-terrain  combinations  studied.  The  complex  matched  filters  that 
achieved  these  s/c's  were  recorded  in  dichromated  gelatin. 

'  t  >  . 

It  was  shown  that  more  than  5  matched  filters  could  be  non-coherentlv  added  and 
15  could  be  coherently  added  in  the  same  memory  position.  In  addition,  it  was 
demonstrated  that  carrier  multiplexing  could  be  used  to  achieve  .10  matched  filters 
per  memory  position.  It  is  also  shojurn  that  the  theoretical  limit  for  memory  den¬ 
sity  of  matched  filters  is  2670  MF/cm^.  This  is  equivalent  to  1,000,000  bits  of 
digital  data.  Various  modificationsof  the  basic  optical  matched  filter  were 
studied  for  the  purpose  of  improving  s/c  ratio  and  increasing  the  storage  density. ' 

A  mathematical  model  of  the  optical  matched  filter  and  its  variations  was  developed 
based  on  the  two-dimensional  Fourier  Transform.  This  model  was  used  to  simulate 
the  optical  device  digitally.  The  output  of  me  digital  simulation  showed  excellent 
agreement  with  the  optical  device,  using  both  visual  and  quantitative  measures. 
The  simulation  allowed  for  the  parametric  study  of  the  optical  device  under  con¬ 
trolled  conditions  and  in  much  less  time  than  could  have  been  achieved  in  a  labora¬ 
tory  setting.  Several  digital  studies  aimed  at  optuni/.ing  the  detection  capability 
of  the  correlator  are  presented,  and  the  results  are  compared  with  those  obtained 
from  the  optical  correlator. 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  f»AOC(W*w>  Dmm  CnlorwO 


Grumman  Research  Department  Report  RE- 634 


INVESTIGATION  OF  LARGE  CAPACITY  OPTICAL  MEMORIES 
FOR  CORRELATOR  APPLICATIONS 


Kenneth  G.  Leib  and  Jay  Mendelsohn 
Research  Department 
Grumman  Aerospace  Corporation 
Bethpage,  New  York  11714 


October  1981 


Final  Technical  Report  1  April  1978  -  31  March  1981 


Approved  for  Public  Release 
Distribution  Unlimited 


Prepared  for 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 
Bolling  Air  Force  Base 
Washington,  DC  20332 

AIR  FORCE  OFFT  (TE  OF  SCIENTIFIC  RESEARCH  ( AFSC ) 
NOTICE  OF  «VJ*S»ITTAI.  TO  DTIC 
This  technics!  report  Ins  been  revlev*d  ?nd  is 
approved  for  p ic  ~ele use  IAW  AFR  130—12. 
distribution  is  unlimited. 

ItATTHBf  J.  KEFPSR  / 

Chief,  Technical  Information  Division 


TABLE  OF  CONTENTS 


Section  Page 

1  INTRODUCTION  . . 1-1 

2  DESCRIPTION  . 2-1 

2.1  Experimental  Setup  . 2-1 

2.2  Simulation  Validation  . 2-1 

2.3  Matched  Filter  Formation  and  Playback  . 2-11 

3  TECHNIQUES  FOR  RECORDING  MEMORY  ELEMENTS . 3-1 

3.1  Geometrical  Configurations  . 3-1 

3.2  Split  Spectrum  . 3-4 

3.3  Matched  Filter  Additions  . 3-7 

3.3.1  Noncoherent  Addition  . 3-15 

3.3.2  Coherent  Addition  . 3-15 

3.3.3  5X  COH  and  Non-COH  Comparisons  . 3-24 

3.3.4  Multiplexed  Carriers . 3-41 

4  MEMORY  RECORDING  MEDIA  . 4-1 

5  OPTICAL  MEMORY  CAPACITY . 5-1 

5.1  Capacity  Based  Upon  Frequency  Content  of  MF . 5-1 

5.2  Capacity  Based  Upon  Film  Limits . 5-3 

6  SUMMARY . 6-1 

7  REFERENCES . 7-1 

8  PROFESSIONAL  ACTIVITIES  .  8-1 

9  PERSONNEL . 9-1 

Appendices 

A  Use  of  Dichromated  Gelatin  for  Matched  Filters . A-l 

B  Construction  of  Multicapacity  Matched  Filters  Using  Birefringent 

Object  Film . B-l 

C  Kowarschik's  Equations  for  Sequentially  Stored  Gratings  . C-l 

D  The  Fourier  Transform  of  Coherently  Added  Images  . D-l 

E  Additional  Simulations  .  . . E-l 

F  Analysis  of  an  Optical  Matched  Filter  Guidance  System  for 

Tactical  Housing  Applications  . F-l 


v 


LIST  OF  ILLUSTRATIONS 


Figure  Page 

1-1  Optical  Matched  Filter  Image  Correlator . 1-2 

1-2  Details  Of  An  Offset  Memory  . 1-4 

1-3  Terrain  Scenes  Used  In  Program  . 1-5 

1- 4  Cumulative  Distribution  Function  Of  Cross  Correlation 

(Tank  vs.  Terrain  Scenes)  . 1-6 

2- 1  Arrangement  For  Matched  Filter  Fabrication,  S/C  And 

Correlation  Measurements . . . 2-2 

2-2  Matched  Filter  Holder  . 2-3 

2-3  Analog  Version  of  (Fourier  Transforms)2  At  Two  Exposures  .  .  .  2-4 

2-4  Digital  Version  of  (Fourier  Transforms)2  At  Two  Exposures  .  .  .  2-5 

2-5  M-60  Power  Spectrum  Scan  Perpendicular  to  Major 

Tank  Axis  (Ref  8)  . 2-6 

2-6  (FT)2  of  Digitized  M-60  Tank  Image  . 2-6 

2-7  Power  Spectrum  of  Three  Digitized  Terrain  Scenes . 2-8 

2-8  Reimaging  of  M-60  Matched  Filters  (Analog) . 2-12 

2-9  Reimaging  of  M-60  Matched  Filters  (Digital) . 2-13 

2-10  Normalized  Fringe  Spacing  vs.  Spatial  Frequency 

Along  Axis  M-60  Tank  Axis  . 2-14 

2- 11  Fringe  Spacing  At  v  -  0  For  Different  Object  Position 

During  MF  Fabrication . 2-16 

3- 1  Techniques  For  Recording  Memory  Elements . 3-2 

3-2  Geometrical  Placement  of  Memory  Elements  . 3-3 

3-3  Comparison  of  COH/Non-COH  M-60  Tank  (Fourier 

Transform)2  . . 3-8 

3-4  Intraelement  Interference  In  A  4X  Tank  Image  MF  . .  3-9 

3-5  Autocorrelations  of  M-60  Tank  From  5x5  Array 

(All  Identical)  . 3-13 

3-6  Scan  Along  Major  Tank  Axis  In  Correlation  Plane  Using 

5  x  5  Hi  Transform  Lens  . 3-14 

3-7  Non-Coherent  Addition  At  Constant  Carrier  Frequency  . 3-17 

vii 

TWBHM3  PAOI  M  VffUflT 


LIST  OF  ILLUSTRATIONS  (contd) 


Figure  Page 

3-8  Non-Coherently  Added  MF  in  Dichromated  Gelatin  . 3-18 

3-9  Impulse  Response  of  5X  Image  MF  in  Dichromated  Gelatin  . 3-19 

3-10  Scan  Across  Correlation  Plane  Using  5X  Non-COH 

Matched  Filter  (Dichromated  Gelatin)  . 3-20 

3-11  Coherent  Arrays . 3-23 

3-12  Autocorrelation  of  Single  M-60  Tank  Playing  Through 

5X  Non-COH  MF  . 3-25 

3-13  Autocorrelation  of  a  Single  M-60  Tank  Playing  Through 

5X  COH  MF  . 3-26 

3-14  Autocorrelation  Planes  When  Entire  5X  Image  Plays  Through 

Non-COH  and  COH  Made  Matched  Filters  . 3-27 

3-15  Correlation  Plane  of  a  5X  Non-COH  Matched  Filter  Using 

SE  Asian  Terrain  Image  Input  Upon  Playback  . 3-28 

3-16  5X  MF  Translational  Sensitivity  . 3-32 

3-17  X-Positional  Sensitivity  of  IX  and  10X  Matched  Filters  in  Thick 

and  Thin  Emulsions  . 3-33 

3-18  Y-Positional  Sensitivity  of  IX  and  10X  Matched  Filters  in  Thick 

and  Thin  Emulsions  . 3-34 

3-19  Azimuth  Angular  Response  of  Thin  and  Thick  Emulsions . 3-36 

3-20  Rotational  Sensitivity  of  Matched  Filter  . 3-37 

3-21  Correlation  Plane  of  10X  COH  MF  for  Two  Terrains . 3-38 

3-22  A  S/C  Ratio  of  15X  MF  in  Dichromated  Gelatin  Using  3 

Terrains . 3-43 

3-23  15X  Matched  Filter  in  Dichromated  Gelatin  . 3-44 

3-24  Limits  on  Fringe  Angles  in  Thick  Media,  N  =  1.53 . 3-47 

3-25  Kowarschik  Formulation  Applied  to  Or  =  +10°  ’ . 3-49 

3-26  Response  of  Multicarrier  Gratings  and  MF  in  DCG  . 3-51 

3-27  Azimuth  Sensitivity  of  MF  with  Multiplexed  Carrier  . 3-52 

3-28  5X  MF  in  DCG  with  One  and  Two  Carriers . 3-53 

3-29  Response  of  15X  COH  with  Two  Carrier  Frequencies . 3-54 

3- 30  S/C  Ratio  of  a  Multiplexed  Carrier  15X  MF  in  DCG  Using  SE 

Asian  Scene . 3-55 

4- 1  Memory  Recording  Mechanisms  for  Media . 4-2 

4-2  Thickness  Criteria  .  . . 4-3 

4-3  Scan  with  Microdensitometer  of  Higher  Orders  of  Matched  Filter  .  .  4-4 

viii 


/ 


LIST  OF  ILLUSTRATIONS  (contd) 


Figure  Page 

4-4  200X  SEM  of  5X  COH  MF  in  Dichromated  Gelatin . 4-7 

4-5  Efficiency  of  Recording  Medium,  0g  =  0  . 4-14 

4- 6  Efficiency  of  Recording  Medium,  fls  =  ®R  . 4-15 

5- 1  Fringe  Plane  Angles  vs.  Focal  Length  6  g  =  0,  £  /F  =  0  . 5-4 

5- 2  Effect  of  Finite  Resolution  of  Media  on  Matched  Filter 

Weighting  Function  H  . . 5-7 

6- 1  Digital  Version  of  Multiple  Array  Matched  Filter  . 6-3 

6-2  2x2  Array  and  Correlation  Results  . 6-4 

6-3  Optical  Memory  Requirements . 6-5 


LIST  OF  TABLES 


Table  Page 

2-1  Validation  of  Simulation  Through  S/C  Ratio  Measurement . 2-9 

2- 2  Autocorrelation  and  S/C  Ratio  for  M-60  in  Three  Terrains  as  a 

Function  of  the  Spatial  Frequency  Cutoff  for  3rd  Order  MF  ....  2-10 

3- 1  Split-Spectrum  Simulation  Results  . 3-5 

3-2  Accuracy  of  Correlation  Location  with  the  Split  Spectrum . 3-6 

3-3  Autocorrelations  from  Regularly  Placed  IX,  4X  and  8X 

Input  Objects  . 3-11 

3-4  Effect  of  Interference  on  Correlation  and  S/C  Ratio  . 3-12 

3-5  Performance  of  Non-Coherently  Added  MF  in  Dichro¬ 
ma  ted  Gelatin . 3-16 

3-6  Positional  Distribution  of  Autocorrelations  of  Coherent  MF  Using 

Thin  Silver  Halide  . 3-22 

3-7  5X  COH  and  Non-COH  MF  Performance  . 3-29 

3-8  S/C  Ratio  for  COH  and  Non-COH  5X  Matched  Filters  . 3-30 

3-9  Orientational  Sensitivity  of  Coherent  Addition  MF  . 3-39 

3-10  S/C  Ratio  for  Multiple  COH  Matched  Filters  Using  SE 

Asian  Terrain  Imagery . 3-41 

3- 11  Range  of  Carrier  Frequencies  for  a  Spatial  Frequency  Spread 

of  ±9.5  c/mm  . 3-46 

4- 1  Positional  Distribution  of  Autocorrelations  of  Coherent  Matched 

Filters  Using  Thick  and  Thin  Emulsions  . 4-5 

4-2  Matched  Filter  Media  Comparison . . . 4-8 

4-3  Parameters  for  Equations  4-2,  -3,  -4  4-11 

4-4  Average  Efficiency  for  nj  =  .04;  Index  Interval:  1.52-1.56  ....  4-12 

4- 5  Average  Efficiency  for  nj  =  .04  at  N  =  1.55  for  Various  . . 4-13 

5- 1  Unit  Capacity  for  Some  Practical  Cases  . 5-2 


x 


I 


I 


1  -  INTRODUCTION 

There  are  a  number  of  potential  Air  Force  applications  for  which  the  matched 
filter  optical  correlator  could  provide  a  significant  improvement  in  performance. 
Typically  these  might  include  missile  guidance,  smart  sensors,  and  target  identifi¬ 
cation  for  reconnaissance.  The  techniques  available  in  Fourier  optics  and  holo¬ 
graphic  technology  make  the  implementation  of  two-dimensional  matched  filters 
(MF)  for  correlators  a  feasible  approach  to  the  general  problem  of  multiple  pattern 
or  target  recognition  in  any  Air  Force  program  concerned  with  real-time  target 
recognition.  Appendix  F  describes  in  some  detail  one  such  Air  Force  application 
within  the  overall  context  of  this  program. 

The  purpose  of  this  program  is  to  investigate  large  capacity  optical  memories 
for  correlator  applications  (Ref  1).  The  basic  memory  element  in  the  correlator  is 
the  matched  filter. 

The  essential  elements  of  correlators  with  coherent  optical  memories  are 
illustrated  in  Fig.  1-1  (Ref  2).  The  Grumman  innovation  for  storing  large  memo¬ 
ries  (Ref  3)  is  to  utilize  combinations  of  multiple  holographic  lenses  and  MF  arrays 
or  memories.  High  repetition  rate  input  transducers  for  correlators  might  soon  be 
available  as  indicated  by  recent  progress  with  such  devices  as  microchannel  light 
modulators  (Ref  4),  CCD  array  light  modulators  (Ref  5),  and  thermoplastic  trans¬ 
ducers  (Ref  6)  upon  which  targets  have  been  imaged. 

However  derived,  the  Fourier  transform  of  the  object  image  in  Fig.  1-1  is 
replicated  by  a  multiple  holographic  lens  addressed  to  a  coordinated  optical 
memory  bank  of  matched  filters.  An  inverse  transform  lens  and  detector  complete 
the  correlator.  Our  concern  is  with  the  matched  filter  optical  memory  as  part  of, 
and  as  measured  and  evaluated  in  the  optical  correlator. 

The  experimental  approach  taken  in  this  program  is  one  in  which  various 
models  of  a  matched  filter  are  constructed,  and  the  parameters  affecting  their 
performance  determined. 

For  example,  a  MF  (not  limited  by  optical  constraints)  has,  in  principle,  an 
unlimited  spatial  bandwidth  which,  in  any  packing  arrangement,  would  give  rise 
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to  crosstalk  with  contiguous  matched  filters  in  the  memory  bank  (see  Fig.  1-2). 

Yet,  one  can  predict  and  verify  that  a  limited  filter  bandwidth  will  yield  an 
acceptable  signal /clutter  (S/C)  ratio  with  little  or  no  crosstalk  (Ref  7).  This 
optimization  then  gives  rise  to  a  high  packing  density.  Other  parameters  affecting 
filter  performance  interact  with  the  detection  process.  In  this  program,  the  posi¬ 
tion  is  taken  that  large  capacity  optical  memories  are  not  isolated  entities,  but 
large  only  when  the  associated  correlation  signals  can  be  properly  detected  and 
used.  This  is  considered  an  important  aspect  of  Air  Force  applications. 

In  this  program,  the  basic  target  of  interest  is  the  US  Army  M-60  tank,  used 
in  an  overhead  view,  and  used  at  different  times  in  various  terrains.  The  latter 
are  shown  in  Fig.  1-3.  Figure  l-3a  represents  a  "SE  Asian"  scene  with  competing 
vehicles;  Fig.  l-3b,  c  and  d  represent  woodland,  village,  and  roadside  scenes 
respectively.  Other  scenes  are  occasionally  referred  to  and  are  described  when 
used.  All  terrain  images  were  obtained  from  the  US  Army  400/1  scaled  simulation 
model  located  at  the  Night  Vision  Laboratory  at  Fort  Belvoir. 

Figure  1-4  shows  the  result  of  determining  the  cumulative  distribution  of  the 
cross  correlation  between  the  M-60  tank  and  the  terrain  scenes  illustrated  in  Fig. 
1-3,  b-d.  The  cross  correlation  is,  of  course,  a  measure  of  clutter  presented  by 
each  of  the  terrains.  The  abscissa  expresses  clutter  in  terms  of  the  ratio  of  the 
cross  correlation  to  the  autocorrelation.  For  any  particular  value  of  clutter,  the 
corresponding  ordinate  is  the  probability  that  the  clutter  is  less  than  that  value. 
Therefore,  the  more  steeply  rising  curve  present  the  least  clutter.  In  the  figure 
we  see  that  for  the  tank  the  woodland  has  the  least  clutter  and  the  village  the 
most.  This  is  as  expected  since  the  village  has  many  more  straight  edges. 

The  matched  filter  is  the  optimum  linear  filter  for  signal  detection  in  white 

gaussian  noise.  Optically,  the  transfer  function  of  the  matched  filter  is  S*(v  ,v  ), 

x  y 

the  complex  conjugate  of  the  Fourier  transform  S(rx,  i>y)  of  the  desired  signal: 


S(v, 


'V  = 


co 

« 


s(x,y)e 


-i2ir(  v  x 
,ye 


+  v  v) 


(1-1) 


dxdy 


j3u)x 

When  interfered  with  a  unit  reference  beam,  e  ,  it  forms  the  FT  hologram  or 
matched  filter,  MF.  The  size  sensitivity  of  the  filter  is: 

S  =  v  IK  =  (AF)  1  (cycles/mm)  per  mm  (1-2) 

x  ,y 

X  is  the  wavelength  and  F  the  focal  length  of  the  transform  lens  in  Figure  1-1. 
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Fig.  1-2  Detail*  of  an  Offset  Memory 
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Fia.  1-4  Cumulative  Distribution  Function  of  Ctom  Correlation  (Tank  v*.  Terrain  Scenes) 
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This  report  is  an  overall  review  of  the  main  elements  of  the  program.  During 
the  program,  the  philosophy  was  to  give  many  presentations.  This  extensive 
effort  is  summarized  in  Section  8.  Many  of  the  details  of  these  presentations  have 
been  presented  as  papers  during  the  course  of  the  program  and  are  referred  to  at 
appropriate  times  in  the  report.  Some  are  included  as  appendices  in  those  cases 
where  publication  might  otherwise  have  been  restricted.  Much  data  is  cited  to 
illustrate  trend  of  the  results. 

Section  2  is  an  extensive  description  of  the  actual  system  used  and  the 
matched  filter.  Our  approach  to  the  simulations  is  also  described  and  repre- 
sentative  results  presented.  The  various  techniques  for  recording  memory 
elements  are  then  presented  in  Section  3.  Representative  details  for  a  5x 
coherently  added  image  matched  filter  are  given.  Details  of  multiplexed  carriers 
are  also  given  in  the  absence  of  coverage  in  other  publications. 

Storage  media  experiences  are  described  in  Section  4.  A  detailed  description 
of  one  such  experience  comprises  Appendix  A. 

An  optical  memory  capacity  formulation  is  developed  in  Section  5  based  upon 
the  definition  of  a  memory  element  and  the  various  techniques  investigated.  The 
limitations  of  media  resolution  upon  storage  capacity  are  also  given. 

In  Section  6,  an  extensive  summary  of  the  results  are  presented  and  all 
references  given  in  Section  7.  Sections  8  and  9  summarize  the  professional 
activities  conducted  during  the  program,  and  the  professional  people  involved 
respectively. 

A  number  of  appendices  are  included  for  added  description,  supportive 
material,  and  essential  papers.  Publications  on  evolving  progress  during  the 
program,  and  others  well  documented  are  not  included  in  the  appendices. 

We  conclude  this  introduction  with  a  bit  of  what  we  believe  are  the  key  accom¬ 
plishments  of  this  research  program.  We  emphasize  that  while  each  individual  item 
is  important  in  itself,  the  total  program  has  made  a  significant  step  forward  in  the 
application  of  optical  correlator  technology  to  problems  of  target  detection,  and 
reconnaissance. 
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1.  Developed  techniques  for  constructing  large  capacity  memories. 

2.  Developed  techniques  for  utilizing  DCG  to  generate  large  memories. 

3.  Demonstrated  that  multiple  holographic  lens  -  matched  filter  (MHL-MF) 
arrays  can  successfully  utilize  large  capacity  storage  techniques. 

4.  Applied  large  capacity  memory  techniques  to  DCG  and  SIL  HAL,  and  found 
that  the  technology  for  photopolymers  needs  further  development.  Thus, 
amplitude  and  index-of-refraction  change  phase  matched  filters  were 
successfully  constructed  in  thin  and  thick  media. 

5.  Developed  a  digital  simulation  of  the  optical  correlator  which  was  validated 
by  both  visual  and  quantitative  comparisons  with  the  device  output.  The 
simulation  allows  the  optimization  of  device  parameters  under  controlled 
conditions  and  will  significantly  reduce  the  development  time  of  various 
applications  of  the  correlator. 

6.  Demonstrated  matched  filter  memory  techniques  with  large  signal-to-clutter 
(s/c)  ratios  using  a  variety  of  terrains. 

7.  Demonstrated  through  s/c  ratio,  a  compact  MHL-MF  array  optical  correlator. 

8.  Determined  the  parametric  sensitivities  of  large  memory  elements,  information 
usable  in  establishing  large  memory  banks. 

9.  Developed  a  formulation  for  capacity  which  shows  that  large  (>2000)  arrays 
of  MF/cm2  are  theoretically  possible. 

10.  Demonstrated  that  5  SMF  non-eoherently  added  and  15X  MF  coherently  added 
,  memory  banks  can  be  made.  These  are  not  limits  but  there  are  indications 

that  these  are  close  to  the  upper  limits. 

11.  Demonstrated  that  two  15X  coherent  MF  memory  elements  can  be  carrier  multi¬ 
plexed  to  yield  30  memory  elements  in  one  memory  position  -  a  number  which 
is  greater  than  the  memory  elements  needed  for  360°  coverage  of  a  tank  in  one 
range  bin. 

12.  Demonstrated  that  the  validated  model  can  be  applied  to  large  coherent 
addition  memory  elements  and  used  with  a  variety  of  terrains. 
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2  -  DESCRIPTION 


2.1  EXPERIMENT AL  SETUP 

The  experimental  test  setup  for  fabricating  and  evaluating  the  matched 
filters  is  shown  in  Fig.  2-1.  Using  an  acetate  based  film  version*  of  the  M-60  tank 
placed  in  the  rotating  film  gate,  appropriate  matched  filters  are  constructed  at 
the  matched  filter  plane.  The  latter  are  recorded  on  plates  held  in  a  specially 
designed  MF  mount  which  has  X-Y-Z,  and  0(azimuth)  and  (rotation)  adjustment 
freedom.  Figure  2-2  is  a  photograph  of  this  device.  The  output  correlation  plane 
was  measured  using  a  scanning  fiber  optic  probe  (50pm)  and  observed  at  times 
with  a  sensitive  television  camera  and  monitor.  The  probe  represents  18%  of  the 
tank  width  when  projected  back  to  the  image  plane,  a  resolution  greater  than  5 
lines  per  millimeter. 

Input  imagery  was  placed  in  the  rotating  liquid  gate  and  properly  oriented  as 
required.  Additional  mirrors,  like  M3,  were  placed  along  the  reference  beam  line. 
When  multiplexed  carriers  were  investigated,  the  sequence  in  recording  was 
usually  dictated  by  minimizing  mirror  adjustment. 

Except  for  the  laser,  the  arrangement  in  Fig.  2-1  was  maintained  through¬ 
out  the  program. 

2.2  SIMULATION  VALIDATION 

A  number  of  simulations  were  conducted  to  augment  the  empirical  aspects  of 
the  program.  Thus  it  was  necessary  to  validate  the  model  used  for  our  simulation. 
The  M-60  tank  was  digitized  and  the  (Fourier  transform)^  obtained  in  an  FFT 
program.  Figures  2-3  and  2-4  illustrate  the  analog  and  digital  versions  of  the 
spectra.  Plotting  the  results  along  an  axis  perpendicular  to  the  major  tank  axis, 
and  comparing  them  to  a  scan  of  an  analog  version  shows  a  remarkable  similarity. 
These  are  illustrated  in  Fig.  2-5  and  2-6.  Each  of  the  maxima  in  the  scan 


•Polyester  based  film  is  most  commonly  available  and  it  was  necessary  to  occa¬ 
sionally  use  it.  However,  it  is  more  birefringent  than  acetate  based  film  and  some 
precautions  are  necessary  in  MF  fabrication.  This  problem  is  detailed  in  the 
paper  contained  in  Appendix  B. 
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Fig,  2-1  Arrangement  for  Matched  Filter  Fabrication,  S/C  and  Correlation  Measurement* 
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perpendicular  to  the  major  axis  of  the  tank  is  designated  an  "order".  A  matched 
filter  optimized  for  visibility  in  the  center  of  that  spatial  frequency  band  is  called 
a  filter  of  that  order  in  the  analog  case.  For  the  digital  simulation,  the  cutoff 
threshold  is  set  at  the  maximum  within  that  band,  and  the  filter  designated  as  one 
of  that  particular  order.  The  validity  of  our  simulations  were  verified  further  by 
determining  a  signal-to-clutter  (S/C)  ratio.  In  this  measurement,  various  terrain 
scenes  are  played  through  various  orders  of  matched  filter  and  the  S/C  ratio 
measured.  Prior  to  this,  of  course,  was  the  need  to  digitize  the  terrain.  This 
was  done  and  the  (FT)2  along  the  major  axis  of  those  scenes  is  shown  in  Fig.  2-7 
to  the  same  scale  as  the  tank  spectrum  of  Fig.  2-4. 

Using  these  terrains t ,  the  auto  correlation  and  S/C  were  measured  for  vari¬ 
ous  orders  of  matched  filter.  The  results  are  shown  in  Table  2-1  and  are  quite 
similar  to  those  obtained  on  a  prior  program  (Ref  8).  Thus,  the  model  used  in  our 
simulations  was  initially  considered  validated  and  additional  assurance  sought. 

Use  was  also  made  of  the  simulation  in  order  to  determine  (a)  the  auto¬ 
correlation  variation  and  (b)  the  signal-to-clutter  (S/C)  ratio  for  the  third  order 
filter  as  the  matched  filter  is  cut  off  at  increasing  spatial  frequencies.  The 
results  are  shown  in  Table  2-2. 

Through  experiments,  the  physical  size  of  the  memory  element  was  estab¬ 
lished  as  the  restricting  spatial  frequency  filter  radius  which  allowed  95%  of  the 
correlation  signal  to  be  determined.  For  the  M-60  tank  this  was  found  to  be  8.25 
c/mm,  or  2  1/2  times  the  optimizing  spatial  frequency.  From  Table  2-1,  it  can  be 
seen  that  the  same  correlation  signal  occurs  at  a  slightly  lower  spatial  frequency, 
suggesting  that  a  smaller  MF  spatial  frequency  radius  is  possible  thus  leading  to  a 
larger  memory  capacity.  Although  the  simulation-experimental  difference  is  not 
large,  it  must  be  considered  that,  in  the  analog  case,  the  recording  medium 
actually  provides  a  low  frequency  cutoff  through  film  saturation  not  present  in 
the  simulation.  Thus,  when  the  overall  energy  passed  by  the  actual  MF  is 
considered,  the  95%  energy  point  will  be  at  a  higher  frequency  and  the 
experimental  and  simulation  results  will  agree. 


f  In  the  empirical  work,  the  woodland  scene  (Fig.  l-3b)  was  occasionally  the  only 
analog  common  to  the  simulations.  This  is  one  of  these  occasions. 


2-7 


/ 


i 


RELATIVE  INTENSITY 


Fig.  2-7  Power  Spectrum  of  Three  Digitized  Terrain  Scenes 
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Table  2-1  Validation  of  Simulation  Through  S/C  Ratio  Measurement 


TABLE  2-2  AUTOCORRELATION  AND  S/C  RATIO  FOR  M-60  IN  THREE  TERRAINS  AS  A  FUNCTION 
OF  THE  SPATIAL  FREQUENCY  CUTOFF  FOR  3RD  ORDER  IMF 
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The  S/C  ratios  are  seen  to  closely  parallel  those  previously  obtained  (Ref  8). 
In  the  woodland  scene,  95%  of  the  maximum  S/C  ratio  occurs  at  the  same  point  as 
the  autocorrelation  signal,  while  in  the  case  of  the  village  it  occurs  at  a  higher 
frequency.  The  roadside  is  unusual.  Ninety-five  percent  of  the  maximum  S/C 
ratio  occurs  just  about  at  the  optimizing  frequency  and  then  shows  a  decline  to  a 
lower  "saturation"  value. 

The  results  of  this  simulation  illustrate  a  valid  modeling  procedure.  A  full 
comparison  is  made  in  the  last  paper  presented  in  this  program  (Ref  30). 

A  final  check  on  the  validation  discussed  in  that  paper  is  illustrated  in  Fig. 
2-8  and  2-9.  In  these,  the  M-60  tank  is  reimaged  through  zero,  third,  and  sixth 
order  MF.  The  similarity  in  results  is  once  again  quite  good. 

2.3  MATCHED  FILTER  FORMATION  AND  PLAYBACK 

Formation  of  an  M-60  matched  filter  (MF)  using  the  nonlinear  silver  halide 
recording  characteristics  has  been  addressed  before.  See,  for  example.  Ref  9  for 
the  general  consideration  and  Ref  7  for  application  to  the  M-60  in  particular. 

Some  insight  to  the  matched  filter  can  be  obtained  by  examining  charac¬ 
teristics  in  the  matched  filter  pis  ie  along  the  axis  of  symmetry.  In  the  program, 
the  object  image  is  perpendicular  to  the  MF  plane.  The  reference  beam  is  offset 
but  in  the  MF  plane  of  symmetry  which  is  perpendicular  to  the  MF  plane  itself. 
This  uniqueness  does  not,  of  course,  exist  in  the  simulation. 

Expressions  have  been  derived  for  the  fringe  separation  d  of  the  MF  along 
this  axis  in  the  plane  of  symmetry  using  the  simple  geometrical  model  described 
above.  In  the  program  in  all  cases  the  signal  beam  angle,  Q  s,  was  always  zero, 
although  comparisons  with  an  equal  beam  case,  i.e.,  0'8  =  or,  are  considered. 

In  those  cases  we  have  the  condition  e'9=  6r  =  es/2.  The  resulting  equation  is: 


dX~l  = 

SinOD  +  Sin 

Tan'YSines  ' 

V  Cos68  +  C\ 

tv 

\  Cose 

+  V  Sine  / 

_ 

'  S 

s  >. 

with  F  being  the  Fourier  transform  lens  focal  length,  and  V  =  H/F,  and  C  =  UF 
the  relative  object  and  spatial  frequency  locations,  respectively.  H  =  0  =  K 
represents  an  on-axis  object  at  zero  spatial  frequency.  Figure  2-10  illustrates 
the  results  of  fringe  spacing  for  the  object  on-axis  case.  The  results,  up  to  a 
spatial  frequency  of  ±15c/mm,  are  most  meaningful  in  terms  of  the  memory  element 
size  definition  given  earlier.  It  is  clear  that  the  matched  filter,  as  constructed. 
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Fig.  2-6  Raimaging  of  M-60  Match  ad  Fitters  (Analog) 


Fig.  2-9  Reimaging  of  M-60  Matched  Filter*  (Digital) 
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is  a  chirp  grating  similar  to  that  discussed  by  Livanos,  et  al.  (Ref  10).  When  the 
object  position  is  off-axis,  although  similar  results  are  obtained,  the  =  0  cross¬ 
over  point  shifts.  These  results  are  more  evident  in  Fig.  2-11.  The  limits  on  V 
were  derived  from  the  assumption  of  a  35  mm  format,  F  =  50  mm.  Looking  at  the 
results  we  find  asymmetry.  At  0  r  =  5°  and  0  §  =  0,  the  ratio  d/  X  is  greater  than 
the  corresponding  curve  for  0  s  =  0  R  =  2.5°  by  less  than  1%.  However,  as  the 
reference  beam  axis  angle  increases,  say  to  60°,  the  variation  in  d/  X  is  propor¬ 
tionately  larger  than  the  corresponding  one  for  0  s  =  0R  =  30°. 

The  results  of  Fig.  2-10  are,  of  course,  independent  of  the  transform  lens 
focal  length  for  the  correlator  system  parameter.  S  =  (XF)-1  accounts  for  the  size 
but  does  not  alter  the  spatial  frequency  characteristics  of  the  pattern. 


One  extension  of  Eq  (2-1)  is  the  application  to  multiple  holographic  lenses  in 
correlators  with  large  memory.  In  this  case  the  equation  still  holds  except  that 


6  =  Tan 
s 


ND  ±  a) 


(2-2) 


where  N  is  the  integral  number  of  offsets  of  memory  position  in  the  corresponding 
matched  filter  plane,  and  D  the  MF  separation. 

The  fringe  spacing  equation  holds  with  this  adaptation.  Note  that  the 
"chirping"  becomes  different  for  each  matched  filter  across  such  an  array. 

The  results  of  Eq  (2-1)  can  also  be  used  to  establish  the  total  number  of 
matched  filters  between  any  two  reference  beam  angles  for  the  case  in  which  no 
two  filters  have  the  same  carrier  frequency.  This  would  be  illustrated  in  Fig. 
2-10  if  the  upper  transform  spatial  frequency  limit  of  one  filter  would  be  adjacent 
to  the  lower  limit  for  the  adjacent  filter  at  some  other  reference  beam  angle.  For 
example,  if  the  minimum  (physically)  realizable  reference  beam  angle  were  5°,  the 
adjacent  reference  beam  angle  would  be  the  one  whose  lower  transform  spatial 
frequency  starts  at  a  normalized  carrier  frequency  of  9.4  (  v  =  30  c/mm)  and 
decreases.  Examination  of  Fig.  2-10  indicates  that  at  least  eight  such  filters  are 
possible.  For  a  reduced  transform  bandwidth  more  are  possible,  and,  in  combina¬ 
tion  with  changes  in  r,  a  large  array  of  MF  with  no  common  carriers  is  possible. 

On  playback,  the  object  may  occur  anywhere  in  the  input  plane  including,  of 
course,  its  position  of  construction.  The  former  is,  in  fact,  the  more  probable 
case. 
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Fig.  2-1 1  Fringe  Spacing  at  v  *  0  for  Different  Object  Position  During  IMF  Fabrication  (<?s  -  01 


The  MF  playback  equation  is  given  by 

-l/Sin  0  -  V0  Cos  +  c 


±  Sin  R  =  Sin  oD  +  Sin 
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Cos  o  +  VD  Sin  e 
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1/  Sin  0  -  V  Cos  o  +  c 
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\  Cos  6  +  V  Sin  9 
'  s  s 


(2-3) 


where  VQ  is  HQ/F,  the  relative  object  position  during  MF  construction  and  V,  the 
relative  position  during  playback.  The  sign  for  (Sin  R)  depends  upon  the  relative 
values  for  H  and  c  .  Obviously,  playback  occurs  at  0  r  when  construction/ 
playback  are  the  same. 

Equation  (2-3)  can  be  used  to  determine  the  output  position  of  the  object  in 
the  correlation  plane  (along  the  axis  defined  before)  for  any  position  in  the  object 
field.  When  H  =  H0,  the  last  two  right-hand  terms  are  not  equal.  The  diffraction 
angle  R  differs  from  9  r.  If  we  call  the  difference  (9r  -R)  the  deviation,  i.e., 

6  =  ( 0 r  -R),  the  input-correlation  plane  movement  is  related  by: 


y 

— ■  =  s  =  Tan  6  =  Tan  (0D  -  R)  (2-4) 

Ft 

where  Ft  is  the  inverse  transform  lens  focal  length.  If  Ft  is  equal  to  the 
transform  lens  focal  length,  we  can  get  a  normalized  deviation  as  the  limits  of  a  35 
mm  film  format  is  traversed  (i.e.  +42  mm). 

Equation  (2-4)  has  been  computed  for  the  usual  two  cases: 

0$  =  0,  ®R=e‘*  ®R  =  0/2.  The  deviation  increases  with  increasing  total  angle  but 
with  decreasing  input  field  coverage.  Symmetry  is  best  for  the  smallest  angle. 


When  we  consider  a  holographic  lens,  MF  (square)  array  in  a  memory  system, 
the  center  of  any  MF  is  located  Nd+  from  the  optic  axis,  where  N  =  0,  1,  2,  ...n. 
Thus, 


0  = 

s 


(2-5) 


Using  this  equation  in  the  previous  equations  extends  them  to  the  one-dimensional 
"array"  of  matched  filters. 


+  Remembering  the  model  covers  only  the  axis  in  the  correlator  plane  of  symmetry 
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3  -  TECHNIQUES  FOR  RECORDING  MEMORY  ELEMENTS 


There  are  a  number  of  approaches  to  multiple  storage  for  matched  filters 
which  might  be  considered  to  increase  memory  density,  with  the  choice  dependent 
upon  the  exact  correlator  configuration,  memory  requirement,  and  specific  applica¬ 
tion.  The  techniques  are: 

a.  Geometrical  configurations 

b.  Split  Spectrum 

c.  Matched  filter  additions 

Non-coherent 

Coherent 

d.  Multiplexed  carriers 

Figure  3-1  illustrates  the  various  techniques. 

3.1  GEOMETRICAL  CONFIGURATIONS 

If  a  memory  element  is  considered  to  be  as  defined  in  Fig.  1.2,  then  one  can 
establish  a  regular  matrix  of  elements  whose  upper  frequency  cutoff  limits  are  con¬ 
tiguous.  When  the  array  is  of  the  same  target  (and  matched  filter),  the  elements 
can  be  offset  with  an  increase  in  the  packing  density  of  (2/^/3).  Obviously,  the 
combination  of  several  objects  whose  upper  cutoff  frequencies  differ  will  lead  to  a 
wide  variety  of  packing  and  require  an  individual  analysis  to  optimize  the  effi¬ 
ciency  in  packing. 

Figure  3-2  shows  the  two  packing  options,  regular  and  offset,  although  for 
the  former,  they  are  not  contiguous.  By  way  of  illustration,  a  five-coherently- 
added-image  matched  filter  has  been  recorded  in  each  of  the  16  MF  positions  (by 
using  a  4x4  multiple  holographic  lens).  Thus,  there  are  basically  5x16  =  80 
memory  positions  for  the  M-60  tank.  In  the  figure  they  are  all  the  same,  when,  in 
fact,  each  would  be  modified  by  an  orientation  rotation  or  scale  change  to  provide 
360°  coverage  of  the  target.  Past  results  have  shown  that  only  23  whole  filters 
are  required  to  cover  the  M-60  with  a  £  19%  range  bin.  Thus,  for  the  memory  bank 
illustrated  in  Fig.  3-2(a),  3  1/2  range  bins  could  be  covered. 
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Fig.  3-1  Techniques  for  Recording  Memory  Elements 
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(A)  REGULAR  ARRAY 
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(B)  OFFSET  ARRAY 


Fifl.  3-2  Geometrical  Placement  of  Memory  Element* 
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In  Figure  3-2b  a  reimaging  of  the  AF  Tribar  target  is  used  to  illustrate  offset 
packing.  The  target  was  used  to  evaluate  the  resolution  of  several  multiple 
holographic  lenses  prepared  in  a  new  technique. 

3.2  SPLIT  SPECTRUM 

From  the  Fourier  transform  (FT)  equation  of  any  two-dimensional  distri¬ 
bution,  it  can  be  seen  that  any  axis  through  F(£=  0,  n  =  0)  is  an  axis  of  symmetry. 
One  half  of  the  FT  hologram  or  matched  filter  then,  should,  in  principle,  provide 
sufficient  information  about  the  object  image  to  generate  autocorrelation  signals 
and  acceptable  S/C  ratios.  In  practice,  as  used  in  the  analog  system  in  this  inves¬ 
tigation,  the  axis  of  symmetry  is  not  arbitrary  because  of  the  offset,  but  is 
uniquely  determined  by  the  signal  and  reference  beam  vectors.  The  reason  for 
this  uniqueness  is  the  reference  beam  offset  which  gives  rise  to  an  asymmetric 
chirp  in  the  matched  filter  which,  in  itself,  is  symmetric  above  and  below  the 
aforementioned  plane. 

The  practical  difficulties  of  implementing  a  split  spectrum  experiment  far  out¬ 
weighed  the  anticipated  benefits,  so  resort  to  our  verified  simulation  was  made. 

The  digital  analysis  of  the  split  spectrum  evaluation  was  performed  by 
zeroing  the  redundant  portion  of  the  two-dimensional  Fourier  transform  of  the 
tank,  the  object  image.  The  resulting  half-transform  was  then  used  to  construct 
the  digital  MF  and  various  high  pass  versions  of  the  MF.  The  woodland  scene 
(Fig.  l-3b),  previously  digitized,  was  used  to  play  through  the  filter. 

The  results  of  the  split  spectrum  test  show  that  the  autocorrelation  signal  is 
halved  and  the  S/C  ratio  somewhat  less  than  this  as  shown  in  Table  3-1. 

Clearly,  the  autocorrelation  is  as  expected,  but  the  cross  correlation,  and 
the  S/C  ratio  derived  from  it,  are  not.  The  deviation  in  the  l“.st  two  are  computa¬ 
tional  artifacts  of  the  simulation  attributed  to  the  inclusion  of  the  imaginary  terms 
which  do  not  integrate  out  when  the  split  spectrum  is  used. 

In  a  second  simulation,  four  tanks  were  randomly  though  squarely  placed  in 
the  woodland  scene.  Then  the  correlation  plane  was  searched  and  the  five  largest 
signals  observed  as  being  either  a  tank  or  a  false  alarm.  In  both  the  cases,  i.e., 
full  and  half  spectrum  matched  filters,  all  four  tanks  were  picked  out  first.  This 
is  summarized  in  Table  3-2. 
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FULL 

SPECTRUM 

MF 


HALF 

SPECTRUM 

MF 


TABLE  3-1  SPLIT  SI 


SIMULATION 

SUPPRESSION 

LEVEL 

1.394 

.026 

.034 

.034 

.010 

.001 

1.394 

.026 

.034 

.034 

.010 

.001 


X283-Q24(a)(T) 


SIMULATION  RESULTS 


NORMALIZED 

CORRELATION 

CAUTO 

CROSS 

1.000 

1.571 

■ 

.068 

.100 

.too 

.100 

.038 

.023 

.011 

1 

.500 

1.110 

1 

.091 

.037 

.105 

.054 

.105 

.054 

.050 

.021 

. 

.012 

.005 
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TABLE  3-2  ACCURACY  OF  CORRELATION  LOCATION  WITH  THE  SPLIT-SPECTRUM 


MAXIMUM 

VALUE 

LOCATION 

TARGET/ 

NO  TARGET 

MISS 

DISTANCE 

ROW 

COL. 

4.21 

44 

140 

T 

1  PIXEL 

FULL 

4.15 

359 

111 

T 

0 

SPECTRUM 

MF 

4.00 

439 

341 

T 

0 

3.57 

139 

291 

T 

0 

1.81 

32 

187 

NT 

- 

2.11 

T 

1  PIXEL 

HALF 

2.08 

T 

0 

SPECTRUM 

2.00 

(SAME) 

T 

0 

MF 

1.78 

T 

0 

0.90 

NT 

- 

1283-01 7<a)(T) 


It  can  be  seen  that  the  half  spectrum  matched  filter  (a)  provides  locational 
accuracy  comparable  to  the  full  spectrum  MF  (with  the  specified  terrain)  and  (b) 
enables  the  four  targets  to  be  determined  prior  to  the  first  false  alarm. 

As  a  memory  element,  the  half  spectrum  performs  satisfactorily. 

3.3  MATCHED  FILTER  ADDITIONS 


Matched  filters  can  be  combined  in  a  noncoherent  or  coherent  manner 
although  one  must  be  careful  to  distinguish  the  meaning  of  the  addition.  Matched 
filters  can  be  sequentially  overlayed  at  the  same  position  on  the  recording  medium 
and  we  can  properly  refer  to  this  as  the  noncoherent  addition  of  matched  filters. 
Alternatively,  one  can  make  a  single  MF  as  the  result  of  a  coherent  addition  of 
images  in  a  single  exposure.  This  is  referred  to  as  the  coherent  addition  MF,  or 
coherent  MF.  The  distinction  between  the  two  should  be  noted  although  reference 
to  them  will  occasionally  imply  that  MFs  are  coherently  added.  Figure  3-3 
illustrates  cases  for  2x  and  3x  coherent  and  noncoherent  addition.  In  the 
noncoherent  case,  rotation  of  the  single  tank  power  spectrum  is  evident. 


Examination  of  the  noncoherent  and  coherent  addition  of  "n"  matched  filters 
can  be  done  using  a  linear  recording  response  (i.e. ,  T^  =  T0  -  Et  =  T0  -  3  It)  and 
obtaining  for  the  amplitude  transmittance: 


T.  =  -Bt  [R‘  -  T  /6t]  -  Bt  2  (S?+  R*S.  +  RS*  )  -  fit  2  (S.  S *) 

A  o  i=l  1  1  1  i  .3=1  1  ] 


(3-1) 


for  the  coherent  case  where  Rj,  Sj  are  the  reference,  signal  amplitudes;  T0,  3 
film  characteristics;  and  t  exposure  time.  The  noncoherent  result  is  the  same 
equation  without  the  coupled  terms  at  the  right  hand  side.  For  simplicity,  it  was 

assumed  in  the  latter  case  T^  =  Tj  +  T2  + .  Lang,  et  al.,  (Ref  22) 

developed  more  extensive  relationships  for  the  two  cases  and  an  explicit  form  for 
the  interobject  interference  terms  contained  in  the  coupled  term  of  Eq  (3-1).  The 
phenomenon  arising  from  this  interference  is  illustrated  in  Fig.  3-4.  The  latter 
shows  the  spatial  frequency  optimization  region  (called  the  third  lobe  or  order) 
for  MF  of:  (a)  a  single  (lx),  (b)  a  four  symmetric  tank  case  (4x),  and  (c)  an 
eight  symmetric  tank  case  (8x).  The  intraobject  interference  effects  can  be 
clearly  seen  in  the  last  two.  Measurement  shows  the  spacing  (g)  dependence  is 
predictable  from 


g  =  X  (sin  ^) 


(3-2) 
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where  Sin  ^  =  (h  F*-1)  ((hF"l)2  +  l)“l/2  and  ^  and  x  are  the  transform  lens 
focal  length,  image  spacing  and  wavelength,  respectively.  The  carrier  inter¬ 
ference  is  clearly  visible.  Lang,  et  al.,  have  shown  that  (image)  interference 
varies  according  to  Cos  (a^  -  a  ,)  ,  where  the  terms  refer  to  the  phases  of  the 
objects  involved.  Thus,  when  (a  -  a  ()  =  n/2  the  deleterious  interference 
effect  disappears  and  the  correlation  ought  to  assume  its  maximum  value  possible. 
This  occurs  at  some  optimum  spacing.  As  more  than  two  objects  are  involved,  the 
number  of  phase  terms  to  be  considered  rises  rapidly. 

A  question  arises  as  to  whether  these  interobject  interferences  reduce  the 
autocorrelation  or  not,  and  whether  or  not  a  random  placement,  or  at  least  a 
uniqueness  in  the  geometrical  placement,  might  be  advantageous.  Matched  filters 
were  derived  from  a  set  of  images  consisting  of  lx,  4x,  and  8x  tank  placements, 
and  the  autocorrelations  were  obtained  with  the  results  shown  (Table  3-3). 

It  seems  that  the  regular  placement  does  not  detrimentally  impact  the 
autocorrelation.  In  a  second  test,  an  8x  MF  was  constructed  from  a  lOx  image 
object  by  using  the  first  eight  tanks.  In  a  similar  manner,  the  autocorrelations 
were  measured.  The  results  were  (coincidentally)  the  same,  i.e.,  cn  =  .120, 
when  compared  to  a  lx  autocorrelation.  Both  results  were  obtained  with  indepen¬ 
dent  tests. 

A  question  also  arises  as  to  whether  such  intersymbol  interference  arises  in 
the  correlation  plane  since  an  object  is  often  regularly  replicated  in  an  array  when 
there  is  a  multiple  holographic  lens-MF  combination.  An  experiment  was  per¬ 
formed  in  which  the  center  position  of  a  5  x  5  holographic  lens  -  MF  array  -  was 
used  to  obtain  an  autocorrelation  signal.  This  is  shown  in  Fig.  3-5a.  Then,  the 
autocorrelation  signal  for  the  center  column  was  obtained  as  shown  in  Fig.  3-5b. 
The  spacing  between  the  intraimage  interference  fringe  can  be  computed  from  Eq 
(3-2)  by  substituting  h  =  NhOI  N  =  1,2, 3,4.  Finally,  when  all  25  positions  are 
allowed  to  transform  and  the  autocorrelation  obtained,  the  signal  in  Fig.  3-5c 
results.  Clearly  this  is  a  worse  case  situation  in  which  all  MF  are  identical  so 
that,  except  for  spacing  effects,  interference  between  all  elements  is  greatest. 

The  situation  is  similar  to  a  multiple  Young  experiment. 

As  with  all  comparisons,  the  autocorrelation  signal  and  S/C  ratio  were 
measured  with  the  SE  Asian  scene.  The  results  are  summarized  in  Table  3-4.  A 
scan  across  the  autocorrelation  peak  of  Fig.  3-5c  is  given  in  Fig.  3-6.  There 
is  a  reduction  in  autocorrelation  signal. 
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TABLE  3-3  AUTOCORRELATIONS  FROM  REGULARLY  PLACED  IX.  4X  AND  8X 

INPUT  OBJECTS 


OBJECT 

SIZE 

K 

OBTAINED  FROM 
ANGULAR 
POSITION 

NORMALIZED  PEAK 
AUTOCORRELATION 

*n 

AVERAGED 

IX 

0° 

1.000 

1.000 

0° 

.325  1 

4X 

90 

.425  f 

.388 

180 

.470  | 

270 

.333  ) 

0° 

.096 

45 

.101 

90 

.129 

8X 

135 

.136 

.  .120 

180 

.144 

225 

.151 

270 

.118 

315 

.083  t 

1283-018<»)(T) 
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AUTOCORRELATION,  c, 

.60 


S/C  RATIO  (dBI 


MF 

AUTOCORREL 

MF33 

.60 

COLUMN  3 

.65 

ALL  25 

1.00 

1283-019(a)(T) 

4.38 

5.16 

5.00 


CORRELATION  INTENSITY 


SCAN  ALONG  FIELD  - ► 


Fig.  3-6  Scan  Along  Major  Tank  Axil  in  Correlation  Plana  Using  5x6  HL  Transform  Lara 
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3.3.1  Noncoherent  Addition 


Noncoherent  addition  is  a  literal  overlap  or  addition  of  matched  filters  within 
the  limits  of  a  memory  element  defined  earlier.  The  MF  were  optimized  in  terms  of 
fringe  visibility  at  the  desired  spatial  frequency  through  the  R  =  (Reference 
Beam/Signal  Beam)  ratio,  and  the  equal  exposure  requirements  of  the  amplitude 
transmittance-exposure  (t^  -  E)  relationship  for  silver  halide,  or  efficiency- 
exposure  (n~  E)  for  dichromated  gelatin.  Figure  3-7  illustrates  a  few  basic 
configurations  of  noncoherent  addition  including  the  use  of  orthogonal  carrier 
beams  (see  Ref  11).  This  variant  is  not  covered  here  but  does  influence  the 
overall  memory  capacity.  Figure  3-7  shows  the  cases  where  (a)  the  image  is 
rotated,  (b)  the  MF  are  overlapped  at  the  same  carrier  frequency,  and  (c)  the 
case  where  the  MF  recording  medium  is  rotated.  Only  the  first  scheme  (a)  was 
used  in  this  program.  The  overlapped  MF  (b,e)  are,  of  course,  redundant  except 
when  the  carrier  frequency  is  different,  and  this  is  treated  later.  In  the  last 
case,  we  are  aware  of  only  Almeida  (Ref  12)  using  rotated  MF  planes.  It  must  be 
realized,  of  course,  that  playback  for  the  c,f  techniques  generates  independent 
correlation  planes. 

An  interesting  variation  of  noncoherent  addition  occurs  when  a  coherent 
array  of  objects  (tanks)  is  used  individually  to  make  a  noncoherent  array.  This 
is  discussed  below. 

A  series  of  filters  was  made  up  to  5x  and  their  effectiveness,  as  measured  by 
the  S/C  ratio,  determined.  In  addition,  the  rotation  sensitivity  was  also  deter¬ 
mined  to  ensure  an  experimental  similarity  in  bandpass  characteristics.  These 
filters  were  constructed  in  dichromated  gelatin.  The  results  are  summarized  in 
Table  3-5.  Pictures  of  the  2x  and  5x  cases  in  dichromated  gelatin  are  shown  in 
Fig.  3-8,  while  the  impulse  response  of  the  latter  case  is  shown  in  Fig.  3-9a. 
Clearly,  rotation  of  the  image  was  done  with  an  osculating  center.  Using  the  5x 
matched  filter,  the  SE  Asian  scene  was  played  through  and  the  S/C  measured  to  be 
8.87  dB.  The  results  of  a  scan  across  the  tank  in  this  evaluation  is  shown  in  Fig. 
3-10. 

3.3.2  Coherent  Addition 

Coherent  addition  entails  the  preparation  of  an  object  transparency  which 
consists  of  an  array  of  targets,  or  some  target  with  an  array  of  views  or  orienta¬ 
tions.  Then,  a  single  matched  filter  is  made  of  the  array.  As  seen  from  Eq  3-1, 
the  MF  is  more  complex  than  the  addition  of  several  individual  matched  filters. 


/ 
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TABLE  3-5  PERFORMANCE  OF  NONCOHERENT  Y  ADDED  MF  IN  DICHROMATED  GELATIN 


#  OF 

MATCHED 

FILTERS 

AT  POSITION 

THICK 

THIN 

s'c 

db 

Cn 

AVE. 

A0 

DEGREES 

S/C 

db 

Cn 

AVE. 

A# 

DEGREES 

t 

6.87 

i.oo 

t  6.75 

5.00 

1.00 

t  7.6 

2 

8.46 

.482 

*  6.99 

5.60 

.33 

t  6.6 

3 

5.67 

.678 

t  6.70 

5.40 

.97 

i  6.5 

4 

8.81 

.449 

±  6.47 

6.80 

.72 

±  6.6 

5 

6.29 

.968 

t  5.35 

- 

- 

- 

1 283-020(1)  (T) 
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IMAGE  ROTATED  SECOND  MF  SUPERIMPOSED  PLATE  ROTATED 

(D)  (E)  (F) 
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Fig.  3-7  Non- Coherent  Addition  et  Constant  Carrier  Frequency 
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Fig-  3-9  Impulse  Response  of  5x  Image  Matched  Filters  in  Dichromated  Gelatin 


Consider  that  the  targets  are  arrayed  as  in  Fig.  3-11  with  either  (a)  a 
regular  or  (b)  random  array.  Let  E  (x,  y)  be  the  transmittance  of  a  single  target. 
Then,  the  effect  of  rotation  and  translation  on  i-targets  yields  the  Fourier 
transform  (Appendix  D)  of  the  array: 


F(2uvx,  2itVy)  = 


V  'i(w  x.  +  w  y.)  fr  -jx’(o>  Cos 0.  +  w  Sine.) 

he  x  1  y  1  JJ  E(x\  y')e  x  y 

-jy'(-u>  Cos9.  +  u;  Cose.) 
x  l  y  l 

a  J 


(3-3) 


where  the  x' ,  y'  are  the  target  coordinates,  and  the  0[  the  corresponding  angles 
for  i-th  target: 


x'  =  (x-xp  Cos  e.  +  (y-y.)  Sin  e. 
y'  =  -(x-x^  Sin  e.  +  (y-y^  Cos  ei 


(3-4) 


The  complexity  of  the  spectra  of  the  coherent  arrays  suggests  that  some 
organized  method  of  optimization  was  needed  other  than  trying  to  maximize  fringe 
visibility  at  a  particular  spatial  frequency.  While  the  latter  technique  works  well 
for  the  well  defined  single  target  case,  it  does  not  appear  applicable  to  the  multi¬ 
ple  case.  Rather,  the  situation  was  handled  by  using  the  lx  target  for  optimiza¬ 
tion  and  then  making  appropriate  adjustments  in  the  beam  ratio  R  for  targets  of 
increasing  complexity.  This  proved  effective  for  silver  halide  and  dichromated 
gelatin  and  was  consistent  with  the  work  of  Mehta  (Ref  13)  which  called  for  a 
modified  fringe  visibility 

2R*N 

v  = - -  (3-5) 

R  +  In 


where  N  is  the  number  of  images  used. 

Images  from  lx  to  lOx,  in  addition  to  15x,  were  layed  out  and  photographed, 
and  the  negative  used  in  the  subsequent  tests.  Individual  autocorrelations  were 
determined  for  the  lx-lOx  set  and  are  shown  in  Table  3-6.  The  table  brings  out  a 
limitation  in  using  analog  methods  in  the  preparation  of  the  imagery.  Despite  the 
exercise  of  care,  unequal  images  still  result.  While  the  N-averages  show  a 
general  decline,  and  while  a  few  of  the  multiple  images  show  a  fair  degree  of 
equality,  a  number  are  better  than  2:1  in  high/low  ratio. 
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TABLE  3-6  POSITIONAL  DISTRIBUTION  OP  AUTOCORRELATIONS  OF  COHERENT  MATCHED 
FILTERS  USING  THIN  SILVER  HALIDE 


CORRELATION  OF  COHERENT  IMAGERY  | 

IX 

2X 

3X 

4X 

5X 

6X 

7X 

8X 

9X 

10X 

1.000 

.492 

.269 

.266 

ESI 

.058 

.054 

.043 

.038 

.442 

.222 

.160 

.043 

.047 

.029 

.024 

.213 

.167 

ESI 

.083 

.046 

.045 

.030 

.017 

.254 

.298 

.074 

.039 

.027 

.022 

.024 

.271 

.071 

.049 

.032 

.025 

.023 

.071 

.044 

.021 

.025 

.020 

.049 

.032 

.025 

.024 

.047 

.030 

.035 

.014 

JE|M 

N-AVERAGE 

1.000 

.467 

.235 

.212 

.291 

.070 

.048 

.038 

.028 

.025 

HIGH/LOW  RATIO 

1.000 

1.113 

1.263 

1.663 

1.259 

6.867 

1.487 

2.571 

3.071 

2.235 

1283-02l(a)(T) 
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Throughout  the  program,  specific  test  objectives  were  met  by  sets  made 
specifically  for  it  with  subsequent  differences  in  correlation  values.  This  pro¬ 
cedure  was  used  to  isolate  a  test  and  account  for  differences  in  emulsion  batch, 
environmental  factors,  and  sometimes  personnel.  The  5X  was  selected  as  repre¬ 
sentative  of  the  group  as  being  the  largest  size  with  the  most  uniformity,  and 
thus,  subject  to  a  more  comprehensive  set  of  evaluations.  The  set  is  often  com¬ 
pared  to  correlation  results  from  a  5X  noncoherent  addition  set. 

3.3.3  5X  COH  and  Non-COH  Comparisons 

The  coherent  and  noncoherent  5X  matched  filter  performances  were  compared 
on  the  basis  of  individual  autocorrelation  signals,  S/C  ratio,  and  angular  and 
translational  sensitivity.  After  the  input  imagery  was  set  at  the  proper  angle  and 
masked,  the  autocorrelation  signals  for  each  of  the  five  positions  was  determined. 
At  any  one  time,  some  degree  of  response  was  obtained  at  the  other  four  non- 
aligned  positions,  although  at  significantly  reduced  levels.  The  views  of  the 
correlation  plane  at  two  threshold  levels  are  shown  in  Fig.  3-12,  3-13  for  both 
techniques  at  comparable  16:1  exposure  levels.  The  same,  but  otherwise 
arbitrary,  angle  was  chosen  in  each  case.  The  numerical  results  are  shown  in 
Table  3-7  for  the  two  silver  halide  emulsions.  Figure  3-14  shows  the  full  5X  MF 
autocorrelation  under  the  same  conditions  as  those  of  Fig.  3-12  and  3-13.  The 
element  placements  should  be  noted.  The  impulse  response  of  the  5X  COH  array 
was  previously  shown  in  Fig.  3-9. 

Each  of  the  filters  was  tested  against  the  SE  Asian  scene  with  5  vehicles,  and 

the  S/C  ratio  determined  by  scanning  through  the  autocorrelation  signal.  The 

clutter  obtained  is  that  from  the  same  scan  line.  These  results  gave  S/C  ratios 

somewhat  larger  than  might  be  obtained  if  the  whole  field  were  searched  for  the 

peak  clutter.  However,  for  direct  comparisons  the  "standard"  scans  were 

employed.  Other  tests  have  shown  that  an  all-plane  scan  would  yield  somewhat 

lower  ratios  for  this  terrain  scene  but  not  always  because  of  the  same  terrain 

feature  for  the  different  MF.  The  results  of  these  scans  for  each  matched  filter, 

* 

and  for  each  of  the  72°  spaced  images,  are  given  in  Table  3-8.  It  can  be  seen  that 
the  S/C  ratio  for  5X  COH  and  Non-COH  matched  filters  all  show  better  than  4:1 
ratio,  while  the  thick  MF  is  about  1.5  dB  lower  than  the  thin  MF. 

Figure  3-15  illustrates  the  correlation  plane  for  the  5X  Non-COH  matched 
filter  when  the  SE  Asian  scene  with  5  vehicles  is  used.  The  threshold  ratio 
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(a)  AUTOCORRELATION.  C  (b)  C  AT  REDUCED  THRESHOLD 

Fig.  3-12  Autocorrelation  of  Single  M-60  Tank  Playing  Through  5X  Non-Coherent  Matched  Filter 
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I?'  CORRELATION,  C 


(b)  C  AT  REDUCED  THRESHOLD 

1 283-81-028(T) 


Fig.  3-13  Autocorrelation  of  a  Single  M-60  Tank  Playing  Through  5X  Coherent  Matched  FHter 
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(al  NON-COHERENT  ADDITION 


*•  •, * 


(b)  COHERENT  ADDITION 

1283-81-029(T) 

Fig.  3-14  Autocorrelation  Planet  When  Entire  5x  Image  Playt  Through  Non-Co herently 
and  Coherently  Made  Matched  Filters 
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(a)  CORRELATION  PLANE  FOR  TERRAIN  IMAGE  PLAYING 
THROUGH  5X  NONCOHERENTLY  ADDED  MF 


(bl  AT  REDUCTION  TO  20%  EXPOSURE 
1283-ai-030(T) 


Fig.  3-15  Correlation  Plane  of  a  5x  Non-Coherently  Added  Matched  Filter 
Using  SE  Asian  Terrain  Image  Input  Upon  Playback 
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-S'  CORRELATION 
n 

NONCOHERENT 

COHERENT  | 

IMAGE  (°) 

THIN 

THIN 

THICK 

•  144 

0.790 

0.793 

0.731 

-  72 

0.816 

1.000 

1.000 

12  O'CLOCK 

0.803 

0.998 

0.809 

+  72 

0.943 

0.915 

0.657 

+  144 

1283-022<a)(T) 

1.000 

0.825 

0.540 
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TABLE  3-8  S/C  RATIO  FOR  COM  AND  NON-COH  5X  MATCHED  FILTERS 


S/C  RATIO  (dBl 

NONCOHERENT 

COHERENT  j 

IMAGE  IDEG) 

THIN 

THIN 

THICK 

- 144 

8.76 

9.42 

4.56 

-  72 

5.26 

8.51 

8.51 

12  O'CLOCK 

8.87 

8.82 

5.28 

+  72 

6.86 

6.88 

5.71 

+  144 

7.57 

6.30 

6.86 

AVERAGE 

i283-023(a)(T) 

7.45 

7.99 

6.16 

0 
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between  the  two  figures  is  5/1.  Further  reductions  reveal  the  solitary  prattli¬ 
n'*  oe  of  the  autocorrelation  signal.  The  pentagonal  nature  arises  from  the 
aperture  defining  the  input  field,  although  full  illumination  of  the  35  mm  format 
leaves  the  S/C  ratio  unaltered.  Tests  with  the  M-60  near  the  aperture  yield 
comparable  S/C  results  and  only  when  the  tank  becomes  obscured  does  the 
autocorrelation  signal  decrease. 

The  translational  sensitivity  of  the  autocorrelation  signal  was  determined  for 
the  5x  MF  and  compared  to  that  previously  done  (Ref  7)  for  a  single  tank  filter 
(lx).  Single  images  were  played  through  the  three  complex  matched  filters  and 
the  matched  filters  translated  in  the  MF  plane.  In  each  case,  the  filters  were 
tuned  for  autocorrelation  peaks  at  the  outset  of  each  test. 

Figure  3-16  shows  the  autocorrelation  signals  as  a  function  of  translational 
movement  parallel  and  perpendicular  to  the  major  tank  axis,  respectively. 

To  a  close  approximation,  a  single  (1  x  2)  mm  tank  image  has  a  Fourier 
transform  whose  zeros  are  at  the  zeros  of  Fy  =  Sin  y/y  and  Fx  =  Sin  x/x  for  the 
parallel  and  perpendicular  directions  respectively.  Putting  the  wavelength,  focal 
length,  and  image  dimensions  in  expanded  forms  of  these  equations  yields  zero 
intervals  of  114  and  228  u  m,  respectively.  Based  upon  a  simple  overlap  model,  it 
can  be  seen  that  axial  translations  of  57  and  117  ysm  should  yield  50%  frequency 
domain  overlaps,  and  thus  perhaps  a  -3  dB  change  in  correlation.  From  the  data, 
we  can  conclude  that  multiple  image  matched  filters  differ  little  from  the  lx  MF  in 
translation  sensitivity  (Ref  7),  or  from  a  simple  rectangle  model.  The  orthogonal 
translations  are  not  the  same,  but  conform  to  the  results  of  the  rectangular  model. 
This  is  conclusively  seen  in  the  results  for  the  lOx  coherently  added  image  MF 
which,  in  a  single  test,  was  directly  compared  with  a  lx  MF  (Fig.  3-17,  3-18). 


Past  experience  has  shown  the  axial  movement  determination  unnecessary. 
Consider  the  depth  of  focus  equation 


2<j> 

A  =  - (p 

1  +  4> 


(3-6) 


where  <f>  =  F  0r-1ro”2)  wjth  rQ  being  the  input  beam  radius.  For  the  parameters  of 
our  system,  and  allowing  a  10%  variation  in  "focal  spot  size”  (i.e.,  p  =  1.10),  the 
depth  of  focus  becomes  9.54  mm,  yielding  little  impact  upon  the  autocorrelation. 
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TRANSLATION,  Mm 

MF  BELOW  TRANSFORM  MF  ABOVE  TRANSFORM 


ALONG  MAJOR  TANK  AXIS 


PERPENDICULAR  TO  MAJOR  TANK  AXIS 

Fig.  3-16  5X  MF  Trantlttional  Santitivity 
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100  50  TRANSLATION,  nm  50 


1283-81-033(T) 

Fig.  3-18  Y- Positional  Sensitivity  of  IX  and  10X  Matched  Filters  in  Thick  and  Thin  Emulsions 


Since  the  size  of  the  matched  filter  depends  upon  (AF)~1,  it  is  expected  that, 
with  the  simple  model  cited,  translational  sensitivity  will  vary  in  a  similar  manner 
according  to  the  "size  sensitivity." 

Azimuthal  sensitivity  is  determined  by  the  rotation  of  the  matched  filter  or 
gratings  about  an  axis  lying  in  the  emulsion-substrate  plane  and  perpendicular  to 
the  reference-signal  beam  plane,  and  measurement  of  the  first  order  signal  or 
autocorrelation.  It  enables  one  to  determine  (a)  alignment  criticality,  (b)  the 
potential  for  carrier  frequency  modulation  as  a  mode  of  increasing  memory 
capacity,  and  Cc)  diffraction  efficiency  of  the  gratings. 

In  making  these  assessments,  the  matched  filters  were  made  and  the  auto¬ 
correlation  signals  measured  as  the  rotation  is  effected.  Then,  plane  gratings 
were  made  at  the  same  carrier  frequency  and  their  first  order  transmission 
response  was  also  determined  as  a  function  of  azimuth  angle. 

Results  for  the  two  recording  specimens  are  given  in  Fig.  3-19  for  (a)  the 
120-02,  6  pm  emulsion,  and  (b)  the  649  F,  17  ym  emulsion  respectively.  The 
figure  illustrates  the  (arbitrary)  sign  designations.  In  the  sign  designation,  (+) 
is  in  a  direction  which  decreases  the  reference  beam  angle.  The  (-3  dB)  maxi¬ 
mum,  first  order  efficiencies  for  the  two  silver  halide  emulsions  were  4.5  and  5.2% 
respectively.  For  these  two  cases,  the  angular  sensitivities  are  approximately 
±  14.6°  and  ±  8.7°,  respectively  for  the  thin  and  thick  emulsions.  The  grating 
responses  were  envelopes  for  the  filter  responses. 

If  allowance  is  made  for  the  experimental  variations  concomitant  with  the 
correlation  measurement  vs  the  direct  transmission  measurement  for  the  gratings, 
the  5x  matched  filters  are  seen  to  closely  follow  the  grating  response.  The  grat¬ 
ings  were  established  by  maximizing  the  interference  fringe  visibility  and  then 
establishing  an  exposure  bracket.  As  with  the  MF,  the  grating  selected  was  that 
which  gave  the  maximum  signal  upon  playback. 

The  Q-value  of  the  two  emulsions  shows  (Section  4)  thin  and  approximately 
thick  properties  respectively,  for  the  EK  120-02  and  EK-649F  photographic  emul¬ 
sions.  The  former  is  easily  classified  as  a  thin  grating  but,  in  the  latter  case, 
some  properties  of  thick  gratings  start  to  appear.  The  relatively  narrow  tuning  is 
evidence  of  this  although  unequivocally  thick  gratings  are  typically  1-2°  wide  at 
the  -3  dB  point.  The  closeness  of  the  MF  response  might  offer  some  evidence  that 
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-15.2°  +  29.5°  GRATING 
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1283-81-034(T) 


(b)  THICK  GRATING  AND  5X  MF 


Fig.  3-19  Azimuth  Angular  Response  of  Thin  and  Thick  Emulsions 
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the  visibility  of  the  gratings  has  also  been  properly  maximized  since  the  selection 
of  the  MF  on  the  basis  of  the  maximum  autocorrelation  signal  does  not  otherwise 
guarantee  this. 

Kobuta  (Ref  17)  has  shown  that  the  signal  and  reference  beams  are  recorded 
through  the  emulsion  depth  in  a  gradually  diminishing  way.  This,  of  course, 
contributes  to  fringe  plane  formation  in  the  thick  emulsion.  In  Kobuta's  model, 
the  fringe  visibility  is  relatively  constant  throughout  the  medium  and  its  average 
value  is  a  determinant  for  the  measured  diffraction  efficiency  which  is  a  maximum 
of  3.7%.  This  value  is  less  than  that  obtained  here  indicating  that  our  results  do 
not  entirely  fit  Kobuta's  model. 

In  many  applications  a  full  360°  coverage  of  an  object  is  necessary  in  the 
memory  to  account  for  any  possible  object  orientation.  The  angular  or  orienta¬ 
tional  sensitivity  then  dictates  the  number  of  views  necessary.  In  the  present 
case,  the  matched  filters  already  described  were  used  to  determine  this  sensi¬ 
tivity.  The  object,  an  isolated  single  tank,  was  rotated  and  the  autocorrelation 
signal  measured  over  an  interval  to  ensure  a  3  dB  decrease  in  magnitude.  The 
results  of  this  are  shown  in  Fig.  3-20.  Both  coherently  added  image  matched 
filters  display  a  (3  dB)  sensitivity  of  ±4.2°  with  good  symmetry,  while  the 
noncoherently  added  MF  shows  a  sensitivity  of  +6.9°  slightly  less  than  an  isolated 
single  MF  at  ±  8.0°.  Based  upon  these  results,  23  whole  filters  would  be  required 
of  the  single  image  variety  in  contrast  to  either  nine  or  six  whole  filters  needed 
when  the  5x  multiple  image  filters  are  constructed  in  either  a  noncoherent  or 
coherent  manner  respectively. 

Figure  3-21  illustrates  the  rotational  sensitivity  of  a  single  tank  MF  derived 
from  the  simulation  and  should  be  compared  to  (a)  the  sensitivity  of  the  single 
tank  MF  from  Ref  7,  and  (b)  to  the  results  for  the  multiple  image  MF  of  Fig.  3-20. 
It  can  be  seen  that  the  results,  allowing  for  slightly  different  spatial  frequency 
band  structures,  are  quite  similar. 

Highly  sensitive  single  MF  are  high  spatial  frequency  filters.  Thus,  one 
might  also  conclude  that  in  spite  of  fringe  visibility  and  autocorrelation  optimi¬ 
zations,  multiple  image  MF  inherently  tend  to  be  high  spatial  frequency  filters. 
Table  3-9  summarizes  orientational  results  for  coherent  addition  matched  filters. 

It  is  clear  that  they  all  illustrate  high  frequency  filter  characteristics  when 
allowance  is  made  for  experimental  variability. 


3-37 


NORMALIZED 

AUTOCORRELATION 


1.0 


DEGREES 


1283-81-036(T) 


Digital 


normalized 

autocorrelation 


1.0 


1283-81-03S(T) 


Analog 


Fig.  3-20  Rotational  Santitivity  of  Matchad  Filtar 
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TERRAIN:  SOUTHEAST  ASIA  WITH  5  VEHICLES 
WIF:  10X  COHERENTLY  ADDED 


Fig.  3-21  Correlator  Output  with  Terrain  Image  Input 
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TABLE  3-9  ORIENTATION  SENSITIVITY  OF  COHERENT  ADDITION  MATCHED 

FILTERS  (THIN  SIL  HAL) 


SIZE  OF 

COHERENT 

ADDITION 

ORIENTATIONAL 

SENSITIVITY 

IDEG.) 

IX 

17.76 

2X 

16.60 

3X 

17.15 

4X 

17.63 

5X 

14.90 

6X 

16.34 

7X 

15.9 

8X 

17.73 

9X 

15.30 

iox 

15.46 

1283-026(a)(T) 
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During  the  course  of  the  program,  coherent  arrays  continuously  up  to  lOx, 
and  at  15x,  were  used  to  make  matched  filters.  In  one  test,  the  entire  spread, 
lx  -  lOx,  was  used,  and  the  tank  angular  position  of  the  multiples,  which  gave 
the  maximum  autocorrelation  signal,  was  used  to  orient  the  SE  Asian  terrain  film 
upon  playback.  The  results  of  this  test  are  shown  in  Table  3-10.  The  values  are 
clearly  higher  than  an  average  S/C  ratio  for  the  size  of  the  memory.  Earlier  tests 
have  shown  fairly  conclusively  that  the  S/C  ratio  follows  the  autocorrelation  value 
in  the  multiple  case,  so  the  result  is  not  unexpected.  It  does  show,  in  conjunc¬ 
tion  with  other  results  given  earlier,  that  for  good  imagery,  uniform  S/C  ratios 
result,  and  that  high  ratios  are  possible  from  multiple  vehicle-terrain  input 
imagery.  Figure  3-21  illustrates  the  correlation  plane  for  lOx  image  MF  when  the 
woodland  and  the  SE  Asian  scene  are  played  through.  The  figure  shows  two 
threshold  levels,  clearly  illustrating  the  high  S/C  ratios  present. 

Finally,  there  were  a  series  of  tests  in  which  a  15x  image  matched  filter  was 
made  and  three  terrains  played  through  the  filter.  The  same  angular  position  of 
the  MF  was  used  in  each  case  to  orient  the  input  image,  thus  the  position  was 
somewhat  arbitrary.  The  scans  through  the  tank  correlation  are  shown  in  Fig. 
3-22.  Once  again,  system  acceptable  S/C  ratios  were  obtained.  The  filter  used  in 
these  tests,  shown  in  Fig.  3-23,  was  fabricated  in  dichromated  gelatin  and  can  be 
seen  to  be  a  high  pass  filter  in  the  expanded  views. 

3.3.3  Multiplexed  Carriers 

The  sequential  storage  of  holograms  with  different  reference  beam  angles  is 
referred  to  as  multiplexing  carriers  and  angle  codification.  Much  of  the  basic 
understanding  of  a  grating  recorded  in  a  thick  media  traces  to  Kogelnik’s  coupled 
mode  analysis  (Ref  14),  and  the  sequential  recording  to  Kowarschik  (Ref  21). 

The  work  was  later  extended  to  absorptive  media  by  Kubota  (Ref  17). 

The  region  of  interference  between  a  signal  beam  S  and  a  reference  beam  R  is 
three-dimensional  in  nature  so  that  when  a  recording  medium  is  used,  a  small, 
though  wavelength-appreciable,  slice  is  made  of  the  volume  of  interference.  In 
the  following  geometrical  model,  the  interference  is  exhibited  along  the  line  that 
the  recording  plane  makes  with  the  plane  containing  the  reference  and  signal  beam 
axes.  This  plane  is  the  matched  filter  fabrication  plane  of  symmetry. 
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TABLE  3-10  S/C  RATIO  FOR  MULTIPLE  COHERENT  MF  USING 
SE  ASIAN  TERRAIN  IMAGERY 


NO.  OF  IMAGES  USED 
TO  MAKE  MF 

S/C  RATIO 

IdB) 

IX 

13.66 

2X 

12.18 

3X 

10.99 

4X 

10.69 

5X 

7.99 

6X 

11.51 

7X 

12.88 

8X 

9.90 

9X 

7.60 

10X 

I283-027(«)(T) 

9.47 
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SE  ASIAN  TERRAIN.  12.1dB 


WOODLAND  TERRAIN.  9.3dB 


The  geometry  of  recording  is  such  that  the  plane  of  the  object  image  is  either 
on-axis  ( 0g  =  0,  0  r  t  0  =0  )  and  parallel  to  the  recording  plate,  or  at  an  angle  to 
the  plate  normal  so  that 

es  =  6r  =  T 

i.e.,  the  total  angle  is  the  same  in  both  cases. 

The  fringe  angle  e  in  the  medium  is 

e  =  ^RE  +  (3-7) 


where  Irj?,  IgE  refer  to  the  reference,  signal  beams  in  the  recording  media.  From 


Snell's  law  this  becomes 


e  =  (Sin  ^Sin  eR/n)  +  Sin  1(Sin  Ig/n))/2 


(3-8) 


where  the  incident  beam  is  given  as 


Ig  =  Tan  ((Sin  0g-(H/F)  Cos  6g  +  U/F))/(Cos  0g  +  (H/F)  Sin  0g))(3-9) 

The  spacing  d  of  the  fringes  (on  the  surface)  was  given  earlier  and  is 
obtained  from  the  basic  equation 


d  =  A  /(Sin  eR  ±  Sin  eg) 


(3-10) 


which  for  any  spatial  frequency  (or  physical  location  C)  requires  that  Ig  be 
substituted  for  eg.  0  r  is  constant  across  the  matched  filter  plane. 

The  results  were  given  in  Section  2,  Fig.  2-6.  A  sampling  of  results  for  a 
range  of  appropriate  spatial  frequencies  between  ±9.5  c/mm  is  shown  in  Table 
3-11.  These  values  of  spatial  frequency  represent  upper  (practical)  limits  of 
some  target  objects  used,  but,  by  no  means,  the  observable  upper  limits  in  a 
matched  filter. 

The  equations  for  the  fringe  angles  have  been  computed  for  the  two  general 
cases,  6g  =  9r;  0  g  =  0,  Br/0,  and  are  shown  in  Fig.  3-24  for  the  on-axis 
fringe  angle  (i.e.,  at  c  =  C/For  v  =  0).  The  limiting  values  are  also  shown  and 
represent  the  upper  limits  of  the  extreme  (practical)  case,  i.e.,  at  f-#  =  .6,  F  = 
25  mm  and  a  35  mm  input  imagery.  In  other  words,  these  are  the  limits  which 
cover  the  case  when  the  object  input  film  is  located  anywhere  in  the  field  during 
fabrication  for  the  particular  reference  beam  angle,  0r.  The  actual  range  of 
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TABLE  3-11  RANGE  OF  CARRIER  FREQUENCIES  FOR  A  SPATIAL  FREQUENCY 

SPREAD  OF  ±9.5  c/mm 


»S  +  *R 

Au  -  ±9.6  C/MM  AND:  | 

(DEGREES) 

®s“® 

®S  "  ®R 

10  1 

284  -  266 

284  -  266 

20 

S49  -  632 

669  -  641 

30 

800-781 

826-806 

40 

1020-1010 

1086-1075 

60 

1283-028(a)(T) 

1369-1361 

1587-1563 

FRINGE  ANGLE  (e).  DEG 
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Fig.  3-24  Limit!  on  Fringe  Angles  in  Thick  Media,  N  *  1 .53 
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values  for  for  a  matched  filter  would  be  an  increment  about  any  of  the  (H/F) 
lines  and  depend  upon  the  upper  spatial  frequency  band  limit  established  for  the 
filter.  The  specific  increment  would  be  obtained  from  the  equations  by  inserting 
±  values  for  c  =  £/F  for  a  particular  v  =  H/F  value. 

From  Fig.  3-24  it  can  be  seen  that  the  range  of  fringe  angles  is  greater  for 
the  case  0g  =  Or,  than  for  0g  =  0  even  when  allowance  is  made  for  a  total  angle 
comparison.  For  example,  when  0g  =  0r  =  15°,  e  =  9.7°  while  when  0r  =  30, 
flg  =  0,  e  =  9,5°.  Thus,  the  two  cases  are  not  strictly  identical  ones  for  fringe 
angle. 

Two  cases  were  examined  in  this  program,  with  0g  =  0,  ±  10°;  and  Or  = 

+10°,  +27°.  These  yield,  from  the  above  equations,  d  =  2.81,  2.81  pm  and  d  = 
2.81,  1.07  pm  vespectively  for  fringe  spacing.  Also,  we  have  e  =  ±3.1°  and  e  = 

+  3.1°,  +  8.7°  respectively  for  the  fringe  plane  angles  based  upon  Ng  =  1.53  and 

o 

X  =  4880  A.  Computing  the  Bragg  angle  requirement, 

Sine+  =  A/2Ngd  (3-11) 

we  find  rr  =  3°15',  quite  close  to  the  actual  angle  in  the  symmetrical  case 
considering  the  index  uncertainties. 

In  treating  sequentially  (noncoherent)  stored  gratings  in  transmission  volume 
holograms,  Kowarschik  (Ref  21)  gives  the  diffraction  efficiency  for  direct 
coupling  and  sinusoidal  modulation: 


n  =  nx  +  n2 


-  N 

|a  i 


js.l  2  +  Lil 

'  11  IaI 


(K19) 


where  the  Si  are  the  signal  waves  and  the  parameters  are  contained  in  the  Di,  A 
terms.  The  applicable  equations  are  given  in  Appendix  C.  Using  X=  4880  A,  N  = 
1.53,  thickness  =  13  pm  and  a  modulation  ni  =  1.5  x  10“2,  the  two  cases  were 
programmed  and  computed.  The  results  are  shown  in  Fig.  3-25  and  yield  HPFW 
angles  of  ±  15.5,  and  +  15.5,  +  16.2°  respectively  for  the  two  cases. 
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Fig.  3*25  Kowarschick  Formulation  Appliad  to  0R  -  +  10° 
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Gratings  were  prepared  in  dichromated  gelatin  and  their  responses  measured 
as  the  recording  plates  were  rotated.  In  both  cases,  the  detector  is  placed  along 
the  reference  beam  axis  as  the  plate  is  rotated.  These  results  are  shown  in  Fig. 
3-26  for  the  corresponding  cases.  The  angles  are  +13.1,  -12.7  and  +13.4,  +4.3° 
respectively.  The  symmetric  case  is  in  close  agreement  with  Kowarschick,  but  for 
the  asymmetric  case  there  is  considerable  disagreement.  Matched  filters  were 
fabricated  and  the  correlations  measured  under  the  same  conditions.  These 
results  are  shown  in  Fig.  3-27  where  the  HPFW  angles  are  +11.2,  +5.7°  respec¬ 
tively,  in  agreement  with  the  grating  test  results.  In  carrying  out  these  tests, 
experimental  restrictions  required  the  27°  angle  to  be  fabricated  first. 

In  order  to  verify  earlier  thin  emulsion  results,  MF  were  fabricated  at  the 
two  angles  ±  10°  and  the  correlation  measured.  These  results  are  included  in  Fig. 
3- 26 (a)  and  show  the  expected  uniform  response. 

Figure  3-28  illustrates  microscopic  views  of  a  5X  image  MF  in  dichromated 
gelatin  with  (a)  one  and  (b)  two  carriers,  specifically  in  this  case,  ±10°.  When 
the  two  MF  are  perfectly  registered,  or  nearly  so,  there  is  considerable  difficulty 
in  bringing  the  image  in  focus.  When  not  registered,  one  MF  and  then  the  other 
are  able  to  be  observed  by  slight  microscopic  adjustment. 

In  a  final  test  of  the  multiplexed  carriers,  a  15x  image  MF  was  constructed 
and  the  correlation  measured.  Then,  using  the  SE  Asian  scene  with  5  vehicles  in 
the  field,  the  S/C  ratio  was  measured.  In  the  first  case,  the  correlation  result 
(using  one  M-60  tank  position  to  measure  the  correlation)  shows  (Fig.  3-29)  the 
same  result  obtained  for  the  single  MF.  That  is,  the  27°  angle  has  a  higher 
response  and  an  angular  sensitivity  of  5.2°  while  the  10°  response  is  lower  with 
an  estimated  13°  angular  response. 

Secondly,  the  same  15x  matched  filter  was  used  for  the  terrain  scene 
evaluation.  The  result  is  shown  in  Fig.  3-30.  There  is  a  decrease  in  the  S/C 
ratio  as  well  as  an  expected  decrease  in  overall  signal  strength. 

Based  upon  the  foregoing,  one  can  obtain  usable  large  matched  and  multiplex 
two  such  filters,  and  obtain  acceptable  S/C  ratios.  However,  unless  the  MF  are 
symmetric  about  the  normal  signal  beam,  the  angular  response  of  the  filters  is 
uneven.  Furthermore,  the  basic  response  of  the  filters  leads  to  the  necessity  of 
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Fig.  3-27  Azimuth  Sensitivity  of  MF  With  Multiplexed  Cerriert 
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Fig,  3-28  5k  MF  in  DCG  With  One  and  Two  Carriers 
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using  two  correlation  planes.  The  Bragg  angle  requirements  for  maximizing  the 
signals  in  the  two  correlation  planes  can  be  met  simultaneously,  a  significant 
alignment . 
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4  -  MKMORY  RECORDING  MEDIA 


This  study  included  the  examination  of  materials  such  as  dichromated  gelatin 
(DCG),  photographic  silver  halide  emulsions  (SILHAL),  and  photopolymers  (PP), 
not  as  a  materials  study  but  as  a  study  of  recording  since  these  materials  repre¬ 
sent  individual  and  different  recording  mechanisms  (Fig.  4-1). 

Silver  halide  photographic  emulsions  for  amplitude  holograms  are  well  known 
and  are  an  inexpensive,  readily  available,  widely  used  storage  medium.  They 
are,  therefore,  a  valuable  optical  memory  resource.  Two  types  of  Eastman-Kodak 
SILHAL  plates  were  used  in  our  program:  120-02  and  649F.  The  nominal  thick¬ 
nesses  were  6  and  17#m,  respectively.  Depending  upon  wavelength  and  reference 
beam  angle,  the  latter  film  thickness  can  approach  and  display  some  character¬ 
istics  of  thick  films  (e.g..  Ref  14),  but  never  quite  be  classified  thick  by  the 
widely  accepted  Q-criterion  (Ref  15), 

Q  =  (2tt  At)/nQd2  (4-1) 

where  X  ,  t  are  the  wavelength  and  emulsion  thicknesses  respectively  for  a  film 
having  an  index  of  refraction  nQ  with  fringes  recorded  at  a  spacing  of  d.  The 
Q-factor  is  summarized  in  Fig.  4-2  for  the  film  and  fabrication  parameters 
employed  in  the  program.  The  factor  itself  has  recently  been  reassessed  so  that 
the  modulation  level  nj  is  also  included  (Ref  16). 

The  SILHAL  films  used  in  the  program  have  sufficient  dynamic  range  and 
resolution  to  record  wide  band  matched  filters.  Figure  4-3  illustrates  a  scan  with 
a  microdensitometer  of  the  eighth  to  fifteenth  order  of  a  MF  showing  that  wide 
band  spatial  frequency  matched  filters  can  be  recorded,  although,  as  shown  in 
Section  1,  the  memory  element  as  defined  has  about  one  half  of  the  spatial 
frequency  content  illustrated  in  Fig.  4-3. 

Some  differences  were  noted  between  the  thin  and  thick  silver  halide 
emulsions.  As  recording  media,  each  type  must  be  optimized  for  the  specific  set 
of  conditions  on  an  individual  basis.  For  example,  there  are  measurable  dif¬ 
ferences  between  exposure  time  and  smaller  differences  for  beam  ratio  adjustment 
between  thick  and  thin  SILHAL.  Table  4-1  illustrates  these  differences  for  lx, 
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Fig.  4-1  Memory  Recording  Mechanitim  for  Media 
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TABLE  4-1  POSITIONAL  DISTRIBUTION  OF  AUTOCORRELATIONS  OF  COHERENT  MATCHED 
FILTERS  USING  THICK  AND  THIN  EMULSIONS 


ORIENTATION 
ANGLE  <%fi>  DEG. 


CORRELATION  OF  COHERENT  IMAGERY 


IX 


THICK 


THIN 


5X 


THICK 


THIN 


JOX 


THICK 


THIN 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


1.000 


1.000 


.334 

.270 

.219 

.181 

.244 


.340 

.270 

.276 

.298 

.271 


.030 

.024 

.018 

.021 

.016 

.012 

.012 

.021 

.013 

.027 


.038 

.024 

.017 

.024 

.023 

.020 

.024 

.035 

.021 

.027 


N -AVERAGE 


1.000 


1.000 


.250 


.291 


.019 


.025 
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5x,  and  lOx  matched  filters.  In  addition  to  the  variations  between  each  element 
of  a  particular  coherent  image  array  (Table  3-6),  it  can  be  seen  that  there  is  some 
improvement  when  the  thinner  emulsion  is  used.  However,  it  has  also  been 
observed  that  when  constant  conditions  are  used  and  optimized  for  the  thick 
SILHAL  emulsion,  there  is  a  significant  increase  in  the  efficiency  of  the  thick 
emulsion  over  the  thin  due  to  the  Bragg-like  improvement. 

Figure  4-4  shows  a  200X  scanning  electron  microscope  (SEM)  picture  of  a  5x 
coherently  added  image  matched  filter  with  a  blowup  of  one  section.  The  details 
of  the  surface  relief  in  the  DCG  should  be  noted  and  compared  with  the  microscope 
picture  of  the  5X  in  the  previous  section. 

Matched  filters  were  selected  which  were  of  comparable  quality  but 
fabricated  in  silver  halide  and  diehromated  gelatin.  They  were  then  compared  in  a 
common  test  for  their  angular  response  and  signal-to-clutter  ratio.  The  results 
are  shown  in  Table  4-2. 

Basically,  their  performance  is  comparable  in  terms  of  response  quality.  In  a 
separate  test,  the  efficiency  of  DCG  vs  SILHAL  was  determined  and  it  was  found 
that  DCG  was  lOx  more  efficient  than  a  comparable  filter  in  thin  silver  halide. 

DCG  is  most  interesting  because  the  molecular  nature  of  the  medium  has  the 
potential  of  high  efficiency  without  (silver)  grain  scattering.  This,  in  turn, 
should  lead  to  memories  with  the  maximum  extractaole  capacity.  The  mechanism  of 
recording  is  primarily  one  of  producing  a  variation  of  the  surface  contour.  We 
have  drawn  upon  contemporary  AFOSR  sponsored  programs  (Ref  18  -  20)  during 
these  efforts  in  addition  to  extending  some  of  the  technology.  A  detailed  review 
of  our  experiences  in  establishing  appropriate  procedures  is  contained  in 
Appendix  A. 

Initially  it  was  intended  that,  in  addition  to  the  two  materials  cited,  some 
form  of  photopolymer  be  included  in  order  to  extend  the  work  to  media  in  the 
thickness  range  20  ^  th  <”  100  micrometers.  A  search  was  conducted  to  determine 
likely  candidate  materials,  the  most  promising  of  which  were: 

Urethane  aligomers  (Thiokol) 

Acrylic  allylic  (R  and  H) 

Methacryllic  allylic  (R  and  H) 

Methyl  Cyanoacrylates  (Loctite)  (Kef  27) 

Unspecified  monomer-initiator  system  (duPont) 

Acrylate-organic  acid  initiators  (NRC) 
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TABLE  4-2  MATCHED  FILTER  MEDIA  COMPARISON 


■m 

B 

# 

IMAGES 

COH/NON 

THICK/THIN 

DCG/SILHAL 

S/C  IdB) 

ANGULAR  RESPONSE 

A 

IX 

- 

THIN 

DCG 

10.6 

N/A 

IX* 

- 

THIN 

SILHAL 

11.2 

±8.0’ 

B 

5X 

COH 

THIN 

DCG 

7.2 

±6.3° 

5X 

COH 

THIN 

SILHAL 

7.0 

t6.T 

C 

tx 

- 

THICK 

DCG 

6.1 

±9.1° 

12B3-030(»)< 

IX 

r) 

— 

THICK 

SILHAL 

6.5 

±7.8° 
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Several  practical  factors  dictated  preferences  in  choosing  materials.  When 
such  factors  as  activation  wavelength,  fabrication,  and  availability  were 
considered,  only  the  last  two  offered  any  reasonable  promise  within  the  program 
objectives.  Use  of  the  discontinued  duPont  material  (Ref  26)  further  limited  the 
list  when  the  only  material  available  was  found  to  be  several  years  old  and 
unreactive.  Thus,  the  best  choice  was  the  last  material  (Ref  28). 

A  number  of  plates  were  prepared  in  the  usual  way  with  25  urn  and  50  um 
thicknesses,  both  of  grating  and  matched  filter  holograms. 

Serious  problems  were  encountered  while  working  with  the  acrylate-organic 
acid  initiators.  In  all  cases,  the  recordings  were  extremely  noisy,  inadequately 
defined  holographically,  and  were  not  able  to  be  made  permanent  in  their  recorded 
state.  A  number  of  experimental  procedures  and  techniques  were  devised  to 
improve  performance.  After  an  extended  investigation  of  this  material,  it  became 
clear  that  the  prospect  of  obtaining  memory  elements  in  the  best  choice  of  photo¬ 
polymers  with  acceptable  quality  was  not  possible  within  the  scope  of  the 
program,  and  the  effort  was  discontinued. 

Our  experience  with  photopolymers  is  not  surprising.  Interest,  as  evidenced 
by  a  paucity  of  technical  papers,  reports,  and  a  very  limited  commercial  availa¬ 
bility,  infer  poor  performance  of  these  materials  at  the  present  time.  Future 
developments,  or  a  dedicated,  in-depth  study  of  available  materials  may  someday 
result  in  acceptable  MF. 

Based  upon  the  parameters  used  in  the  program,  or  those  postulated  as  limit 
conditions,  the  efficiency  of  the  thick  medium  was  examined  with  the  results  of 
Kogelnick  (Ref  14)  as  a  starting  point.  Alferness  (Ref  23)  and  others  (Ref  24,  25) 
have  examined  the  angular  sensitivity  of  gratings  in  thick  media  in  order  to 
establish  intervals  between  angles  to  provide  angular  separability,  or  isolation. 

The  efficiency  is  given  by 


n  =  (Sin(v2  +  +  $2/v2) 


(4-2) 


where  it  can  be  shown  that 


^  _  [Cos  (<fr~e) 


(  A/A  >  (1/2N)  2ir  d 


[A  Cos  6  -  (X/N)  Cos  4> ] 


(4-3) 
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and 


v  =  (irn,  d/A)  [Cos  6  (Cos  6  -(A/A))(Cos  <j>/N)] ~ 1  (4-4) 

A  number  of  parameters  have  been  previously  defined  but,  for  convenience, 
all  are  given  in  Table  4-3  with  typical  values  used  in  the  computations. 

In  carrying  out  the  computations  we  are  again  concerned  with  the  two  cases, 

0 s  =  0,  0r,  and  0S  =  *  0* 

First,  the  efficiency  was  determined  as  a  function  of  reference  beam  angle 
with  thickness  and  modulation  taken  as  constants.  One  value  for  the  latter  (i.e., 
.01)  was  taken  close  to  the  modulation  value  of  Kowarschik  as  used  in  the  pre¬ 
vious  section.  The  results  for  es  =  0  and  es  =  are  shown  in  Fig.  4-5  and  4-6 
respectively. 

In  the  first  case,  the  influence  of  modulation  for  a  given  thickness  as 
discussed  by  Moharam  (Ref  16)  is  apparent.  In  each  instance  the  thicker  medium 
yields  the  larger  efficiency  with  the  P  =  20,  N  =  .01  combination  almost  constant 
for  any  reference  beam  angle  of  interest.  The  lower  modulation  yields  the  higher 
efficiency  in  each  of  the  two  groups. 

The  second  case  (Fig.  4-6)  was  the  usual  6S  =  0r  where  the  total  angles 
would  be  the  same  as  the  0S  =  0  case  above.  The  abscissa  is  expanded  to  reflect 
this  quality.  It  is  apparent  that  care  must  be  exercised  with  the  symmetrical 
configuration  in  selecting  an  appropriate  reference  beam  if  a  low  pass  MF  is 
desired  (i.e.,  where  v  =  0,  or  near  it). 

Due  to  the  uncertainties  in  DCG  index,  the  efficiency  of  the  MF  for  variations 
in  index  of  refraction,  N,  was  also  examined  at  a  constant  modulation  N  =  .04.  It 
was  found  that,  while  the  efficiency  varied  with  reference  beam,  it  was  fairly  con¬ 
stant  with  index.  Thus,  an  average  value  was  taken  and  is  summarized  in  Table 
4-3.  It  can  be  seen  that,  like  the  results  of  Fig.  4-5,  4-6,  there  is  a  significant 
difference  in  efficiency  between  the  10  and  20  pm  cases,  as  well  as  between  the 
asymmetrical  (S  =  0)  and  symmetrical  (S  =  R)  cases. 

Finally,  in  the  same  group  of  computations,  the  variation  across  a  matched  fil¬ 
ter  (|/F)  was  examined  for  fixed  index  and  modulation.  These  results  are  shown 
in  Tables  4-4  and  4-5  and  reveal  a  significant  change  across  a  filter  for  the  large 
reference  beam  angles,  and  a  larger  variation  for  the  symmetrical  case,  S  =  R. 
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TABLE  4-3  PARAMETERS  FOR  EQUATIONS  4-2,  -3.  -4 


0 

IS  THE  ANGLE  OF  THE  GRATING  VECTOR  K  BUT  ALSO, 

(ff/2-el  WHERE  e  IS  THE  FRINGE  PLANE  ANGLE  GIVEN  EARLIER 

a 

ANGLE  OF  INCIDENCE  IN  THE  MEDIUM; 

9o 

= 

BRAGG  ANGLE  MEASURED  FROM  FRINGE  PLANE,  SIN9BR  -  X/2NT 

&e 

DEVIATION  FROM  BRAGG  ANGLE,  OR  DIFFERENCE  BETWEEN  REQUIRED  BRAGG  ANGLE 

AND  INCIDENT  ANGLE  IN  MEDIUM, 

N 

= 

REFRACTIVE  INDEX  (1.56), 

Nl 

* 

MODULATION  (0.01  -  0.04); 

P 

* 

MEDIA  THICKNESS  (10-20  Mm). 

\ 

= 

FABRICATION  OR  PLAYBACK  WAVELENGTH  AS  REQUIRED; 

c 

= 

SLANT  FACTOR,  THE  RATIO  OF  THE  OBLIQUITY  FACTORS  -  Cp/Cg 

CR 

= 

cosa 

CS 

rosa  -  KCOS  0 

0 

= 

COS0  -*/ON  SINe 

D 

= 

FRINGE  SPACING 

K 

- 

2ir/\ 

0 

2*NI\ 
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.04  INDEX  INTERVAL, 


TABLE  44  AVERAGE  EFFICIENCY  (%)  FOR in,  - 

1.52  - 1.56 


(W  FOR  S  -  R 

1Z83-032(»)(T) 
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TABLE  4-5  AVERAGE  EFFICIENCY  (%)  FOR  n-j  »  .04  AT  N  *  1.55 
FOR  VARIOUS  v 


v  OR  NORMALIZED  DISTANCE.  g/F 


9n 

0.0168 

0 

-0.0175 

10 

78.99 

80.02 

79.88 

30 

55.20 

61.82 

61.05 

60 

0.50 

5.68 

12.42 

(a)  FOR  S  =  0 


v  OR  NORMALIZED  DISTANCE,  g/F 


®R 

0.0175 

0 

-0.0175 

5 

51.94 

62.26 

15 

15.21 

24.72 

30 

0 

0 

1.06 

283-033(A)(T) 


lb)  FOR  S  =  R 
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Fig.  4-5  Efficiency  of  Recording  Medium,  0$  -  0 
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Fig.  4-6  Efficiency  of  Recording  Medium,  6$  =  0R 
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The  net  experience  then  is  that  either  thin  or  thick  silver  halide,  or  dichro- 
mated  gelatin  can  be  satisfactorily  used  to  record  correlator  memory  elements. 

The  trade-offs  between  them  are  not  to  be  based  upon  the  performances  as  used  in 
the  program,  but  upon  other  factors  such  as  cost,  speed,  etc.  Given  sufficient 
resources,  the  dichromated  gelatin  is  to  be  preferred  for  large  memories. 


4-16 


5  -  OPTICAL  MEMORY  CAPACITY 

5.1  CAPACITY  BASED  UPON  FREQUENCY  CONTENT  OF  MF 

In  addition  to  specific  techniques  for  storing  matched  filters  as  described  in 
detail  earlier,  one  is  also  interested  in  the  number  of  such  techniques  which  might 
concurrently  be  used  on  the  optical  memory  capacity. 

Figure  1-2  illustrated  the  basic  storage  element.  The  region  defined  by  the 
element  is  circular  with  a  radius  equal  to  N  times  the  optimizing  spatial  fre¬ 
quency,  v0,  which,  for  the  M-60  tank,  was  3.33  cycles /mm.  It  is  believed  that 
each  vehicle  would  be  somewhat  different.  In  the  case  of  the  M-60,  N  =  2.5  to  3,  a 
spatial  frequency  cut-off  in  which  the  autocorrelation  signal  is  reduced  to  0.27  dB 
of  its  "infinite  bandwidth"  level. 

In  a  regular  array  then,  the  cross  talk  does  not  exceed  6%.  One  matched 
filter  requires  (N  vQ/ S)  mm2  area  where  S  =  (  XF)“1  is  the  "system  constant."  As 
before,  X  =  wavelength  of  the  Fourier  transformation  and  F  the  focal  length.  If 
we  include  G  =  the  geometrical  packing  factor,  S  =  the  split  spectrum  factor,  and 
P  =  the  overlap  or  multiplicity  factor,  then  the  capacity  for  matched  filters  of  a 
single  species  of  target  becomes: 

C.  =  (25sPG)/(Nv  AF)2  matched  filters/cm2  (5-1) 

A  O 

where  vQ,  X  ,  and  F  are  given  in  cycles/mm,  mm,  and  mm  respectively.  Using  the 
realistic  values:  S  =  1,  P  =  5,  G  =  1.16,  N  =  3.0,  X=  0.488  m,  and  F  =  25  mm,  then, 
for  example, 

CA  =  2670  (matched  filters/cm2). 

Table  5-1  illustrates  some  values  of  unit  capacity  for  some  practical  system 
cases.  Any  value  can  be  appropriately  modified  by  multiplying  the  capacity  given 
by  an  overlap  or  multiplicity  factor,  packing  density,  and  so  forth. 

It  is  apparent  that  a  large  number  of  matched  filters  can  be  packed  in  a  small 
area.  It  is  also  apparent  that  when  one  considers  large  memories,  fabrication  goes 
beyond  the  realm  of  manual  processes  and  computer  controlled  or  computer 
generated  memories  are  a  necessity. 


5-1 


TABLE  5-1  UNIT  CAPACITY  FOR  SOME  PRACTICAL  CASES 


FOCAL  LENGTH 

F  (mm) 

WAVELENGTH 

(ANGSTROMS) 

CAPACITY 

CA  (MF/Cm2) 

25 

4880 

1683 

6328 

1001 

200 

4880 

26 

6328 

16 

360 

4880 

8.1 

6328 

4.8 

720 

4880 

2.0 

6328 

1.2 

1000 

4880 

1.05 

128 3*0 34(a)  (T) 

6328 

0.63 

Secondly,  there  is  a  power  division  with  increasing  size  requiring  large  laser 
power.  Earlier  (Ref  7),  it  was  shown  that  the  power  level  employed  at  that  time 
was  10~4  what  might  be  used  before  emulsion  damage  occurred  for  the  media 
employed  (silver  halide,  bleached  silver  halide).  In  the  current  program,  the 
limits  on  dichromated  gelatin  damage  threshold  were  not  established.  However,  it 
was  determined  that,  as  expected,  it  is  well  above  that  of  silver  halide. 

Thus,  it  is  clear  that  there  are  some  choices  available  when  a  specific  memory 
is  designed,  size  requirements  developed,  the  optics  established,  and  a 
wavelength  chosen. 

The  formula  above  can,  incidentally,  be  applied  to  bit-by-bit  memories  and, 
with  acceptable  state-of-the-art  values,  yields  10®  bits/cm^,  a  value  often  quoted 
in  the  literature. 

The  overlay  or  multiplicity  factor  referred  to  above  does  not  take  multiplexed 
carriers  into  consideration.  A  factor  accounting  for  the  impact  of  this  has  not 
been  incorporated  because  of  the  need  to  have  additional  correlation  planes.  Yet, 
the  ability  to  use  multiplexed  carriers  to  increase  the  memory  density  is  real. 

Based  upon  Eq  3-11,  one  can  at  least  determine  the  limits  on  multiplexing  as 
they  affect  the  capacity.  Figure  5-1  illustrates  the  range  of  fringe  angle  planes 
for  the  practical  situation  of  processing  a  35  mm  format  with  transform  lens  focal 
lengths  of  25  to  1000  mm  for  reference  beam  angles  of  5°  and  45°.  The  curves  for 
both  are  similar,  as  are  the  limit  curves  for  any  reference  beam  except,  of  course, 
for  the  asymtotic  limits.  Also,  negative  reference  beam  angles  are  symmetric. 

The  angles  previously  computed  for  the  multiplexed  symmetrical  case  can  be  seen 
to  lie  within  the  ±  10°  reference  beam  curves. 

5.2  CAPACITY  BASED  UPON  FILM  LIMITS 

The  use  of  silver  halide  as  the  recording  medium  raises  the  question  of 
capacity  limits  restrained  by  finite  resolution  of  the  medium. 

An  upper  limit  to  the  number  of  storage  addresses  in  a  holographic  optical 
memory  due  to  the  finite  resolution  capability  of  the  recording  medium  used  to 
store  the  holograms  may  be  established.  Assume  that  the  recording  medium 
response  to  the  source  of  illumination  used  in  the  holographic  fabrication  process 
is: 

R(x,y)  =  Eq  +  &  E(x,y)*m(x,y)  (5-2) 
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where  R(x,y)  denotes  the  spatial  variation  of  the  response,  E0  is  the  average 
illumination,  E(x,y)  is  the  information-carrying  portion  of  the  illumination 
(assumed  to  be  a  small  deviation  about  the  average  illumination),  and  E*m(x,y) 
indicates  a  linear  relationship  with  a  fixed  point  spread  function  m(x,y)  whose 
Fourier  transform  M(fx,fy)  represents  an  equivalent  MTF  of  the  recording 
process.  If  the  storage  medium  yields  an  amplitude  hologram  it  will  be  assumed 
that  the  transmission  through  the  hologram  is  proportional  to  the  exposure  during 
fabrication,  i.e.,  tR(x,y).  If  the  media  yield  a  phase  hologram  it  is  assumed  that 
the  transmission  is  given  as 

t  =  ke]kR(x,y)  -  K(1  +  jkR(x.y))  (5-3) 

which  is  based  on  a  simple  paraxial  treatment  of  a  thin  phase  hologram  with  small 
phase  deviation  justifying  the  resulting  approximation  to  the  exponential.  In 
either  case  (amplitude  or  phase),  the  information-carrying  term  in  the  resulting 
transmission  function  for  the  hologram  is  proportional  to  the  medium's  response 
R(x,y). 


If  a  standard  Vander  Lugt  approach  is  used  to  fabricate  a  hologram  for  the 
Fourier  transform  of  the  matched  filter  function  h(x,y),  then  one  can  express  the 
transmission  of  the  hologram  by  the  following  (a  modified  version  of  the  typical 
expression  due  to  the  presence  of  finite  resolution  of  the  medium  (Ref  29)): 


t<4.&)  =  r„2  +  +—  (He>2™Vm  +—  (5-4) 


where  H(x,y)  is  the  Fourier  transform  of  h  obtained  at  a  wavelength  with  a  lens 
of  focal  length  f,  a  =  sin  0/X  with  e  the  offset  angle  between  the  reference  and 
image  beams  used  in  the  recording  process,  and  &  =  fxXf,y  =  fyXf  are  the  physical 
spatial  variables  expressed  in  terms  of  the  Fourier  transform  frequencies  fx  and 
fy.  If  the  hologram  is  played  back  with  a  beam  that  is  modulated  by  the  image 
y(x,y),  then  the  portions  of  the  inverse  transformed  output  that  corresponds  to 
the  correlations  between  g  and  h  are  given  by 

F-1{G(£,y)[R(&,$)e"^2iray*m(£,y)] } 


=  F  X{e  ^2lT0tyG(x,y)[H*(e^2lTaym)] } 


(5-5) 


=  g(x,y)*[h(-x,-y)M(x,y  +  aXf)] 


where  x,  y  are  the  spatial  coordinates  in  the  output  correlation  plane  with  origin 
taken  at  0  and  a  X  f  relative  to  the  optic  axis  of  the  inverse  transform  lens. 

In  order  to  achieve  the  desired  matched  filter  operation,  the  output  should 
be  g*h(-x,-y),  but,  due  to  the  finite  resolution  capability  of  the  recording 
medium,  it  is  modified  by  the  presence  of  the  MTF  function  (M).  If  we  assume 
that  M  is  in  the  form  of  a  low -pass  filter  with  nominal  bandwidth  Bm,  then,  in 
order  to  avoid  excessive  alterations  of  the  matched  filter  weighting  function,  it  is 
necessary  that  a  <  Bm,  as  is  clear  from  Fig.  5-2  where  we  assume  that  the  spatial 
extent  of  the  matched  filter  h  is  small  compared  to  the  recording  process 
bandwidth. 

It  has  been  tacitly  assumed  that  the  reconstruction  of  the  correlation  plane  is 
achieved  without  distortion,  and  the  well-known  constraint  on  the  reference  beam 
angle  9  r  necessary  to  avoid  overlap  with  the  dc  spot  requires  that 

«XF  >  (3  Wh  +  2  Wg]  (Ref  29) 

where  Wh  and  Wg  are  the  spatial  width  of  the  matched  filter  image  h,  and  the 
input  image  g.  Combining  this  restriction  with  the  previous  one  (that  avoids 
distortion  due  to  medium  resolution)  leads  to  the  inequality  BmXf>aXf >3/2 
Wh  +  Wg. 

Now  the  spatial  width  of  h  is  inversely  proportional  to  the  nominal  bandwidth 
of  the  matched  filter,  i.e.,  =  1/Bn  where  Bn  is  the  matched  filter  bandwidth. 

In  the  optical  system,  the  matched  filter  bandwidth  Bh  corresponds  to  a  physical 
dimension  d  =  Bh  X  f  =  X  f/Wj,.  In  general,  only  a  portion  of  the  entire  spectrum 
of  the  matched  filter  is  recorded  in  the  holographic  system.  If  it  is  assumed  that 
is  is  necessary  to  record  a  total  matched  filter  spectrum  that  is  K  times  its  nomi¬ 
nal  bandwidth ,  then  this  would  require  a  linear  dimension  in  the  hologram  trans¬ 
form  plane  given  by  2  kd  =  2k  X  f/Wh,  which  implies  that  the  number  of  matched 
filters  that  could  be  stored  per  unit  areas  (without  multiple  overlapping  tech¬ 
niques)  is  given  by  (l/2kd)2  =  (W^/2k  X  f)^.  The  effect  of  the  recording  media 
now  becomes  apparent  since  the  previously  derived  inequality  implies  that 

W.  *  B 

_a  < _ m  - 

2Xf  3  +  2  W  /W. 

e  h 

and  when  combined  with  the  expression  for  the  number  of  filters  per  unit  area 
this  leads  to 
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A)  TRANSFER  FUNCTION  OF  THE  MEDIA 
RECORDING  PROCESS 
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B)  SITUATION  OEPICTING  SEVERE  DISTORTION 
OF  MATCHED  FILTER  IF  o>Bm 
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Cl  SITUATION  DEPICTING  LITTLE  DISTORTION  WHEN  a>Bm 
1283-81-053(T) 

Fig.  5-2  Effect  of  Finite  Reiolution  of  Media  on  Matched  Filter  Weighting  Function  h 
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Number  of  matched  filters  per  unit  area  < 


[2 


|K  (3+2  Wg/Wh)j 

where  Bm  is  the  bandwidth  of  the  recording  medium,  and  K  is  the  number  of 
"matched  filtf^  bandwidths"  that  are  to  be  stored  in  the  memory  plane.  Wg  is  the 
spatial  length  of  the  input  image  (or  aperture  size)  and  is  the  length  of  the 
matched  filter  image.  As  an  example  we  take  Bm  =  1000  cycles/mm  and  K  =  4  to 

104 

yield  the  number  of  matched  filters  per  unit  area  <  - g  equal  106  MF/cm2 

3  mrn 

when  Wg/ Wn  <<  2  or  10*  MF/cm2  if  Wg  =  lOW^  (MF  image  size  1/10  of  the  size  of 
the  aperture). 
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6  -  SUMMARY 


It  has  been  demonstrated  that  large  optical  memories  for  correlator  appli¬ 
cations  can  be  achieved  using  various  techniques  for  storing  the  spectral  signa¬ 
ture,  such  as  geometrical  packing,  split  spectrum,  coherent  and  noncoherent 
addition,  and  multiplexed  carriers.  Furthermore,  in  one  memory  position, 
sufficient  memory  elements  can  be  stored  to  cover  all  orientation  angles  of  the 
representative  M-60  tank  target  for  one  to  two  range  bins.  Acceptable  S/C  ratios, 
>  5  dB,  have  been  experienced  under  these  conditions.  Using  realistic  values  for 
focal  length  and  wavelength,  166  matched  filters  could  be  stored  on  a  single  filter 
basis.  The  use  of  any  of  the  efficient  storage  techniques  should  increase  this 
number  by  a  factor  between  1.16  and  30.  However,  as  the  number  of  memory 
elements  increases,  the  laser  energy  output  per  element  decreases  for  a  fixed 
input  power  level. 

To  increase  the  efficiency  of  recording  and  the  possibility  of  large  memories, 
silver  halide,  dichromated  gelatin  and  photo  polymers  were  investigated. 

Silver  halide  was  an  extensively  used  and  well-known  quantity  in  the 
program,  and  often  compared  with  dichromated  gelatin. 

DCG  proved  to  be  a  highly  efficient,  68.9%,  low-noise  medium  for  recording 
matched  filters.  When  tested  with  terrain  scenes,  S/C  ratios  ranged  from  5  to  12 
dB  depending  upon  the  particular  terrain  scene  used  and  the  size  of  the  memory 
element.  DCG  is  susceptable  to  relative  humidity  (>50%)  and  must  be  otherwise 
preserved,  usually  by  sealing  a  cover  plate  on  it.  This  technique  was  used. 

Photo  polymers  recording  is  noisey  and  inefficient,  and  requires  additional 
advances  in  technology  before  the  matched  filter  performance  is  acceptable. 
Therefore,  DCG  as  a  medium  for  optical  correlator  memories  is  recommended. 

The  multiplexed  carrier  technique  has  been  successfully  used  with  the 
semi-thick  (~15#*m)  dichromated  gelatin  recording  medium.  Thirty  memory 
elements  were  stored  in  one  position  with  a  5.3  and  6.0  dB  S/C  ratio  for  each  of 
the  carriers. 
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The  matched  filter  can  be  viewed  as  a  "chirp"  grating  whose  performance  is 
dependent  upon  the  size  of  the  memory  element.  As  the  number  of  coherently 
added  images  to  make  a  matched  filter  increases,  the  low  frequency  region  of  the 
filter  is  saturated.  The  MF  acts  as  a  high  pass  filter.  This  feature  probably 
accounts  for  the  favorable  discriminating  nature  of  large  coherent  arrays.  By 
contrast,  noncoherently  added  MF  have  a  frequency  response  of  a  single  image  MF. 

Noncoherent  additions  up  to  five  and  coherent  additions  up  to  15  have  been 
successfully  demonstrated  but  are  not  the  upper  limits  of  each  technique. 

Rather,  program  and  time  constraints  limited  explorations  beyond  these  values. 
There  are  indications  however,  that  the  limits  are  not  much  higher. 

The  translational  sensitivity  of  large  memory  elements  (e.g.,  5x,  lOx)  is 
comparable  to  that  of  a  single  MF.  Thus,  for  system  considerations,  no  greater 
locational  tolerance  needs  other  than  usual  are  introduced.  Target  azimuthal 
response  shows  that  the  thin  MF  media  (~6  pim)  have  broad  correlation  responses 
and  lower  efficiencies  than  the  more  sharply  tuned  seinithick  media  (~15/xm). 
Further  analyses  show  that,  when  thicker  media  ( — 25  p.m)  are  used,  greater 
efficiencies  should  result. 

Successful  use  of  a  computer  model  has  been  accomplished  in  the  program. 
This  has  enabled  a  number  of  extended  analyses  to  be  performed  in  a  more  timely, 
practical  manner.  For  example,  multiple  target-in-terrain  situations  can  be 
readily  simulated  and  analyzed  in  this  manner.  Also,  some  detailed  examinations 
of  matched  filter  geometries  can  be  performed  in  a  consistent,  repeatable  manner. 
Figure  6-1  illustrates  a  digital  version  of  a  multiple  image  coherent  array  MF  in 
which  the  low  frequency  cutoff  has  been  accurately  established.  The  impulse 
response  of  the  filter  is  also  shown. 

A  significant  advantage  of  an  optical  processor  over  alternate  means  is  size. 
To  illustrate  compactness,  a  demonstration  system  was  put  together  with  existing 
correlator  elements  from  bench  components.  This  consolidated  "system"  was 
assembled  and  tested.  A  2  x  2  multiple  holographic  lens  was  used  to  fabricate  the 
corresponding  2x2  array  of  matched  filters,  which  were  two  5x  and  two  lx  MF. 
The  correlation  plane  and  the  MF  array  are  shown  in  Fig.  6-2.  The  southeast 
Asian  scene  (Fig.  l-3a)  was  the  object  input  film.  A  S/C  ratio  of  7  dB  was 
obtained  from  the  unit. 
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(bl  PHOTO  OF  2x2  MF  ARRAY 
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Fig.  6-2  2x2  Arriy  &  Correlation  Rasults 
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Finally,  in  Fig.  6-3  we  have  tabulated  various  factors  that  can  influence 
optical  memory  systems  for  correlator  applications.  The  list  is  general  in  nature. 
It’s  purpose  is  to  summarize  our  experience  and  design  considerations  throughout 
the  course  of  this  program. 


1  GENERAL 

II  CONSTRUCTIONAL 

•  MODE:  WRITE 

WRITE,  READ 

WRITE.  READ,  ERASE 

•  RECORDING  SPEED 

•  DYNAMIC  RANGE 

•  DAMAGE  THRESHOLD 

•  SIZE 

•  CAPACITY 

•  POTENTIAL  FOR  REPLICATION 

•  WRITING  ENERGY 

•  PACKING  DENSITY 

•  POSITIONAL  ACCURACY 

•  ABILITY  TO  RECYCLE 

•  COST 

•  ABILITY  TO  ADD-ON 

•  DEFECT-FREE 

•  PROCESSING  NEEDS 

III  ENVIRONMENTAL 

IV  OPERATIONAL  t  PERFORMANCE 

•  SHELF  LIFE 

•  ACCESS  TIME 

•  STABILITY 

•  RESOLUTION 

•  INTEGRITY 

•  SENSITIVITY 

•  ARCHIVAL  LIFETIME 

•  LINEARITY  ACROSS  FIELD 

•  DIFFRACTION  EFFICIENCY 

•  TEMPERATURE  RANGE 

OF  OPERATIONISTORAGE 

•  S/C  OR  S/N  RATIO 

•  SHOCK,  VIBRATIONAL  SURVIVABILITY 

•  SPATIAL  FREQUENCY  BANDWIDTH 

1283-81-056<T) 

Fig.  6-3  Optical  Memory  Requirements 
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APPENDIX  A 


USE  OF  DICHROMATED  GELATIN  FOR  MATCHED  FILTERS 


The  use  of  dichromates  as  light  sensitive  activators  of  organic  matters  dates  back  to 
about  1830  (Ref  1,2).  Although  old,  it  has  in  the  last  10  years  been  suggested  and 
examined  as  a  holographic  medium  (Ref  3,4)  and  applied  to  such  applications  as  heads- 
up  displays  (Ref  5).  Interest  in  this  medium  for  large  capacity  optical  memory  stems 
from  its  potential  advantages  in  recording  efficient  matched  filters.  It  has  large  modu¬ 
lation  capability,  essential  when  overlapped  filters  are  required.  Since  the  rp>x;hanism 
is  molecular  in  nature,  it  is  not  bound  by  the  spatial  frequency  limiting  aspects  of  silver 
halide  grain  size,  and  thus,  as  a  medium  it  may  enable  the  limits  on  some  parameters  to 
be  utilized.  The  molecular  nature  of  the  dichromated  gelatin  also  suggests  low  scatter¬ 
ing  and  thus,  the  potential  for  phase  type  matched  filters  is  there  without  the  inherent 
reduction  in  S/C  ratio  when  photographic  emulsions  are  bleached. 

Although  there  are  many  advantages,  the  dichromated  gelatin  recording  medium  suf¬ 
fers  two  drawbacks,  neither  insurmountable.  First,  the  material  has  about  1/10  the 
sensitivity  of  slow  photographic  emulsions  and  requires  a  blue  excitation  (unless  other¬ 
wise  dyed  and  hardened  as  described  in  Ref  6,7).  Secondly,  the  material  like  bromine 
bleached  photographic  emulsions  is  susceptible  to  degradation  due  to  high  relative  humid¬ 
ity.  This  is  surmountable  by  environmental  control  or  the  use  of  an  optical  quality 
cover  plate. 

The  hardening  process  is  a  complex  chemical  reaction  which  after  many  years  is 
still  not  entirely  understood*,  a  situation  not  surprising  when  one  considers  the  organic 
and  variable  origins  of  the  gelatin.  In  the  current  program  it  is  not  essential  that  the 
process  be  understood  or  further  unraveled  but  rather  that  a  consistency  be  maintained 
so  that  at  least  the  process  can  be  duplicated  in  terms  of  holographic  formation.  Our 
other  objective  in  examining  this  medium  is  to  determine  whether  the  phenomenology  of 
surface  relief  inherent  to  the  medium  does  in  anyway,  advance  the  capability  for  large 
capacity  memories. 


♦An  excellent  recent  study  (Ref  8)  has  been  conducted  in  which  through  parametric 
control,  the  mechanisms  arising  in  the  dichromated  gelatin  process  are  explored  in 
depth. 


The  undoped  dichromated  gelatins  have  high  absorptivity  to  blue  and  pv  radiation 
and  negligible  absorption  to  the  red.  This  prescribes  argon  or  helium-cadmium  lasers 
for  fabrication  although  playback  can  be  at  any  wavelength  with  appropriate  compensa¬ 
tions.  The  indices  of  refraction  as  given  by  Shank  off  (Ref  4)  are  as  follows 

Basic  Gelatin  N  =  1.5426  ±  0.0001 

Gel+ADich  N  =  1.5486 

Gel  +  A  Dich  +  Exp.  N  =  1.5572 

Gel  +  A  Dich  +  Exp.  N  =  1.5515 

+  Dev. 

Gel  +  A  Dich  +  Unexp  N  =  1.5488 

+  Dev. 

From  these,  we  can  see  that  when  properly  processed  holograms  should  have  a 
differential  in  index 

N  ~  1.5515  -  1.5150 
=  0.037  ±  0.005 

with  the  uncertainty  arising  from  the  tolerance  in  the  glass  index. 

In  carrying  out  the  empirical  aspects  of  the  program  dealing  with  the  phenomenology 
of  materials,  it  was  decided  that  with  regard  to  dichromated  gelatin,  commercial  sources 
would  be  employed.  It  was  believed  that  the  consistency  of  base  materials  would  be 
superior  to  that  which  could  be  otherwise  immediately  achieved. 

A  large  number  of  base  gelatins  are  available  but  those  generally  of  interest  to 
other  aspects  of  the  program  are  the  Eastman-Kodak  649F  and  120-02.  Discussions  with 
E-K*  indicated  that  the  emulsions  have  various  initial  degrees  of  hardness  when  made 
available  as  silver  halide  emulsions .  In  decreasing  states  of  hardness  they  are  120-02 
and  649F  with  thickness  typically  6  and  17  pm,  respectively.  After  fixing  however, 
there  is  a  significant  reduction  in  thickness  which  depends  upon  the  silver  halide 
"doping"  that  is 'typically  25%.  Thus,  for  the  same  state  of  humidity,  one  might  expect 
final  thicknesses  to  be  4.5  and  13.5pm,  respectively. 


"‘Eastman-Kodak,  R.  Anwyl,  Sales  Representative 
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There  are  many  processing  combinations  both  for  preparation  and  development 
(e.g.,  Ref  9-12,  5)  which  appear  to  work.  For  the  present  program,  the  primary  con¬ 
cern  is  not  a  general  formulation  but  one  which  specifically  influences  the  creation  of  an 
optical  matched  filter  memory.  On  this  basis,  only  a  few  steps  were  examined  which 
might  prove  vital  to  matched  filter  (MF)  performance.  The  first  of  these  steps  is  the 
pre-doping  fix  where  the  use  of  Eastman  Kodak  standard  mixtures  were  used  on  one  of 
the  plate  types  previously  cited .  In  one  set  of  plates ,  the  hardener  was  omitted ;  in  the 
other  it  was  included.  The  prepared  plates  were  used  to  fabricate  MF's  which  were 
single  and  then  5x  noncoherent  overlays.  Upon  playback  the  results  clearly  show  that 
649F  emulsion  film  ought  not  be  pre-hardened .  This  can  be  seen  in  Table  A-l  where  in 
both  types  of  MF's,  the  autocorrelation  signal  is  decisively  less  in  the  prehardened 
case.  Based  upon  these  results,  plate  preparation  processing  thereafter  eliminated 
prehardening  so  that  the  procedure  consistently  used  throughout  the  program  was  that 
given  in  Table  A-2. 

In  studying  references  concerned  with  processing,  it  was  decided  that  at  least  some 
examination  of  the  developing  procedures  was  also  in  order  if  for  no  other  reason  than 
that  a  standard  be  followed.  Development  was  narrowed  down  to  a  choice  between  the 
more  simple  approach  of  Chang  (Ref  13  and  14)  and  others,  or  the  more  recent  but  more 
involved  technique  of  Graube  (Ref  8).  The  procedure  adapted  to  local  conditions  is 
given  in  Table  A-3 . 

Since  it  was  indicated  in  the  literature  that  the  time  history  of  the  dichromated 
gelatin  was  an  important  influence  on  the  generated  gelatin  quality,  two  tests  were 
established  to  determine  the  influence  upon  matched  filters.  In  one,  the  time  history 
after  removal  from  the  alcohol  wash  was  determined  and  in  the  second,  the  efficiency  of 
gratings  made  at  various  times  after  emulsion  doping,  were  made.  Exposures  in  the 

latter  followed  the  Fibonacci  series,  1,  1,  2,  3,  5,8,  13,  21,  34 . after  the  first 

one.  These  exposures  were  arrived  at  after  a  study  of  earlier  efficiency  results  where  in 
neither  binary  nor  linear  distributions  in  exposure  time  seemed  adequate. 

The  results  of  the  first  test  are  shown  in  Fig.  A-l  and  show  the  time  history  start¬ 
ing  with  the  removal  of  the  dichromated  gelatin  plate  from  the  final  (100%)  alcohol  wash. 

In  another  case  a  plate  was  washed  twice  as  long  as  that  plate  whose  results  are  shown 
in  Fig.  A-l.  This  concern  arose  because  of  the  need  to  delay  final  alcohol  dry  on  the 
one  plate  until  the  test  on  the  other  was  complete.  The  results  show  that  the  two 
plates,  otherwise  similarly  processed,  are  the  same.  Degradation  of  a  filter  was  also 
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TABLE  A-1  EFFECT  OF  PREHARDENING  UPON  CN  IN  M9F-BASED 
AMMONIUM  DICHROMATED  GELATIN  FILM 


STATE 

( f 

ORIENTATION 

INCOHERENT 

OVERLAY 

ir 

64“ 

72“ 

UNHON 

1.00 

IX 

UNHON 

0.816 

0.785 

0.140 

0.433 

1.204 

5X 

HDN 

0.335 

IX 

HDN 

0.224 

0.032 

0.029 

0.252 

0.064 

5X 
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TABLE  A-2  PREPARATION  OF  DICHROMATED  GELATIN  EMULSIONS 
FROM  E-K  120-02  AND  649F-BASED  FILMS 


1.  GENERAL 

A)  USE  CHEMICALS  WHICH  HAVE  BEEN  AT  T  -  20-21°C  FOR  24  HOURS 

B)  DISPOSE  OF  ALL  CHEMICALS  IN  ACCORDANCE  WITH  ESTABLISHED  ENVIRONMENTAL 
CONTROL  PROCEDURES 

Cl  PLATE  EXPOSURE  SHOULD  BE  AT  RH  <  50% 

2.  FIX  WITHOUT  HARDENER 

A)  STANDARD  EASTMAN-KODAK  MIX 

B)  FOR  10  MINUTES  IN  THE  DARK 

Cl  AGITATE  FOR  10  SECONDS  ONCE  PER  MINUTE 

3.  WASH  IN  RUNNING  DEIONIZED  WATER 
A)  FOR  10  MINUTES 

4.  FLUSH  WITH  PHOTO-FLO 
Al  FOR  30  SECONDS 

5.  DRY  IN  DUST  FREE  ENVIRONMENT 

A)  FOR  ABOUT  15  HOURS  AT  ROOM  TEMPERATURE  <~  20°  Cl 

B)  IN  RH  20%  -  30%  (E.G.,  DRIED  COMPRESSED  AIR  ATMOSPHERE! 

6.  IDENTIFY  GELATIN  SIDE  BY  NUMBERING  PLATES  IN  CORNER 

7.  SOAK  IN  5%  AMMONIUM  DICHROMATE  SOLUTION 

A)  FOR  10  MINUTES  IN  THE  DARK 

B)  AGITATE  10  SECONDS  ONCE  PER  MINUTE 

8.  USE  LENS  TISSUE  TO  WIPE  NON-EMULSION  SIDE  CLEAR  OF  EXCESS  DICHROMATE.  WIPE  DRY 

9.  DRY  AS  IN  STEP  5 

A)  VERTICALLY  ORIENTED 

B)  KEEP  IN  DARK  UNTIL  DEVELOPED 

C)  USE  WITHIN  5  DAYS 
1283-010(a)(T) 
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TABLE  A-3  DEVELOPMENT  OF  DICHROMATED  GELATIN  PLATES 

1.  ESTABLISH  SYSTEM  CONDITIONS  FOR  FRINGE  VISIBILITY,  ORIENTATION,  ETC. 

A)  EXPOSE 

B)  KEEP  IN  DARK  AFTER  EXPOSURE 

2.  DEVELOPMENT 

A)  IN  0.5%  AMMONIUM  DICHROMATE  SOLUTION 

B)  FOR  5  MINUTES  IN  RED  LIGHT  OR  DARK 
Cl  AGITATE  FOR  10  SECONDS  EVERY  MINUTE 

3.  FIX 

A)  IN  STANDARD  EASTMAN-KODAK  FIXER  WITH  HARDENER 

B)  IN  RED  LIGHT  FOR  5  MINUTES 

4.  WASH 

A)  IN  RUNNING  DEIONIZED  WATER 

Bl  FOR  10  MINUTES  (CAN  BE  IN  ROOM  LIGHTING  IF  OTHER  CONDITIONS  PERMIT) 

5.  WASH 

A)  IN  PHOTO-FLO  FOR  30  SECONDS 
&  WASH 

A)  IN  DEIONIZED  WATER  FOR  2  MINUTES 

7.  DEHYDRATE 

A)  FOR  3  MINUTES  IN  60/50  ISOPROPANOL/WATER  MIX 

8.  DEHYDRATE 

A)  FOR  3  MINUTES  IN  100%  ISOPROPANOL 

9.  DRY 

A)  FOR  15  HOURS  AT  ROOM  TEMPERATURE  (20° C) 

B)  IN  DUST  FREE  AREA  AT  RH  20%  -  30% 

12S3-014UHT)  _ 


considered.  We  established  an  acceptable  performance  in  efficiency  for  a  plane  grating 
and  then  retained  that  grating  in  a  humid  atmosphere  during  its  test  life.  The  RH  was 
67%  ±  3%.  Efficiency  was  periodically  measured  until  the  efficiency  remained  constant. 
The  results  are  shown  in  Fig.  A-2.  It  should  be  pointed  out  that  the  test  on  the  grating 
was  initiated  immediately  after  the  grating  was  removed  from  a  20%  RH  storage  area. 

Individual  gratings  were  made  at  various  times  after  the  plates  were  doped,  and  the 
efficiency  measured  at  the  indicated  times .  The  start  times  followed  the  Fibonacci 
sequence  given  above.  It  appears  that  within  reasonable  variation  the  efficiency  re¬ 
mains  relatively  constant  at  the  value  it  initially  had  after  processing.  In  order  to  sort 
out  the  vertical  distribution  better,  the  best  plate  of  the  two  plates  at  each  initial  start 
time  was  used,  and  efficiency  plotted  as  a  function  of  Time  after  Doping.  The  basis  for 
this  choice  was  that  the  "best  of  the  two  plates"  would  be  closer  to  the  maximum  capa¬ 
bility  at  the  given  "age."  With  this  rational  we  have  the  results  of  Fig.  A-3.  It  would 
appear  that  there  are  some  times  that  are  more  optimum  than  others  with  one  clear 
choice  -  immediate  exposure. 

The  qualification  on  the  thickness  of  the  medium  as  it  influences  efficiency  and  the 
readout  from  storage  is  discussed  earlier.  In  that  section  and  its  references  it  has  been 
shown  that  thickness  is  a  significant  factor  which  influences  efficiency.  One  direct 
empirical  way  to  determine  the  influence  of  thickness  is  through  an  efficiency  measure¬ 
ment  as  the  grating  in  question  is  azimuthally  rotated . 

This  was  done  with  the  two  plane  gratings  whose  results  are  shown  in  Fig.  A-4  for 
thick  and  thin  representatives,  respectively.  As  a  reminder,  these  were  fabricated  from 
120-02  and  649F  based  plates.  In  conducting  this  assessment,  the  two  wavelengths  6328 
and  4880  were  used.  Also,  the  efficiency  was  determined  by  scanning  from  15  degrees  on 
one  side  of  the  beam  -  plate  normal  condition  to  an  equal  angle  on  the  other  side  with  an 
additional  precision  used  near  plate  normality. 

o 

At  the  fabrication  wavelength  (4880A)  the  azimuthal  sensitivity  is  ±4.5  degrees 
(and  ±4.8°  for  A  =  6328A)  for  the  thick  plate  and  essentially  independent  of  rotation  in 
the  thin  case.  Part  of  the  "skirt"  drop  off  can  be  attributed  to  the  beginning  of 
vignetting  by  the  test  beam. 

Figure  A-5a  shows  a  phase  contrast  photograph  of  a  5x  non-coherently  added  MF 
illustrating  the  surface  relief  while  Fig.  A-5b  illustrates  the  quality  of  a  2.81  pm  grating 
in  the  thick  (649F  based)  DCG. 
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Fig.  A-2  Degradation  of  a  Djchromated  Gelatin  Grating  in  High  RH  Atmosphere 


A-9 


/ 


DIFFRACTION  EFFICIENCY.  * 


A-3  Efficiency  of  Cretins)  (f  =  356  mm1)  in  Dichromated  Gelatin  at  Function  of 
Tima  of  Generation  After  Doping 
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Fig.  A-4  Response  of  Thin/Thick  DCG  Gratings  to  Azimuthal  Rotation 
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APPENDIX  B 


CONSTRUCTION  OF  MULTICAPACITY  MATCHED  FILTERS  USING 
BIREFRINGENT  OBJECT  FILM 


Fabrication  of  matched  filters  for  a  large  capacity  memory  for  correlator 
applications  requires  the  use  of  an  object  film.  The  growing  tendency  in  the  field  is 
toward  polyester  based  film  which  is  significantly  more  birefringent  than  the  older 
acetate  type.  Thus,  care  must  be  exercised  in  using  such  film  in  the  formation  of 
matched  filters,  particularly  when  the  same  object  is  used  and  rotated  in  order  to  have 
360°  object  coverage. 

The  paper  in  this  appendix  addresses  this  problem. 
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Construction  of  muiticapacity  matched  filters  using 
birefringent  object  film 


Robert  A.  Bondurant  and  Kenneth  G.  Leib 


Target  objects  used  to  construct  matched  filters  for  correlator  applications  are  generally  recorded  on  photo¬ 
graphic  film.  Most  commercial  emulsions  of  interest  are  coated  on  birefringent  substrate.  The  authors  ad¬ 
dress  the  effects  of  this  birefringence  upon  the  interference  pattern  recorded  as  a  matched  filter  in  coherent 
optical  correlators.  Results  are  then  expanded  to  include  several  types  of  multiple  matched  filters.  Finally, 
techniques  are  discussed  that  reduce  the  deleterious  effects  of  input  film  birefringence. 


In  coherent  optical  correlators,  a  collimated  laser 
beam  illuminates  a  frame,  or  portion  of  an  input  film, 
which  is  the  subject  of  the  optical  analysis.  This  film 
is  usually  located  one  focal  length  in  front  of  a  Fourier 
transforming  lens,  which  performs  the  transform  op¬ 
eration  and  superimposes  the  transform  upon  a  stored 
Vander  Lugt  type  of  matched  filter  located  in  the  back 
focal  plane  of  the  transform  lens.1  The  signal  emerging 
from  the  filter  is  the  product  of  the  incident  and  stored 
transfer  functions. 

Behind  the  matched  filter,  a  second  lens  takes  the 
inverse  transform  of  the  filter  output  and  displays  the 
squared  modulus  of  this  as  the  peak  correlation  signal 
if  the  film  contains  a  target  corresponding  to  the  one 
stored.  Otherwise,  a  general  clutter  signal  is  ob¬ 
tained. 

Thus,  the  correlator  has  two  basic  operations: 
matched  filter  fabrication  and  correlation  measurement. 
A  matched  filter  is  the  product  of  signal  and  reference 
beam  interference,  thereby  requiring  similarly  polarized 
waves. 

Demonstration  systems  such  as  the  Grumman-de- 
veloped  OMFic2  incorporate  some  means  of  increasing 
the  system  storage  capacity  to  handle  input  scene 
variables.3  In  addition  to  the  multiple  channels  of  the 
multiple  holographic  lens-matched  filter  combination 
in  the  OMFIC,  the  capacity  of  a  correlator  can  be  in¬ 
creased  by  storing  multiple  image  matched  filters  at 
each  address.  Such  complex  matched  filters  in  which 
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images  are  added  at  various  orientations  in  a  coherent 
and  noncoherent  fashion  are  being  investigated.4  Co¬ 
herent  addition  involves  a  multiple-target  input  image 
and  one  exposure  to  construct  the  matched  filter,  while 
noncoherent  addition  involves  multiple  exposures  of 
one  target  at  a  time  to  construct  the  matched  filter. 
Generally  orientation  is  the  parameter,  and  a  single 
input  image  is  rotated  in  an  iso-optical  gate  between 
successive  exposures  for  noncoherent  addition.  An 
example  of  multiple  target  imagery  used  for  preparing 
coherent  addition  matched  filters  is  shown  in  Fig.  I. 

Photographic  emulsions  of  input  imagery  for 
matched  filter  memories  are  generally  prepared  with 
polyester  (polyethylene  terephthalate)  substrates  as 
these  are  commonly  used  by  the  aerial  reconnaissance 
community  due  to  the  high  film  footage  per  pound  and 
dimensional  stability  possible.  Molecular  orientation 
in  the  polyester  sheet  gives  rise  to  birefringence,5  which 
affects  the  polarization  of  the  signal  beam  in  the  cor¬ 
relator,  This  effect  does  not  occur  when  cellulose  ace¬ 
tate  film  base  is  used.  However,  acetate  is  used  with 
decreasing  frequency.  Figure  2  shows  the  birefringence 
of  several  photographic  films  of  interest  as  the  signal 
( through  an  analyzer)  vs  film  angle.  The  analyzer  was 
rotated  for  minimum  output  with  a  vertically  polarized 
laser  input,  and  then  the  test  film  was  inserted  in  front 
of  the  analyzer  and  rotated  (0°:  film  is  vertical). 

From  the  experimental  data  in  Fig.  2,  the  polyester 
based  films  have  a  fast  axis  advance  u>  of  89.2-89.4° 
as 


Calculated  values  of  u>  were  obtained  from 

u>  =  |(2»)/A)(n,  -  n,  W  (21 

Thickness  (  is  0.103  mm  for  t  he  polyester  based  films, 
and  reported  indices5 are  n,  =  1.66and  n,  =  1.64.  The 
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Fig.  1.  Seven-target  input  image  fur  coherent  addition. 
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Fig.  2.  Birefringence  effects  in  commercial  (Kastman  Kodak) 
emulsions. 


result  is  uj  =  91.9°  at  633  nm,  and  thus  the  experimental 
and  calculated  values  for  w  are  in  close  agreement. 

Transmission  differences  at  633  nm  between  649F 
and  SO- 1 73  films  due  to  sensitizing  dye  differences  ac¬ 
count  for  the  disparity  between  the  two  relative  peak 
signals  with  polyester  substrates  shown  in  Fig.  2. 

The  impact  of  birefringence  of  the  input  emulsion  on 
correlator  performance  is  not  significant  in  the  playback 
mode,'’  and,  with  thin  emulsion  matched  filters,  all  input 
films  may  he  used  unless  polarization  sensitive  detection 


equipment  is  used.  Some  effects  may  be  noticed  with 
thick  emulsion  matched  filters  due  to  grating  sensi¬ 
tivity. 

However,  in  fabricating  matched  filters,  the  phase 
information  contained  in  the  interferometric  fringes 
comes  from  the  reference  and  signal  beam  interference 
optimized  for  either  contrast  or  fringe  visibility.  This 
optimization  is  done  at  the  optical  spatial  frequency  of 
interest,  where  the  beams  are  set  typically  for  a  beam 
ratio  of  R  =  1,  where 

reference  beam  intensity  f 

H =  — ;  ;  r  *  l*>) 

signal  beam  intensity 

It  is  apparent  that  the  spectra  of  targets  of  int  erest  are 
highly  variable  in  amplitude  so  that  fringe  visibility  is 
usually  optimized  at  a  specified  spatial  frequency  along 
one  of  the  two  orthogonal  axes  defining  the  optical 
spectrum.1'  The  measure  of  the  optimization  of  the 
similarly  polarized  inteference  of  the  signal  beam  /.s  and 
reference  beam  In  is  the  visibility  function  v,  which  for 
coherent  beams  becomes 

1'  ~  ^  tminl 

=  20hIk)uVUr  +  Is)  =  2/f1/-/(l  +  Kl  (4) 

The  effect  of  the  input  film  on  the  linearly  polarized 
beam  in  the  signal  channel  is  dependent  on  the  amount 
of  retardation  and  orientation  of  the  fast  axis.  The 
intensity  of  the  beam  emerging  from  the  film  plane  that 
results  from  variations  in  the  parameters  can  be  derived. 
The  result  is  that  the  signal  intensity  varies  as 


hi*  5 


-  oos2u>  cos2(y! 
2  \ 


(f>l 


where  u>  is  as  before,  and  (v  is  the  offset  angle  of  the  fast 
axis.  When  w  =  0  (linear  polarization  and  n  =  ir!2),  the 
beams  are  coplanar  and  equal.  The  variations  also 
become 


/.„  =  </„/•■>> 


/  ■>  sin-«/2. 


(61 


which  we  can  use  to  determine  the  visibility  of  fringes 
when  one  of  two  linearly  polarized  beams  is  rotated 
around  the  propagation  axis.  From  Eqs.  (3)  and  (4) 


2|(/»/2)A/»smW2l|l'i’ 

~  I  +  (V2l/(/n  sin  Wit 
=  2  simv/d  +  sin2a).  (7) 

In  general,  the  visibility  becomes 

n  t  V'- 

2 - cos'2<r  cos2uil 

\2  2  ) 

- .  (SI 

1  I  1 

- cos'2i*  cos2w  t 

2  2  / 

so  that  Eq.  (7)  is  an  upper  bound  on  the  fringe  visibility 
for  small  values  of  <v  only.  Many  combinations  of  n  and 
u>  are  possible. 

Not  only  is  the  fringe  visibility  affected  by  the  rota 
tion  of  one  linearly  polarized  beam  with  respect  to  an¬ 
other,  but  the  noninterfering  part  of  the  rotated  signal 
beam  becomes  a  biasing  term  on  the  photographic  plate 
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used  to  record  the  matched  filter.  Thus  for  o>  =  0  and 
any  angle  between  the  beams,  the  bias  intensity  is 

/b,».  =  UJ'2)  sin-(90°  -  <t)  (9) 

for  a  total  signal  beam  intensity  of  IJ2. 

The  net  result  of  the  change  in  fringe  visibility  and 
increase  in  prebias  on  the  matched  filter  is  a  reduction 
in  the  optimization  of  the  matched  filter.  Both  corre¬ 
lation  signal  and  signal-to-clutter  ratio  are  reduced. 
The  multiple  exposure  (noncoherent)  addition  tech¬ 
nique  is  particularly  affected  as  the  input  image  is  ro¬ 
tated  for  each  exposure,  changing  fringe  visibility.  Also 
the  prebias  terms  add  with  each  successive  exposure. 
Figure  1)  shows  a  comparison  between  noncoherently 
added  matched  filters  with  good  fringe  visibility  and 
birefringence-reduced  fringe  visibility  when  tested  by 
playing  a  known  reconnaissance  scene  containing  the 
recorded  target  (M-60  tank)  in  a  terrain  background 


(»)  With  Birefringence-Reduced  Carrier  Fringe  Visibility 


(b)  With  Good  Carrier  Fringe  Visibility 


Fig.  3.  Comparison  of  correlation  planes  using  noncoherently  added 
matched  filters  without  (a)  and  with  (b)  optimum  fringe  formation 
due  to  birefringent  input  film. 


with  false  targets  (other  vehicles).  The  signal-to-clutter 
ratio  is  quite  poor  in  (a),  which  was  made  with  bire¬ 
fringent  input  imagery,  while  the  target  correlation 
stands  out  well  in  (b),  which  was  made  with  good  fringe 
visibility. 

Techniques  must  be  employed  to  minimize  the 
problems  caused  when  working  with  birefringent  im¬ 
ages.  Copying  the  input  film  onto  acetate  based  ma¬ 
terials  overcomes  the  problem,  but  it  is  not  necessarily 
a  desirable  technique.  A  blank  piece  of  processed  film 
may  be  cut  from  the  same  birefringent  stock  as  the  input 
image,  placed  in  an  iso-optical  gate  in  the  reference 
beam,  and  rotated  as  necessary  with  the  input  (signal) 
image  to  yield  best  fringe  visibility.6  Recently,  multiple 
target  imagery  such  as  that  shown  in  Fig.  1  has  been 
used  to  produce  effective  noncoherently  added  image 
matched  filters.  The  image  is  first  rotated  to  align  one 
strand  axis  (fast  or  slow)  with  the  beam  polarization  to 
maximize  fringe  visibility  with  the  fixed  reference  beam. 
Then  a  mask  is  used  to  uncover  one  target  at  a  time  for 
each  successive  exposure  of  the  matched  filter.  The 
final  result  is  a  multiple  overlay  matched  filter,  nonco¬ 
herently  added  with  no  input  image  rotation  required. 
Four-element  matched  filters  of  this  type  have  been 
successfully  recorded  in  silver  halide  using  this  tech¬ 
nique.  Naturally,  the  properly  aligned  multiple  target 
image  may  be  used  for  preparation  of  high-quality  co¬ 
herently  added  matched  filters  as  well. 

As  suggested  by  several  reviewers,  the  degree  of  bi¬ 
refringence  varies  to  a  certain  extent  over  the  area  of 
each  sample  film  due  to  nonuniform  stretching  during 
manufacture.  Such  variations  make  the  previously 
suggested  corrections  only  first-order  approximations, 
which  do  not  totally  solve  the  problems  introduced  by 
the  birefringence.  However,  the  suggested  techniques 
have  resulted  in  significant  improvements  in  the  per¬ 
formance  of  matched  filters  recorded  from  birefringent 
imagery. 

The  research  was  sponsored  in  part  by  the  Air  Force 
Office  of  Scientific  Research  (AFSC),  U.  S.  Air  Force, 
under  contract  F49620-  '8-C-0051. 
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APPENDIX  C 


KOWARSCHIK'S  EQUATIONS  FOR  SEQUENTIALLY  STORED  GRATINGS 


1.  From  "Diffraction  Efficiency  of  Sequentially  Stored  Gratings  in  Transmission  Volume 
Holograms, "  Optica  Acta,  25,  67-81,  1978 


A  =  -2.  p 
]  pz 


u  -  a2  2 
B  =  6  -  P 

C.  =  2  K  .  6 

l 


Di  =  -2i6iz 
e  =  e2  -  6 2 


(6) 


where  k.  =  tt  n./X 

3  =  2tt  n/X 


(7) 


and  n  =  refractive  index  (mean  value)  of  recording  medium 
nj  =  amplitude  of  spatial  modulation  of  i-th  grating 
X  =  vacuum  wavelength 
P  =  reference  beam  wave  vector 
6  =  signal  beam  wave  vector 
2  =  surface  normal 

2.  The  diffraction  efficiency  (  R  )  is  given  by: 

n  =  nl  +  n2  I 


C-l 


/ 


/ 


where  the  abbreviated  terms  are  as  follows: 


Cj  =  -i(u  +  v), 

h  =  T(U  '  v)’ 

A  3  a'i 

v  =  -  yj(u  +  v)  +  j  — , 


vifr-  • 


2,-2 


"1=  *  (“l  *^)+E>2+E22’ 

H3=  t(^2+  Sj^+V-V*ei>. 


(C22  -  jCjV) , 


E,  E_  B 

v  =  -  *  4-  — +  — 

1  DX  D2  A  ’ 


01 2  =  — - - ^EiBD2  "  C1  D2  +  AE1E2  +  BD1E2  "  C2  Dl^  ’  (10) 

AD1D2 


1—  (Cj^j  +  C^j-BEjEj). 


adxd2 


A  1 

Pl  =  U  +  v  -  — 


p2  =  -  y(u  +  v)  +  j  T  (u  -  v)  -  —  ,  l  (ll) 


A  1  /3  “l 

p3  =  --(u  +  v)  -  j  T(U-  v)  -  T* 


1 


U  =  y  [-q  +  /  (q2  +  P3)] , 

v  =  ^  [-q  -  /  (q2  +  P^)l » 


o)2  <*>! 


_  /M3  “j^2  ^3_ 

"  \  3  /  "  6  2 


The  signal  terms  Sj  are  given  by : 


•(Z)|2MC22(^22-v2)V1{Fi2C22Coos  [C«] 


U  z])2  +  (C/jSin  [vz]  +  G.sin  U2z])  >. 
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APPENDIX  D 


THE  FOURIER  TRANSFORM  OF  COHERENTLY  ADDED  IMAGES 

RECTANGULAR  COORDINATE  SYSTEM  / 

This  appendix  deals  with  the  Fourier  transform  properties  of  an  N-element  input 
array  which  constitutes  the  input  image  for  the  coherent  addition  of  images  to  form  a 
single,  complex  matched  filter.  The  array  is  formed  by  considering  the  basic  image 
pattern  (e.g.,  an  M-60  tank)  as  being  translated  and/or  rotated  in  a  plane. 

Let  E(x,y)  represent  the  amplitude  transmittance  of  the  single  pattern.  The 
shaded  area  in  Fig.  D-l  represents  schematically  this  2-D  transmittance  pattern.  In 
this  experimental  program,  the  origin  is  chosen  as  the  center  of  the  regularly  placed 
images  of  constant  radii.  Random  arrays  are  random  in  radii  only.  (See  Fig.  3-12).  The 
Fourier  transform  distribution  for  this  single  input  is  labeled  F(  u>x,  Wy),  in  terms  of 
the  spatial  frequency  variables  wx,  u>y 

Let  S-p  (x,y)  be  the  amplitude  transmittance  of  a  translated  pattern,  as  illustrated 
in  Fig.  D-l.  If  xj,  yT are  the  coordinate  distances  through  which  the  pattern  has  been 
translated  (without  rotation),  then  S-p  can  be  expressed  as 

ST(x,y)  =  E((x-xT)  ,  (y-yT))  (D-l) 

In  Fig.  D-l,  the  input  pattern  has  been  rotated  (without  translation)  through  an 
angle  0  r.  The  amplitude  transmittance  Sr  (x,y)  is  given  by 

SR(x,y)  =  E(x'.y')  (D-2a) 


where 


x'  = 


y’  = 


x  cog9R 
-  x  sin0„ 


+  y  sin0R 
+  y  cos0_ 


D-l 


/ 


Or  equivalently 


x  =  x'  cos0R  -  y'  sin0R 


y  =  x'  sin0R  +  y'  cos0R 


(D-2b) 


An  expression  for  the  amplitude  transmittance  (x,y)  for  an  input  pattern  with 
combined  translation  and  rotation  as  in  Fig.  D-2  is  now  sought.  First,  accounting  for 
translation,  one  has 

STR(x,y)  =SR((x-xT),  (y-y,f)). 


From  Eq.  (2),  one  finds  that 

STR(x,y)  -  E(x,y) 


(D-3a) 


where 


x  =  xT  +  x  cos 0R  -  y  sin0R 
y  =  yT  +  x  sin0R  +  y  cos0R 


(D'3b) 


Now  consider  the  amplitude  transmittance  T(x,y)  for  a  collection  of  N-elements 
which  are  translated  and  rotated  replications  of  the  particular  input  pattern .  Let  S- 
(x,y)  be  the  transmittance  pattern  of  the  i-th  element;  i.e., 

®j(x»y)  =  (STR  (x,y))i 


Thus 

N 

T(x,y)  =  £  S.(x,y) 
i=l  1 


he  Fourier  transform  of  this  sum  is 


F(  U!  ,  w  )  = 


N 
£ 

i=l  J- 


S.(x,y)  e 


-j(w  x+u)  y).  , 
jv  x  y'^dxdy 


(D-4) 


CD-5) 


D-3 
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Using  Eq.  D-3,  this  reduces  to 


F(  to  ,  to  ) 
x  y 


N  -jto  X  -jto  Y 

2  e  x  e  y  1 


00 

yYi  \\  E(X.,Y.). 

-  CO 

e-ix.((oxcose.  +  o)ysine.), 

e'iyi^wxsin0i  +  \cosQi>' 


d^dy. 


(D-6) 


where 


*i=<xTVyi  =  <yTVei  =  ‘VV 


To  simplify  the  notation  set 


a.  =  to  cos0.  +  to  sin0. 
i  x  l  y  i 

0  i  =  -  to  sin0.  +  w  cos0. 
x  i  y  i 


Finally,  one  has 


Ja  -j(to  X.  +  to  Y.) 

P<W=I  6  y‘F(aj,B1) 


(D-7) 


(D-8) 


relating  the  Fourier  transform  of  the  array  to  that  for  a  single  input  pattern.  (Note 


that  at  <ox  =  0=  toy 


F(O.O)  =  ^  F(0,0)  =  N*F(0,0) 


D-5 


2 

Hence,  the  corresponding  DC  light  intensity  or  exposure  is  proportional  to  F*  (0,0)  =  N 
FF*(0,0) ,  where  F*  is  the  complex  conjugate  of  F.  If  the  N-input  images  in  the  array 
are  processed  in  sequence  rather  than  simultaneously,  the  resulting  DC  light  exposure 
is  the  sum  ^ 

F*(0,0)  =  ^  FF*(0,0)  =  N*FF*(0,0).  ) 

i=l 


Symmetries  in  the  input  transmittance  pattern  carry  over  to  the  Fourier  transform 
function.  For  symmetry  about  the  y-axis,  one  has  E(-x,y)  =  E(x,y)  and  consequently 


F(-w  ,oj  )  =  F(;d  , 
x  y  x 


u 

y 


) 


Similarly,  for  a  pattern  that  is  symmetric  about  the  x-axis  one  has  E(x,-y)  =  E(x,y), 
and  then 

F«V  '“y)  =  F  (wx’“y) 


Finally,  for  symmetry  with  respect  to  both  axes,  (i.e.,  with  respect  to  the  origin),  one 
has  E  (-x,-y)  =  E  (x.y)  and 

F(- u  ,  -u  )  =  F(oj  ,w  ) 
x  y  x  ’  y' 


Tables  D-l  to  D-3  illustrate  the  specific  parametric  values  and  the  transform  for 
N  =  2,  4,  8  in  rectangular  coordinates.  The  transform  expressions  are  omitted  for  N  =  10 
in  Table  I>3  for  simplicity. 

POLAR  COORDINATE  SYSTEM 

For  input  functions  that  possess  radial  symmetry,  E(x,y)  =  E(  p  ),  where  p  is  the 
radial  polar  coordinate.  In  this  case,  one  is  only  concerned  with  translational  effects 
since  the  transmittance  function  is  insensitive  to  rotation;  i.e.,  Sjjfx.y)  =  E(  p  )  as 
well.  From  Eq.  (D-l) 


S-p  (x,y)  =  E(  p  ) 


(D-9a) 


TABLE  D-1  EXAMPLES  OF  SYMMETRICALLY  ORIENTED  PATTERNS  WITH 
RECTANGULAR  SYMMETRY 


N 

x' 

D 

Q 

“i 

0j 

F  <^x.wyl 

2 

mm 

D 

0° 

wx 

"v 

EXP  (ju>y  d)  •  F  (wx.  «y)  + 

Si 

D 

180 

'“x 

•wv 

EXP  1-jWy  d)  •  F  <-wx,  -wy> 

3 

0 

-d 

0 

wx 

U) 

V 

EXP  (j  wy  d)  •  F  (wx,  wy)  + 

v/3  d/2 

d/2 

30 

l/3«K*uy|/2 

<-wx  +  v/3  Wy)/2 

EXP  <-jd(v/3<->x  +  wy)/2)- 

-v/3  d/2 

d/2 

120 

l-v/3  <*»x  +  <*>v)/2 

<-ux-v/3  «y)/2 

F  1  [  v/3  wx  +  wy ]  12,  l-wx  +  v/3  wy |  /2)  + 

EXP  (-jd  (-/3ux+tjy)/2i- 

F  (|-v/3u>x+wy)/2,-tcJx+  v/3wy)/2) 

4 

0 

-d 

0 

u> 

X 

wv 

EXP  IjWy  d)  F  (wx,  wy)  + 

d 

0 

90 

“v 

•wx 

EXP  l-j  wx  d)  F  |uy,  -ujx)  + 

0 

d 
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•wv 

EXP  l-j  Uy  d)  F  l-wx,  -Wy)  + 

-d 

0 
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V 

wx 
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8 

0 

-d 

El 

u? 

X 

“v 

EXPljuiy  d)  F(u>x  Wy)  + 

d/v/2 

-d/v/2 

45 

(u>x  +oiy|/v/2 

(-wx  +  u»y)/v/2 

EXP  H  wxd)  EXP  IjWyd)  F  ( <*>x  +  wy>  •  u>x  +  wy ) 

d 

0 

90 

WV 

■"x 

v/2  v/2  v/2  /2 

d/v/2 

d/v/2 
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(-Wx  +  Uyl//2 

-<WX  +Wy)//2 

+  EXPI-i  wx  d)  F  (wy  -wx)  + 

0 

d 

180 

■"k 

■"v 

EXP  (-j  <*>x  d  |  EXP  t-i«yd)  F  (-wx  +  u>y  -(wx  +  wyM 

v/2  v/2  v/2  v/2 

+  EXP  l-j  w„  d)  F  Ujx  -Wy)  + 

-d/v/2 

-d 

d/v/2 

0 

225 

270 

-|uix  +U)y)/J2 

-Wy 

(WX-Wy)//2 

wx 

-d/2 

-d/v/2 

315 

(WX-Wy)//2 

(u>x  +  Cvjy  )  /  v/ 2 

EXP  |juxd)  EXP  (-jWyd)  F  Kwx  +  u>yl,u)„ -wy ) 
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|  1283-01 1(»)(T) 

D-7 


/ 


i 


TABLE  D-2  EXAMPLES  OP  SYMMETRICALLY  ORIENTED  PATTERNS  WITH 
RECTANGULAR  SYMMETRY  FOR  N  -  10 


N 
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*i 
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d 
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V 
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-0.8090  ui  .  0.5878  w 
x  V 
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|  1283-01 2(1  )(T)  | 
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TABLE  D-3  EXAMPLES  OF  SYMMETRICALLY  ORIENTED  PATTERNS 
WITH  RADIAL  SYMMETRY 


N 

eii 

*Ti 

F  <w,  D/FM 

2 

d 

OP 

d 

180P 

2coi  ludcos  () 

3 

d 

0 

/  \  /  \ 

d 
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d 
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4 

d 

0 

d 

90 

2  col  fud  cos  {)  +  2co»(u)d»in 

d 
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d 
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8 

d 

0 

d 

45 

2  co«  (cud  cm  {)  +  2  cm  (ok 1  «in 

• 

90 

+  4  co»  (ojd  cm  tl\f2)  •cdelhidein 

• 
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• 
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where  p  ,  tp  are  polar  coordinates  centered  about  the  point  (x-p.  Yx^»  and 


x-x^,  =  P  COS  l|/ 


y-yT  =  p  sin  ^ 


x_  =  pT  cos  , 


(D-9b) 


Figure  D-3  shows  these  coordinate  variables. 

The  amplitude  transmittance  for  a  collection  of  N-elements  is  given  by  Eq.  (D-4). 
The  corresponding  Fourier  transform  of  this  sum  for  N- radially  symmetric  inputs  is 


N 


F(“x,Wy) 


-j(u)xX  +  wyY) 


dxdy 


£  \\  E(x'xi’  y'yi)  6 

i-1  -oo 

V  -i<%xi  +  “vYi)  ft 

=  Z  e  ))  E(IV  diV*  V 

-  oo 

e']^i(u)xC08  ^i  + 


i=l 
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using  the  i-subscript  notation.  With  the  substitutions  for  spatial  frequency  variables  in 
polar  form 


u 


X 


u>  cos  E, 


w. 


to  sin  £  , 


the  last  exponential  term  in  Eq.  (D-10)  becomes 
exp(-jp.  to  cos(ip.  -  5)). 


The  integral  over  the  -variable  is  a  Bessel  function  identity 

2r 


Jc(p,  »>  =  { « 'ipt  “  co8(Js  •  °  d? 


(D-ll) 
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Thus,  Eq.  (D-10)  becomes 
N 


F(w 

Because 

thus 


i  8  'll“*xi  +  W  /e(0 i>  Si  v-l  “)  dSj 

i=l 

of  radial  symmetry,  the  integral  is  identical  for  all  values  of  1; 
C  -JwpT.(cos  5  cos  ^Ti  +  sin  5  sin  tTi> 


F(u>,  O  =  F(  u>)  2-' 
i=l 


or  finally 

xV  -jwP-ijCOSff;  -  !/>-,•) 

F(0),  o  =F(U))2-  e  T1  T1 


N 


i=l 
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Table 


D-4  illustrates  the  specific  parametric  values  aud  the  transient,  for  N  =  2-4. 


8  with  polar  symmetry. 
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APPENDIX  E 


ADDITIONAL  SIMULATIONS 

During  the  course  of  this  contract  we  developed  a  digital  simulation  of  many  aspects 
of  the  OMFIC  device.  This  digital  simulation  allowed  parametric  studies  to  be  carried 
out  rapidly  and  with  repeatability  and  control  not  achievable  in  the  laboratory.  This 
appendix  contains  two  papers  which  describe  in  detail  two  such  studies.  The  first 
paper  is  concerned  with  the  sensitivity  of  matched  filters  to  different  but  typical 
terrains.  It  shows  the  possibility  of  quantifying  the  motion  of  clutter  using  cross¬ 
correlation  as  a  measure  of  clutter.  In  the  second  paper  a  study  of  some  digital 
preprocess  techniques  is  carried  out.  Also  the  rotational  sensitivity  of  the  matched 
filter  combined  with  various  preprocessing  methods  is  examined. 
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Digital  analysis  of  the  effects  of  terrain  clutter  on  the 
performance  of  matched  filters  for  target 
identification  and  location* 

J.  Mendelsohn,  M.  Wohlers,  end  K.  Leib 

Research  Department 
Grumman  Aerospace  Corporation 
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Abstract 

A  digital  analysis  of  the  effect  of  terrain-target  interactions  on  the  signal-to- 
clutter  (S/C)  ratio  of  the  output  of  a  matched  filter  is  made.  Three  terrain  types  (wood¬ 
land,  village,  and  roadside)  were  included  in  the  investigation.  A  photographic  image  of 
an  M-bO  tank  used  as  the  target  was  digitized  (512  x  512  pixels)  and  a  matched  filter  con¬ 
structed  digitally,  using,  a  two-dimensional  FFT  algorithm.  The  matched  filter  was  then 
used  to  obtain  the  correlation  between  the  target  and  the  various  terrain  types  to  obtain 
the  si pnal-to -clutter  ratio. 

Various  high  frequency  versions  of  the  matched  filter  were  systematically  investigated 
for  possible  optimization  in  the  S/C  ratio  in  an  attempt  to  optimize  the  design  which  was 
eventual lv  realized  in  analog  optical  form.  The  computational  aspects  of  the  problem, 
from  the  point  of  view  of  doing  digital  image  analysis  on  a  minicomputer  with  limited 
memorv ,  is  also  discussed. 


Int  roduct ion 

Optical  matched  filtering  tcchnioues  have  been  under  investigation  for  a  number  of 
years  as  a  means  for  detecting  known  targets  in  various  terrain  backgrounds.  The  feasi¬ 
bility  of  this  approach  rests  in  part  on  the  ability  to  store  many  different  targets  (and 
different  views  of  the  same  target)  in  a  holographic  memory  and,  equally  important,  on  the 
ability  to  search  this  memory  rapidly.  The  system  depicted  in  Figure  1  and  patented  by 
A  Or ume t  allows  for  both  the  storage  of  a  large  number  of  matched  filters  and  the 


Holographic  Lenses 

simultaneous  interrogation  of  all  the  matched  filters  This  is  accomplished  by  means  of  an 
arr.iv  of  holographic  lenses  which  take  its  Fourier  t ransf orm ,  repl ica te  it  for  identification 

This  work  was  partially  supported  hv  the  Air  Force  Office  of  Scientific  Research  under 
font ract  K49620- 78-0-005 1 


190  SPIE  Vo!  186  Digital  Processing  of  Aerial  Images  11979 ) 


DIGITAL  ANALYSIS  OF  THE  EFFECTS  OF  TERHAIN  CLUTTER  ON  THE  PERFORMANCE  OF 
MATCHED  FILTERS  FOR  TARGET  IDENTIFICATION  AND  LOCATION 

and  position  the  replicated  transforms  at  corresponding  positions  in  the  matched  filter  or 
memory  plane.  The  correlation  between  the  object  and  the  images  scored  in  the  memory  is 
then  obtained  from  the  inverse  Fourier  transform  lens.  The  question  <f  the  capacity  of  this 
type  of  system  becomes  a  question  of  the  physical  dimensions,  in  ft-  holographic  memory, 
that  the  matched  filter  occupies,  which  depends  upon  the  wavelength  of  the  optical  system, 
the  focal  length  of  the  system,  and  the  spatial  frequency  bandwidth  of  the  matched  filter. 
Techniques  that  would  increase  the  number  of  matched  filters  that  could  be  stored  in  any 
fixed  sized  holographic  memory  were  previously  discussed . (2 ) 

The  purpose  of  this  paper  is  to  present  the  results  of  an  investigation  that  was 
directed  toward  the  question  of  how  much  space  must  be  allocated  to  each  single  matched 
filter.  The  study  was  conducted  using  a  digital  simulation  of  the  analog  optical  holo¬ 
graphic  system  and  expresses  the  "size"  of  the  matched  filter  in  terms  of  the  spatial 
bandwidth  of  the  filter  needed  to  achieve  a  given  level  of  target  detectability  in  the 
presence  of  the  clutter  signals  that  are  generated  through  the  matched  filter  by  the  back¬ 
ground  terrain.  Thus  the  results  presented  will  be  phrased  in  terms  of  the  S/C  ratio  of 
the  peak  values  in  the  output  correlation  plane  of  the  system,  and  the  interest  is  in 
utilizing  the  smallest  matched  filter  bandwidth  to  achieve  a  given  S/C  ratio. 

Description  of  the  Digital  Simulations  Utilized  in  the  Study 

The  holographic  technique  used  for  the  storage  of  the  matched  filters  is  the  usual 
Vander  Lugt  one,  (3)  where  the  intensity  of  the  exposure  of  the  hologram  is  assumed  to 
be  given  by  - 

j  2  n  a  x 

I£  =  I (F  +  Re  ) | 

with  R  the  amplitude  of  the  reference  beam,  ;<  =  sin •.•/'■  with  t>  the  offset  angle  between  the 
reference  and  signal  beam,  X  the  wavelength,  x  the  physical  coordinate  variable  in  the 
transform,  plane,  and  F  the  spatial  Fourier  transform  of  the  stored  image  (i.e. ,  the 
matched  filter).  The  medium  used  to  store  the  hologram  is  then  modeled  by  assuming  that 
its  optical  transmission ,  t,  is  a  linear  function  of  the  illuminating  intensity,  Ip, 
given  by  t 


where  E  is  a  constant.  The  transmission  function  of  the  hologram  will  contain  all  the 
information  in  the  illuminating  beam  if  Ip,  •  Ec ,  otherwise  there  will  be  a  "clipping"  of 
the  i’luminating  beam  and  therefore  distortion  in  the  stored  hologram.  This  "clipping" 
is  in  fact  used  to  select  a  band  of  spatial  frequencies  that  will  be  stored  in  the  hologram. 
With  the  simple  model  described  above  we  have  no  "clipping"  if  max  (Ip)  (|Fi  +  R)‘>  E  , 
some  clipping  if  (|F|  +  R)^  •  E  ,  and  complete  clipping  if  ( i F |  -  R)‘  'E  Now  if  we  c 
select  the  reference  beam  amplitude,  R,  so  that  R  =  lF(f  ) |  =  T,  where  fc  is  some  specific 
spatial  frequency,  and  in  addition,  adjust  the  exposure  Pime  so  that  (R  ¥  T)2  =  E  ,  then 
we  get  no  clipping  or  distortion  for  spatial  frequencies  such  that]F(f).  |F(f  5  and 
complete  clipping  if  (|F|  -R)2  (|F|  -  T)^  '  E  =  (T  +  R ) 2  =  4t2  or  !F|  3T.  We  see  that 

there  is  a  transition  region  between  T  and  3T  where  there  will  be  a  nonlinear  distortion 
of  the  stored  transform  F.  The  model  which  was  adopted  neglects  this  transition  region 
and  assumes  that  the  transmission  of  the  hologram  is  directly  proportional  to  F  for  all 
spatial  frequencies  such  that  |F(f)l  *■  T  =  R,  and  that  the  transmission  is  zero  for  all 
other  spatial  frequencies.  Thus  when  another  image  is  transformed,  say  to  produce  G(f) 
and  then  used  to  illuminate  the  hologram,  the  resulting  correlation  component  is 
obtained  from  the  the  product  G*(f)  F(f)  where  F(f)  is  a  modified  version  of  the  transform 
of  F(f)  as  indicated  above;  specifically,  it  is  a  band  pass  version  of  F(f)  in  that  only 
those  spatial  frequencies  for  which  |F(f)|  •  | F ( f  )|  will  yield  nonzero  outputs.  Figure 

2  illustrates  typical  results  obtained  from  such  §  model  applied  to  a  digitized  image  of 
an  M-60  tank  (taken  on  a  512  x  512  pixel  format  with  approximately  36  pixels  across  the 
tank  width).  These  images  correspond  to  the  reimaged  versions  of  modified  transforms  of 
the  tank  (set  G*(f)  =  1).  The  three  different  filters  were  obtained  by  selecting  three 
different  frequencies  at  which  the  threshold  T  was  placed.  Some  of  these  thresholds 
correspond  to  the  amplitude  of  the  minor  lobes  in  the  Fourier  transform  of  the  tank  along 
the  frequency  axis  corresponding  to  the  small  dimension  (width)  of  the  tank.  A  plot  of 
the  transform  along  this  axis  is  shown  in  Figure  3.  The  correspondence  between  these 
images  and  those  obtained  in  an  analog  optical  system  is  excellent.  Figure  4  then  contains 
the  corresponding  avit ocorre lat ion  functions  obtained  with  these  matched  filters 
(setting.  G*  =  F- )  .  Again,  the  correspondence  with  the  analog  system  is  excellent;  parti¬ 
cularly  the  dominant  diamond  shape  at  the  center  of  the  correlation  signal  for  all  filters. 
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c.  5TH  LOBE  FILTERED  TANK  IMAGE  d.  VERY  HIGH-FREQUENCY  LOBE  FILTERED 

TANK  IMAGE 

Fig.  2  Filtered  Tank  Images  Obtained  From 
the  Simulation  of  the  Holographic 
Matched  Filter 


Fig.  I  Digital  Simulation  of  the  Amplitude  of  the  M-60  Tank 
Spatial  Fourier  Transform  Along  the  Frequency  Axis 
Corresponding  to  the  Small  Dimension  of  the  Tank 
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c.  5TH  LOBE  FILTERED  TANK  d.  VERY-HIGH  FREQUENCY  FILTERED  TANK 


Fig.  A  Autocorrelations  Obtained  from  the 
Simulated  Matched  Filter  of  Fig.  2 

The  results  shown  in  the  previous  figures  were  obtained  from  a  tank  image  that  was 
digitized  on  a  512  x  512  pixel  format;  subsequently,  all  the  two-dimensional  Fourier 
transform  operations  were  carried  out  using  a  two-dimensional  Fast  Fourier  Transform  (FFT ) 
algorithm.  The  simulation  was  implemented  on  a  minicomputer  system  consisting  of  a  Data 
General  Eclipse  S/200,  with  32K  memory,  a  10  megabyte  disc,  magnetic  tape  for  long  term 
storage  of  the  images,  and  a  scan  converter/TV  monitor  output  from  which  the  images  that 
appear  in  the  figures  were  photographed. 

Simulations  of  Target  Detection  in  Various  Terrain  Backgrounds 

The  performance  of  various  modifications  of  the  matched  filter  were  then  studied 
using  both  the  previously  discussed  thresholding  and  by  restricting  the  total  spatial 
frequency  band  used  (by  simply  eliminating  all  spatial  frequencies  above  some  given  cutoff 
frequency).  Three  different  terrain  views  were  utilized  to  ascertain  the  amount  of 
clutter  or  false  target  response  that  would  be  produced  in  the  various  matched  filters. 
Figure  5  shows  the  variation  of  the  magnitude  of  the  Fourier  transform  of  these  three  scenes 
(again  obtained  with  a  512  x  512  discrete  FFT).  Note  the  village  scene  and  roadside  scenes 
have  greater  high-frequency  content  than  the  woodland,  as  might  be  expected  since  these 
scenes  have  greater  numbers  of  sharp  edged  discontinuities  in  contrast.  Various  matched 
filters  of  the  M-60  ’■an1,  were  simulated,  and  the  scenes  were  passed  through  these  filters 
to  determine  the  cross  correlation.  T ' gure  6  shows  a  cumulative  histogram  of  the  magnitude 
of  the  two-dimensional  cross  correlation  between  the  three  scenes  and  a  single  3rd  lobe 
matched  filter  with  the  output  normalized  to  the  maximum  value  of  the  autocorrelation  of 
the  tank.  (The  peak  crosscorrelation  on  any  one  of  the  three  scenes  was  0  6535  of  the 
autocorrelation  )  The  curves  demonstrate  that  the  correlation  values  or  clutter  signals 
are  higher  statistically  across  the  village  scene  than  the  other  two  and  that  the 
woodland  scene  has  the  lowest  clutter  of  the  three.  As  an  aside,  we  note  that  the  cross- 
correlation  with  the  M-60  tank  is  a  quantitative  measure  of  the  clutter  in  the  scenes 
relative  to  the  tank  as  a  standard,  and  this  quantification  could  be  of  value  to  those 
trying  to  obtain  nonphysiological  measures  of  such  conepts  as  clutter. 

Table  1  contains  a  compilation  of  the  peak  clutter  signals  as  well  as  their  relation¬ 
ship  to  the  peak  autocorrelation  values  from  the  M-60  tanks  for  a  variety  of  matched 
filters.  In  all  cases,  the  matched  filters  were  thresholded  to  the  indicated  suppression 
level  (given  as  a  fraction  of  the  value  of  the  magnitude  of  the  zero  frequency  component 
of  transform  of  the  original  tank),  but  otherwise  all  frequencies  were  present.  We  note 
that  the  zero  order  filter  (unsuppressed)  yielded  unacceptable  performance  in  all  cases 
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Fig.  5 


DISCRETE  SPATIAL  FREQUENCY  (255  MAX) 

Distribution  of  the  Magnitude  of  the 
Fourier  Transform  of  Three  Different 
Scenes  Used  in  the  Study 


Fig.  6  Cumulative  Histogram  of  the  Two- 

Dimensional  Correlation  (C)  Signal 
Obtained  Between  a  Third  Lobe  Matched 
Filter  of  an  M-60  Tank  and  the  Three 
Scenes  (Normalized  to  the  Peak  Auto¬ 
correlation  of  the  Tank) 


Table  1. 

Digital  Simulation  of 
Order  Matched  Filters 

the  Auto  and  Cross  Correlation 
of  an  M-60  Tank  with  the  Three 

of  various 

Scenes 

Terrain 

Lobe 
Or  der 

Suppression 

Level 

Peak 

Auto 

Peak 

Clutter 

S/C 

(DB) 

0 

1.00000 

168. 783 

265.462 

-3.93347 

1 

1.85138  E-02 

30.6461 

11.3396 

8.63555 

loodland 

2 

2.47506  E-02 

35.5385 

16.8633 

6.47514 

3 

2.47455  E-02 

35.4695 

16. 8008 

6.49053 

4 

7. 34271  E-03 

16.9125 

6.45378 

8.36788 

5 

6.90004  E-04 

3.90986 

1.78171 

6.82648 

6 

1.23787  E-03 

5.73741 

2.24171 

8.16271 

0 

1.00000 

168.783 

225.261 

-2.50712 

1 

1. 85138  E-02 

30.6461 

15.3904 

5.98251 

Roadside 

2 

2.47506  E-02 

35.5385 

20.6673 

4.70832 

3 

2.47455  E-02 

35.4695 

20. 7427 

4.65980 

4 

7.34271  E-03 

16.9125 

10.1173 

4.46287 

5 

6.90004  E-04 

3.90986 

1.71455 

7.16022 

6 

1.23787  E-03 

5.73741 

2.91552 

5.87997 

0 

1.00000 

168.783 

248.933 

-3.37507 

1 

1.85138  E-02 

30.6461 

20.2929 

3.58065 

'i  1  lage 

2 

2.47506  E-02 

35 . 5385 

23.2104 

3.70033 

3 

2.47455  E-02 

35.4695 

23.1797 

369496 

4 

7. 34271  E-03 

16.9125 

12.7060 

2.48399 

5 

6.90004  E-04 

3.90986 

2.53437 

3  76583 

6 

1.23787  E-03 

5.73741 

4.21209 

2.68436 

(i.e.,  the  peak  clutter  exceeded  the  peak  tank  signal)  and  that  other  suppressions  yielded 
various  levels  of  acceptable  performance.  Figure  7  displays  typical  performance  variations 
obtained  on  the  village  scene  (as  we  additionally  imposed  a  frequency  cutoff  on  the  filter). 
It  contains  three  of  the  previously  mentioned  M-60  matched  filters, as  well  as  one  additional 
filter  having  a  much  lower  suppression  level,  which  tends  to  pass  only  much  higher  frequen¬ 
cies.  We  see  that  the  very  low  suppression  level  filter  produces  the  best  performance 
(lowest  peak  c lut ter- to-peak  signal),  but  at  the  expense  of  requiring  a  much  higher 
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DIGITAL  ANALYSIS  OF  THE  EFFECTS  OF  TERRAIN  CLUTIER  ON  THE  PtREORMANCE  Of 
MATCHED  FILTERS  FOR  TARGET  IDENTIFICATION  AND  LOCATION 


a.  RESULTS  FOR  A  5TH  LOBE  MATCHED  FILTER 


b.  VERY  HIGH  FREQUENCY  FILTER 


Fig.  7  Performance  of  Various  Matched  Fig.  8  Location  of  Peak  Matched  Filter  Out  - 

Filters  of  an  M-60  Tank  (Used  put  of  Two  M-60  Tank  Filters  on  a 

on  the  Village  Scene)  with  Village  Background  Containing  Four 

Various  Cutoffs  Tanks 

frequency  cutoff;  i.e.,  the  1st,  3rd  and  5th  order  filters  could  have  been  cutoff  at  a 
normalized  frequency  of  140  and  still  produce  essentially  the  same  performance  ratio  as 
would  have  been  achieved  at  a  cutoff  of  255,  but  the  very  high  frequency  filter  requires 
the  full  255  cutoff  to  achieve  the  desired  imnrovement  . 

Finally,  a  simulation  was  conducted  of  the  performance  of  a  variety  of  M-60  matched 
filters  in  locating  tanks  in  the  three  scenes.  Figure  8  and  Table  2  contain  the  results 
of  these  studies  for  the  village  scene  In  Figure  8,  the  location  of  certain  output  cor¬ 
relation  plane  peaks  (local  maxima)  are  shown  in  relation  to  four  tanks  that  have  been 
"digitally”  placed  into  the  scene.  All  peaks  are  shown  whose  amplitude  exceeds  the 
smallest  value  needed  to  locate  all  four  tanks.  A  peak  located  a  tank  if  it  fell  on  the 

tank,  and  once  a  peak  was  found  at  some  point  in  the  scene  all  other  peaks  that  fell 

within  a  region  centered  on  that  point  with  dimensions  of  the  tank  were  not  recorded 
as  additional  peaks.  These  results  demonstrated  that  the  very  high  pass  filter,  which  we 
saw  in  Figure  7  had  the  least  clutter  ratio  and  was  thus  considered  superior  to  the  others, 
did  not  perform  as  well  as  anticipated.  The  reason  for  this  was  that  the  results  quoted 
in  Figure  7  (as  well  as  in  Table  1)  used  the  autocorrelation  of  an  isolated  tank  (white 
background)  to  form  the  ratio  of  signal -to-c lut t er  used  as  the  performance  criterion  The 
alternative  would  have  required  a  computation  to  be  made  with  the  tank  placed  in  a  large 
number  of  positions  in  the  scene  --  clearly  an  excessive  burden  Now  when  the  targets 
are  positioned  in  the  scene,  as  was  done  to  obtain  the  results  given  in  Table  2  and  Figure 

8,  the  effects  of  the  nearby  scene  around  the  tank  are  felt  in  the  peak  correlation  signal 

on  the  tank.  This  is  borne  out  by  the  impulse  responses  of  the  various  matched  filters  of 
the  tanks  (see  Figure  2).  One  sees  that  the  filtered  tank  images  "spread  out"  hevond  t  he 
boundaries  of  the  original  tank,  and  thus  the  nearbv  terrain  will  make  a  contribution  to 
the  correlation  peak  on  the  target.  The  verv  high  frequency  filter  snreads  out  well  hevond 
the  tank  outline  anti  in  a  wav  that  is  very  sensitive  to  the  nearbv  terrain. 
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Table  2.  Simulated  Performance  of  Various  M-60  Matched  Filters 
in  Locating  Four  Tanks  in  Various  Positions  in  the 

_ Village  Scene _ _ _ 

_ _ _ _ Vil  tage _ _ _ _ 

Very  High  Low  Frequency 

1st  Lobe  5th  Lobe  Pass  Filter  5th  Lobe 

Filter  ’  Filter  3  4  _  Filter  3 _ 


Range  of 

Correlation  0.88-0.50 

Signals  on 
Targets  1 

Max.  Clutter  0.47  (.44) 

1,6 


1.31-0.75 

0.45  (.46) 


0.41-0.20 

0.38  (.204) 


1.29-0.60 

1.15 


Error's  in  1/36,  1/36, 
Target  1/36,  1/36 
Locat ions“ 


1/36,  1/36  1/36,  9.1/36  0/36,  1/36 

1/36,  1/36  26/36,  28/36  1.41/36,  1.41/36 


Notes 

1  Expressed  as  the  squared  ratio  of  the  peak  of  the  correlation  signal  on 
the  target  to  the  peak  of  the  filtered  autocorrelation  signal  of  the  target 
on  a  hlanck  background. 

2  Error  in  location  of  the  four  targets  in  the  field  of  view  expressed  as 
fractions  of  the  maximum  width  of  the  target. 

3  1st  and  5th  lobe  filters  were  limited  to  a  frequency  band  that  was  less 

than  8.45  times  the  characteristic  frequency  corresponding  to  the  target  width. 

4  The  very  high  pass  filter  suppresses  all  frequency  components  larger 

in  magnitude  than  1.35  x  10"°  times  the  magnitude  of  the  zero  frequency 
component,  but  it  had  no  frequency  cutoff  as  did  (3). 

5  5th  lobe  filter  with  a  normalized  frequency  cutoff  of  4.23  times  the 
characteristic  frequency  corresponding  to  the  target  width. 

6  Numbers  in  parenthesis  are  the  squares  of  the  ratio  of  peak  clutter 
signals  with  no  targets  present  to  the  peak  of  the  filtered  autocorrelation 
signal  of  the  target  on  a  blanck  background. 

Summary 

Some  selected  results  of  a  digital  s imulat ion , us ing  two  dimensional  FFT  algori thms , of 
an  analog  optical  matched  filter  correlation  system  have  been  presented.  These  results 
demonstrate  that  excellent  correspondence  can  be  obtained  between  simulation  and 
analog  system  by  means  of  a  linearized  model  of  the  holographic  image  storage  medium.  This 
in  turn  allows  a  simple  Fourier  transform  model  of  the  process  to  be  employed  that  can  be 
dealt  with  efficiently  on  a  minicomputer  system  using  a  two-dimensional  FFT  algorithms. 

The  resulting  simulation  enabled  a  study  to  be  made  of  various  technioues  that  would  yield 
a  reduction  in  the  physical  space  allocated  for  the  storage  of  a  holographic  matched 
filter  This  work  is  part  of  a  study  of  means  to  achieve  large  capacity  optical  memories 
for  correlator  applications  and  demonstrates  that  simulations  can  be  effectively  used  in 
achieving  that  objective. 
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Abstract 

A  digital  simulation  of  an  optical  matched  filter  is  described.  Two  important 
aspects  of  the  optical  device  are  investigated  using  the  digital  techniques.  The  first 
study  involves  the  preprocessing  of  the  image  prior  to  performing  the  analog  matched 
filtering  operation  and  the  second  the  sensitivity  of  the  device  to  rotation  of  the  tar¬ 
get.  These  studies  are  carried  out  using  a  military  target  in  common  scenes  such  as 
villages,  woodlands,  and  roads. 

Introduction 

For  the  past  several  years,  the  Research  Department  of  Grumman  Aerospace  Cor¬ 
poration  has  been  investigating  optical  matched  filtering  as  a  means  of  detecting  tar¬ 
gets  in  a  variety  of  military  situations'"^.  The  system  depicted  in  Fig.  1  and 
patented  by  A.  Grumet  allows  for  storage  of  a  large  number  of  matched  filters  and  their 
simultaneous  interrogation.  The  speed  of  such  a  system  makes  it  feasible  to  overcome 
the  traditional  problems  of  the  sensitivity  of  the  matched  filter  to  scale  and  ori¬ 
entation  by  storing  many  versions  of  the  matched  filters  in  various  sizes  and  ori¬ 
entations.  However,  the  optimization  of  the  analog  devices  with  respect  to  its  many 
parameters  (number  of  filters,  spacing,  bandwidth,  etc.)  is  tedious  and  time  consuming 
because  of  the  laboratory  time  consumed  in  constructing  a  single  matched  filter,  and  the 
difficulty  of  precisely  controlling  this  procedure.  With  this  in  mind,  a  digital 
simulation  of  the  optical  device  was  undertaken.  During  the  course  of  the  development 
of  this  digital  simulation  and  subsequent  investigations  it  became  evident  that  digital 
preprocessing  of  the  image  could  possibly  enhance  the  performance  of  the  optical 
devices.  In  a  previous  paper 3,  we  discussed  the  studies  performed  to  modify  the  basic 
matched  filter  and  to  optimize  the  bandwidth  of  the  resulting  matched  filters.  In  this 
paper  we  will  concentrate  on  two  areas  of  investigation:  digital  preprocessing  of  the 
imagery  and  rotational  sensitivity  of  the  filter. 


a. 


MATCHED 

FILTER 


ONLY  TWO  POSITIONS 


DIGITAL  METHODS  IN  THE  DESIGN  OE  OPTICAL  MATCHED  FILTERS  FOR  TARGET  RECOGNITION 


paper. 3  The  basic  tool  used  in  the  simulation  is  a  512  x  512  two  dimensional  Fast 
Fourier  Transform  (FFT) .  Photographs  of  scale  models  of  three  typical  backgrounds 
(woodland,  village,  and  roadside)  and  an  M-60  tank  were  digitized  in  a  512  x  512  pixel 
format  and  stored.  (See  Fig.  2.)  Then  the  FFTs  of  the  stored  arrays  were  taken  and  the 
results  conjugated  to  form  an  approximation  to  the  matched  filter  used  in  the  optical 
method.  Agreement  between  the  digital  matched  filters  constructed  in  this  way  and  the 
analog  filters  were  excellent.  Once  the  digital  representation  of  the  object  have  been 
obtained,  it  is  a  simple  matter  to  perform  controlled  studies  of  the  effect  of  filter 
modifications.  The  limiting  factor  becomes  the  time  to  perform  the  calculation 
(primarily  the  2D  FFT).  In  the  analog  world,  the  construction  of  a  single  matched 
filter  may  take  several  hours  and  it  is  a  relatively  difficult  matter  to  maintain  pre¬ 
cise  control  of  the  process.  The  advantages  of  the  digital  simulation  is  clear,  since 
with  even  moderate  computational  capability  one  can  construct  five  digital  matched 
filters  (512  x  512)  in  an  hour,  with  complete  reproducibility. 
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Criteria  for  evaluation 

In  performing  digital  studies  two  criteria  were  used  in  the  evaluation  of  the  re¬ 
lative  worth  of  a  filter.  The  first  was  the  peak-signal-to-peak-clutter  ratio  (S/C) 
which  was  obtained  by  finding  the  maximum  intensity  in  the  cross-correlation  plane  be¬ 
tween  the  target,  and  the  scenes  without  the  target  in  them,  and  then  comparing  that 
value  to  the  maximum  cross  correlations  between  the  target  on  a  clear  background  and  the 
matched  filter. 

The  second  method  was  designed  to  provide  an  operational  criteria  for  a  particular 
matched  filter.  To  obtain  such  a  measure,  four  tanks  were  placed  at  random  in  a  scene 
(See  Fig.  3).  The  matched  filter  under  study  was  then  cross-correlated  with  that 
scene, and  then  the  resulting  cross-correlation  plane  was  scanned.  A  hit  was  recorded  if 
the  point  of  maximum  intensity  in  the  plane  landed  on  one  of  the  tanks.  A  lethal  radius 
(approximately  half  the  dimension  of  the  tank)  was  then  imposed  so  that  in  the  next  scan 
this  area  was  excluded.  The  process  then  was  repeated  until  all  four  tanks  were  hit  (or 
20  weapons  were  expended).  We  now  proceed  to  the  two  principle  investigations  of  this 
paper;  digital  preprocessing  and  rotational  sensitivies. 


Fig.  3  M-60  tanks  in  three  typical  background  scenes 
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Digital  preprocessing 


It  has  been  established  that  a  matched  filter  thresholded  in  such  a  way  as  to  supress 
amplitude  values  which  are  greater  than  the  value  at  the  5th  lobe  of  the  Fourier  trans¬ 
form  of  the  target,  perform  relatively  well^.  Such  a  filter  can  be  realized  in  an  en¬ 
tirely  analog  fashion  by  using  photographic  techniques  to  suppress  the  lower  lobes.  We 
did,  however,  have  evidence  that  processing  the  scene  before  applying  the  matched  filter 
could  improve  on  this  performance  without  compromising  the  speed  of  the  analog 
filter^.*.  The  important  requirement  is  that  the  preprocessing  could  be  performed  in 
real  time.  A  3  x  3  convolution  such  as  the  Laplacian  represented  by  the  mask 

0-10 
-1  4  -1 

0-10 

was  chosen  as  a  candidate  after  screening  several  other  possibilities.  The  procedure 
was  to  apply  the  mask  to  the  M-60  tank  and  to  the  scenes  and  then  take  the  Fourier 
transform  of  the  resulting  images.  Using  the  Laplaced  tank  as  the  basis  for  the  matched 
filter,  the  two  previously  discussed  criteria  were  applied.  Two  versions  of  the 
Laplacian  preprocessing  technique  were  evaluated.  In  the  first  method,  the  Fourier 
transform  of  the  absolute  value  of  the  Laplacian  image  was  used.  This  we  designate  as 
the  magnitude  Laplacian.  The  second  technique  did  not  attempt  to  change  the  negative 
values  of  the  Laplaced  image  and  we  denote  it  as  the  signed  Laplacian.  The  magnitude 
Laplacian  is  of  interest  because  it  presents  a  somewhat  easier  problem  of  implementation 
to  the  analog  device. 

In  Table  1 ,  the  results  for  three  filters  on  each  of  three  scenes  is  given  with  re¬ 
spect  to  the  S/C  criteria.  In  general,  the  results  indicate  that  the  woodland  presents 
the  least  problem  with  respect  to  clutter  and  the  village  the  most  severe.  The  only 
contradiction  is  the  signed  Laplacian  which  does  equally  well  in  the  woodland  and 
village.  In  any  case,  the  signed  Laplacian  does  uniformily  better  than  the  other  two 
filters  and  the  magnitude  Laplacian  is  the  worst  of  the  three. 

_ _ Table  1  The  Results  of  Three  Filters  on  Three  Scene 


Signal/Clutter 


The  signal  is  the  cross  correlation  between  the  5th  lobe  filtered 
tank  and  the  unfiltered  tank. 

(2)  The  signal  is  the  autocorrelation  of  the  Laplaced 


Table  2  shows  how  each  of  the  techniques  performed  against  the  scenes  containing  the 
4  targets.  We  see  that  the  signed  Laplacian  and  the  5th  lobe  filter  perform  equally. 


No .  of 
We 


Woodland 


Avg . 


No.  of 


No .  of 


5th  Lobe 
Magnitude 

Laplacian(3)  14  0  8  1.0  4  1.0 

Signed 
Laplacian (3 


(1)  The  number  of  weapons  that  would  be  required  to  score  hits  on  all  4  targets. 

(2)  Miss  distance  in  pixels  averaged  over  the  4  targets.  The  targets  dimensions  were 
36  pixels  wide  and  80  pixels  long. 

ace  opera t 


In  all  three  scenes,  the  tanks  were  hit  with  the  first  four  weapons  with  the  miss 
distance  averaging  around  1  pixel.  The  miss  distance  refers  to  the  Euclidean  distance 
in  pixels  between  the  desired  aim  point  (in  our  case  a  point  on  the  turret  of  the  tank) 
and  the  point  on  the  tank  at  which  the  hit  occurred.  The  performance  of  the  magnitude 
Laplacian  is  again  inferior  to  the  other  filters.  Fourteen  scans  (weapon  firings)  of 
thecross-correlation  plane  were  required  to  score  hits  on  the  four  targets  in  the 
woodland  scene  and  8  scans  were  required  for  the  Roadside.  However,  the  magnitude 
Laplacian  does  as  well  as  the  other  two  in  the  village.  This  is  quite  surprising  since 
the  village  scene  yields  the  smallest  S/C.  While  the  overall  performance  of  the  5th 
lobe  and  signed  Laplacian  filters  were  better  than  the  magnitude  Laplacian,  the  behavior 
of  the  latter  in  the  viLlage  scene  warrants  further  investigation.  For  the  remainder  of 
this  paper,  only  the  first  two  filters  will  be  studied.  To  this  point,  only  the  S/C 
allows  us  to  discriminate  between  the  remaining  two  filters,  and  with  respect  to  that 
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criterion  (Table  1)  the  Laplacian  clearly  dominates.  If  it  necessary  to  threshold 
the  correlation  plane  at  some  intensity  such  that  returns  le- ;  than  that  intensity  will 
be  considered  clutter,  then  the  filter  with  the  higher  S/C  should  perform  better.  We 
now  turn  our  attention  to  the  problem  of  the  rotational  sensitivity  of  the  two  filters 

Rotational  Sensitivities 


The  sensitivity  of  the  matched  filter  to  a  change  in  the  target  orientation  has  been 
a  problem  in  the  application  of  this  technique.  The  optical  device  attempts  to  overcome 
this  by  storing  many  matched  filters  of  the  target  in  different  aspects  and  performing 
multiple  simultaneous  matched  filtering  operations.  However,  the  fewer  the  filters 
required  the  more  cost  effective  the  overall  device,  so  that  the  question  of  rotational 
sensitivity  remains  one  of  interest. 

We  first  examined  the  effect  of  rotating  the  filtered  tank  and  correlating  it  with 
the  unrotated  tank.  this  is  in  effect  the  correlation  that  would  take  place  if  the 
target  tank  had  been  rotated  in  the  scene.  In  the  case  of  the  Laplacian,  we  considered 
two  situations;  the  first,  in  which  both  the  reference  tank  and  the  target  tank  are 
filtered  by  the  Laplace  operator  and  the  second  in  which  only  the  reference  tank  is 
filtered. 

Since  the  candidate  filters  tend  to  act  like  high  pass  filters,  we  include  in  the 
comparisons  the  results  for  the  all  pass  or  DC  filter.  It  is  clear  from  the  results  in 
Table  3  that  the  modified  filters  suffer  significant  losses  due  to  the  rotations.  The 
5th  lobe  filter  and  the  second  case  of  the  Laplace  filter  are  very  similar.  The  first 
Laplace  filter  (Laplace  reference  vs  Laplace  target)  is  significantly  worse  than  the 
other  two.  All  the  filters  are  of  such  low  total  energy  and  high  frequency  that  it 
could  be  expected  that  misalignment  of  the  target  would  cause  significant  degradation  of 
the  correlation  signal.  It  should  be  emphasized  that  the  DC  filter  is  included  here 
only  for  comparison  purposes.  Previous  studies  have  shown  its  performance  to  be 
unsatisfactory. 


Table  3  Cross-Correlation  of  the  Rotated  and  Unrotated  Tank  (dB) 


Rotation 

(deg) 

5th  Lobe 

Signed  Laplacian 

_ CD _  C2) _ 

DC 

0 

0.0 

0.0 

0.0 

0.0 

1 

-0.66 

-  2.4 

-6.9 

- 

0.06 

2 

-1  .3 

-  4.4 

-1  .3 

- 

0.11 

3 

-2.4 

-  7.1 

-2.3 

• 

0.19 

4 

-3.4 

-  9.0 

-3.3 

- 

0.30 

5 

-4.3 

-10.7 

-4.2 

- 

0.39 

6 

-5.0 

-1  1.2 

-5.1 

- 

0.51 

7 

-5.8 

-12.2 

-6.1 

- 

0.63 

8 

-6.4 

-12.7 

-7.0 

- 

0.76 

9 

-6.8 

-13.5 

-7.9 

- 

0.87 

10 

-7.5 

-14.5 

-8.9 

0.88 

(1)  Rotated 

Laplace 

tank  cross-correlated 

with  unrotated 

Laplace  tank. 

(2)  Rotated 

Laplace 

tank  cross-correlated 

with  unrotated 

tank. 

We  next  examined  how  the  various  filters  performed  in  locating  the  tanks  in  the 
scenes.  To  investigate  this  we  made  up  the  modified  versions  of  the  matched  filters 
from  tanks  rotated  through  1-10  degrees  and  correlated  these  matched  filters  with  the 
scenes  containing  the  unrotated  tanks.  We  then  searched  the  correlation  plane  in  the 
manner  described  in  the  section  on  the  evaluation  procedures. 

Table  4  gives  the  results  for  the  woodland  scene.  The  superiority  of  the  second 
Laplace  filter  is  clear.  Through  all  10  degrees  of  rotation  it  remains  quite  stable, 
never  requiring  more  than  5  weapons  to  hit  the  4  tanks  and  holding  the  miss  distance  to 
less  than  3  pixels.  The  5th  lobe  filter  and  the  first  Laplace  perform  about  equally. 

The  first  Laplace  is  unsatisfactory  at  5°  and  the  fifth  lobe  at  7°. 

In  the  village  scene  both  Laplace  filters  outperform  the  5th  lobe  filter  which  is 
unstable  even  at  the  very  small  angular  rotations.  The  performance  of  the  5th  lobe 
filter  in  the  village  could  be  expected  because  of  the  very  low  S/C  reported  in  Table  1. 

The  second  Laplace  is  quite  a  bit  better  than  the  first  Laplace  with  respect  to  the 
average  miss  distance,  although  at  the  very  h'gh  rotations  it  does  require  more  weapons 
to  hit  the  4  tanks.  However,  for  rotations  1  ;ss  than  8  degrees  it  is  better  with 
respect  to  both  miss  distance  and  weapon  requirements. 
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DIGITAL  METHODS  IN  THE  DESIGN  OF  OPTICAL  MATCHED  FILTERS  FOR  TARGET  RECOGNITION 


Table  4 

Rotational 

Sensitivity  in 

Woodland  Scene 

b  m  m 

Signed  Laplace(l) 

Average 

No",  o  l 

Average 

No.  of 

Average 

Weapons 

Miss  Dist. 

Weapons 

Miss  Dist 

Weapons 

Miss  Dist. 

- - 

- 775 - 

4 

.  5 

4 

1.57 

4 

.25 

4 

.  5 

4 

2.42 

4 

.25 

5 

.75 

4 

2.42 

4 

.25 

5 

9.7 

4 

2.65 

4 

0.5 

10 

7.4 

4 

2.65 

6 

6.8 

12 

16.8 

4 

2.65 

10 

.5 

9 

17.3 

4 

2.42 

15 

9.5 

10 

14.6 

4 

2.65 

12 

22.15 

8 

10.2 

5 

2.65 

12 

19.12 

12 

16.1 

5 

2.65 

■ J  J  i  x  a-*. 

Table  5 

Rotational 

sensitivity  in 

Village  Scene 

■—HFTTT.  UTirBEWif 

No  .  of 

Average 

No .  Of 

Average 

No .  of 

Average 

Weapons 

Miss  Dist. 

Weapons 

Miss  Dist 

.  Weapons 

Miss  Dist. 

- - 

“T.3 

4 

1.1 

4 

1  .4 

5 

1  .0 

4 

1.1 

4 

1  .4 

15 

1  .0 

4 

1  .4 

4 

1  .4 

7 

1  .0 

5 

11  .1 

4 

1  .4 

20*(2) 

1  .4 

16 

8.6 

4 

1  .4 

20* (1 ) 

41  .4 

6 

23.7 

6 

1  .4 

20* (1 ) 

1  .4 

14 

14.4 

9 

1  .3 

20*(0) 

— 

9 

12.7 

12 

1  .4 

20*  ( 1 

22.0 

17 

27.6 

20*(3) 

1  .4 

20*0 

12.0 

1  1 

15.3 

20*(3) 

1  .4 

(1)  Laplace  operator  appLied  to  reference  target  and  scene. 

(2)  Laplace  operator  applied  to  scene  only. 


*  Search  program  concluded  after  20  weapons  were  fired.  Numbers  in  parenthesis 
_ represent  number  of  hits  at  that  time. _ 

Conclusions 

Various  modifications  of  matched  filters  have  been  examined.  A  Laplacian  was  found 
to  yield  excellent  results  with  respecc  to  S/C  and  rotational  sensitivities.  Results  of 
a  simulation  show  that  the  Laplace  filter  performs  best  when  applied  only  to  the 
reference  target  and  not  to  the  scene.  Our  investigation  of  the  optical  matched  filter 
and  its  modifications  for  target  cueing  and  location  are  continuing. 

Acknowledgement 

This  work  was  partially  supported  by  the  Air  Force  Office  of  Scientific  Research 
under  contract  F49620-78-C-0051  . 

References 

1.  Grumet,  A.,  U.S.  Patent  9  3,779,492,  "Automatic  Target  Recognition  System,"  December 
1973. 

2.  Leib,  K.,  Bondurant,  R.,  Wohlers,  M. ,  "Match  Filter  Correlator  Memory  Techniques  and 
Storage  Capacity,"  Proceedings  of  Soc.  of  Photo-Optical  Instrumentation  Engineers,  Op¬ 
tical  Information  Storage,  Vol.  177,  April  1979. 

3.  Mendelsohn,  J.,  Wohlers,  M.,  "Digital  Analysis  of  the  Effects  of  Terrain  Clutter  on 
the  Performance  of  Matched  Filters,"  Proceedings  of  Soc.  of  Photo-Optical  Eng.,  Digital 
Processing  of  Aerial  Images,  Vol.  186,  May  1979. 

4.  Wohlers,  M. ,  Mendelsohn,  J.,  "Tartet  Detection  Using  Hybrid  Digital-Analog  Cor¬ 
relation  Techniques,"  Published  in  proceedings  of  MICOM  Workshop  on  "Imaging  Trackers 
and  Autonomous  Acquisition  Applications  for  Missile  Guidance,"  Nov.  19-20,  1979,  Hunts¬ 
ville,  Ala. 

5.  Pratt,  W.  K.,  "Digital  Image  Processing"  John  Wiley  &  Sons,  New  York,  1978. 


SPIE  Vol  21  8  Devices  ond  Systems  for  Optical  Signal  Procotsing  (1 980)  /  I  S3 

E-14 


APPENDIX  F 


ANALYSIS  OF  AN  OPTICAL  MATCHED  FILTER 
GUIDANCE  SYSTEM  FOR  TACTICAL 
HOMING  APPLICATIONS 


This  paper  was  presented  at  the  AIAA/NASA  Conference  on  "Smart"  Sensors,  NASA 
Langley  Research  Center,  Hampton,  VA,  November  1978.  It  represents  an  extension  of 
the  initial  concepts  of  large  capacity  optical  memories  developed  in  this  program  to 
tactical  homing  applications . 
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30.  ANALYSIS  OF  AN  OPTICAL  MATCHED 
FILTER  GUIDANCE  SYSTEM  FOR 
TACTICAL  HOMING  APPLICATIONS  (U) 

By:  Kenneth  G.  Leib,  M.  Ronald  Wohlers,  Hans  K.  Hinz,  and  Gardner  Moyer+ 


(U)  Abstract  -  This  paper  describes  an  Optical  Matched  Filter  Image  Cor¬ 
relator  (OMFIC)  that  can  identify  targets  and  their  number,  and  provide 
information  about  their  bearing,  positional  coordinates,  range  interval, 
and  velocity.  Thus,  it  has  potential  for  the  missile  or  projectile 
guidance  application.  The  OMFIC  system  can  be  implemented  using  holo¬ 
graphic  techniques  to  achieve  large  capacity  optical  memories,  and  lasers 
to  achieve  coherent  signal  processing.  Real  time  implementation  is  pos¬ 
sible  using  electro-optic  incoherent-to-coherent  image  transducers.  The 
optical  memory  size  and  real  time  transducer  update  requirements  are  two 
of  the  most  important  aspects  of  the  OMFIC  system  when  used  for  guiding 
a  projectile.  An  analysis  of  these  factors  is  given  which  shows  that 
the  memory  bank  requirements  are  manageable  and  within  current  technology 
capabilities.  Other  aspects  of  the  OMFIC  memory  bank  and  transducer  re¬ 
quirements  are'  developed  in  the  paper  including  a  system  simulation  yield¬ 
ing  CEP  performance  information  based  upon  T-6  THAD  test  missile. 


1.  Introduction  (U) 

(U)  The  advent  of  the  laser  has  given  impetus  to  the  development  of 
various  electro-optical  approaches  for  missile  or  projectile  guidance. 

One  of  the  first  was  a  semiactive  guided  projectile  using  a  laser  des¬ 
ignator  (Ref.  1.).  More  recently,  artillery  shells  containing  surveil¬ 
lance  television  systems  have  been  demonstrated  in  a  simulated  flight 
(Ref.  2).  These  efforts  and  others  between  and  since  show  the  serious 
interest  in  projectiles  which  ultimately  lead  to  "fire  and  forget"  type 
of  weapons.  This  is  confirmed  by  the  issue  of  a  specification  for  at 
least  one  such  projectile  (Ref.  3). 

(U)  The  Optical  Matched  Filter  Image  Correlator  (OMFIC)  is  a  system  under 
development  in  which  targets  can  be  identified  and  selected  and  the  re¬ 
sulting  signals  used  to  provide  information  for  terminal  guidance.  Con¬ 
sidering  the  amount  of  information  that  the  system  can,  in  principle, 
reveal  about  one  or  more  targets,  it  would  seem  to  be  a  logical  candidate 
for  an  electro-optic  missile  guidance  system.  (See  Table  l).  In  this 
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(U)  Table  \  -  Information  That  OMFIC  System  Can 

Reveal  About  Targets  (U)  _ 

•  POSITIVE  IDENTIFICATION  OF  TARGETS 

•  NUMBERS  AND  TYPES  OF  TARGETS 

•  BEARING  INFORMATION  *■.  tadpft 

•  POSITIONAL  COORDINATES  OF  EACH  TARGET 

•  RANGE  (WITHIN  AN  INTERVAL) 

•  |F  PROVIDED,  SPEED  AND  DIRECTION  OF 

•  OUTPUT  AVAILABLE  FOR  TRANSMISSION  TO  REMOTE 

STATION 


Figure  1  OMFIC*  System  in  Tactical  Homing  Applications 
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paper  we  summarize  the  results  of  an  analysis  on  two  of  the  chief  element 
of  the  system,  present  some  results  to  date,  and  indicate  several  tech¬ 
nology  areas  under  investigation.  These  all  bear  upon  the  use  of  the 
OMFIC  system  for  tactical  homing  applications.  (Fig.  l) 

2.  Basic  Elements  of  the  OMFIC  System  (U) 

(U)  The  OMFIC  system  has  been  described  in  detail  elsewhere  (Refs  U  and 
5),  but  the  basic  elements  are  shown  in  Fig.  2.  A  noncoherent  image  (the 
target  in  the  terrain)  is  put  on  a  transducer  that  functions  to  convert 
the  image  into  a  coherent  replica  on  a  laser  beam.  One  such  device  is  a 
liquid  crystal  transducer,  although  solid  crystalline  cnes  are  also  avail 
able  (Ref.  6).  A  coherent  beam  is  then  directed  through  the  transducer 
to  pick  up  the  image  and  transfer  it  to  a  multiple  holographic  lens. 

This  lens  performs  a  Fourier  transformation  in  replication  and  directs  it 
to  corresponding  positions  in  the  matched  filter  plane.  Many  views  of 
one  or  more  tragets  are  stored  with  high  density  in  the  matched  filter 
plane,  and  in  some  cases  several  are  stored  in  each  post  ion.  The  input 
scene  is  compared  with  the  stored  information  and  when  a  target  coinci¬ 
dence  occurs,  the  corresponding  correlation  signals  can  be  obtained  by 
taking  a  second  (but  inverse)  Fourier  transformation.  Under  any  circum¬ 
stance,  signals  are  observed  in  the  correlation  plane  and  through  base¬ 
line  subtraction,  or  thresholding,  the  target  correlations  can  be  detec¬ 
ted. 

(U)  A  priori  information  about  the  matched  filter  and  the  subsequent 
detection  of  a  target  give  rise  to  target  identification  and  positional 
information.  This  is  the  case  since  there  is  an  almost  linear  relation¬ 
ship  between  input  and  output  location  and  movement.  Thus  the  informa¬ 
tion  derived  in  the  detection  processor  can  be  used  to  derive  guidance 
information . 

(U)  From  an  optical  point  of  view,  the  matched  filter  memory  and  input 
transducer  are  the  most  important  elements  and  are  the  subject  of  this 
analysis . 

3.  Requirements  on  OMFIC  System  (u) 

(U)  Figure  3  illustrates  the  change  of  image  size  as  a  function  of  time 
for  a  terminal  guidance  scenario  in  which  the  target  is  to  be  initially 
acquired  at  a  range 

RACQUISITION=1‘5km 
where  the  missile  velocity  is  assumed  to  be 

v  =  210  and  300  m/sec 
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OVERLAPPED  MF 
Figure  2  OMF  1C  System 


t 


Figure  3  Image  Changes  as  a 


Function  of  Time 


V/R  =  14 


(U)  From  this  information,  and  knowing  the  discrete  nature  of  a  memory 
element,  the  number  of  matched  filter  "range  bins"  needed  to  provide 
complete  guidance  coverage  to  some  final  range  to  the  target  can  be  de¬ 
termined. 

(U)  Figure  3  also  illustrates  the  changing  requirement  on  image  smear 
as  the  missile  approaches  the  target.  During  flight  the  duration  of  ex¬ 
posure  for  recording  the  image  on  the  transducer  is  comparable  to  that 
for  photographic  exposures  in  aerial  reconnaissance  cameras  without  the 
benefit  of  image  motion  compensation.  The  percentage  enlargement  smear 
is  given  by 


where  V  =  v/R  and  At  is  the  transducer  write  time. 


(U)  Since  there  is  no  a  priori  knowledge  of  a  targets  aspect,  allowance 
must  be  made  in  the  memory  bank  for  all  views  likely  to  be  encountered. 
Figure  U  illustrates  the  variation  in  matched  filter  response  for  an  M-60 
Tank  target  vs  orientation.  This  parametric  sensitivity  like  any  others 
of  concern,  shows  that  the  parametric  limits  vary  according  to  the  degree 
of  discrimintation  desired.  One  often  accepts  a  (-3dB)  variation.  Thus, 
if  it  were  desired  to  identify  an  M-60  among  other  tanks,  an  orientation¬ 
al  sensitivity  of  ±8°  is  the  maximum  for  normal  discrimination  so  that 
=  23  whole  filters  would  be  required  for  a  full  360°  coverage.  In  a 
similar  fashion,  we  could  determine  the  sensitivity  of  the  matched  filter  to 
scale  change,  we  will  define  a  to  be  the  fractional  linear  scale  change 
to  produce  a  (-3dB)  change  in  correlation.  The  number  of  matched  filters 
can  then  be  determined  in  order  to  provide  both  scale  and  orientation 
coverage  from  the  equation 


where  F  i.  the  total  linear  size  variation  over  the  guidance  path  and 
W  is  the  factor  corresponding  to  orientation,  e.g,  W  =  23  in  the  above. 


!* .  Selected  Results  of  Analysis  (U) 

(U)  Using  Eq.  (2),  the  size  of  the  memory  bank  can  be  computed.  The  re¬ 
sults  are  summarized  in  Table  2  for  various  combinations  of  orientation 
and  scale.  It  is  quite  apparent  that  as  the  degree  of  discrimination  in¬ 
creases,  exemplified  by  decreasing  scale  and  orientation  sensitivities, 
so  does  the  number  of  required  match'  filter  positions. 
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autocorrelation 


Figure  4  Representative  Sensitivity  Response 


<U)  Table  2  -  Memory  Requirments  For  Continuous  Range  Coverage 
1.5  KM-  50  M  (U) 


NUMBER  OF  FILTERS 
TO  COVER  SCALE  AND: 

■E22EBB 

♦ 10  ORIENTATION 

1 12:  ORIENTATION 

•  10% 

397 

311 

259 

•  15% 

271 

213 

177 

•  20% 

208 

163 

- - - - - ■ 

136 
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(ll)  The  image  smear  a:-,  a  percentage  of  the  image  size  when  the  shutter 
is  opened  is  given  by  Kq.  (l)  and  lias  been  computed  for  two  velocity  ra¬ 
tios  V  =  v / H  and  A 1  =  LO  and  ‘>0  milliseconds.  Then  two  exposure  times 
were  chosen  because  they  bracket  near  term  anticipated  and  current  trans¬ 
ducer  technology .  The  results  are  shown  in  Fig.  5.  It  is  not  surprising 
that  shorter  exposures  result  in  less  smear,  nor  that  near  impact  the 
smear  is  greatest. 

(U)  In  carrying  out  our  analysis  it  was  initially  assumed  that  either 
the  transducer  cycling  or  the  guidance  update  requirements  would  dictate 
the  overall  system  cycling  time.  For  the  transducer,  the  cycling  time 
was  set  equal  to  the  time  interval,  T,  that  the  target  is  in  the  "final 
range  bin"  matched  filter  with  nominal  range  R  .  'Hie  result  is 


(U)  Using  t.wo  velocity  rat  ios.  V  and  several  parametric  sensitivity  in¬ 
crements  we  obtain  the  results  of  Table  3. 

(ll)  It  can  be  seen  that,  under  the  conditions  given,  the  most  severe 
transducer  requirement,  called  for  is.  akin  to  TV  frame  rates.  As  is 
shown  below,  the  guidance  update  rate  is  less  severe,  so,  overall,  the 
decision  becomes  more  one  of  properly  selecting  memory  elements  since 
neither  fast  transducer  operation  nor  a  large  number  of  memory  positions 
are  necessarily  required. 

(U)  One  method  used  to  assess  the  target  identification  aspects  of  sys¬ 
tem  performance  war-  to  determine  the  probability  of  detection  as  a  func¬ 
tion  of  false  alarms.  This  war  none  by  using  terrain  scenes  such  as 
those  in  Fig,  o.  Then,  the  ( s i gnu  1 /elut ter )  ratio  in  the  correlation 
plane  was.  determined,  and  the  false  alarm  rate  computed  for  decreasing 
threshold  levels..  A  r>  lat  ionship  between  the  target  detection  probabil¬ 
ity  as  a  function  of  threshold  level  and  number  of  target  parameters 
that  are  varying,  statistically  was  derived  (Ref.  T).  Equating  the  two 
by  eliminating  the  threshold  yields  the  target  detection  probability  as 
a  function  of  the  percentage  false  alarms.  Selected  results  from  this 
analysis,  are  given  in  Table  U. 

.  Results  for  Complex  Terrain  Target  (U) 

(ll)  Although  the  present,  analysis  is  for  the  tactical  scenario  involv¬ 
ing  a  military  target,  one  might  consider  an  entire  scene  as  the  target 
(Ref.  (I)  such  as  a  power  plant  region.  One  scene  obtained  from  a  TV 
monitor  is  shown  in  Fig.  Y  along  with  the  matched  filter  and  correlation 
spot.  The  tabular  insert  gives  the  restrictions  on  the  various  para¬ 
meters-  of  concern.  A  memory  bank  providing  for  continuous  coverage  with- 
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RUNNING  TIME,  SECONDS 


Figure  5  Enlargement  Smear  As  A  Function  of  Flight  Time 


(U)  Table  3  —  Transducer  Update  Requirements  (U) 


VELOCITY 

RATIO 

(v/R) 

IMAGE  SCALE 
SENSITIVITY 
INCREMENT 

(.«) 

REQUIRED 

UPDATE 

INTERVAL 

(MSI 

FINAL 

RANGE  •'BIN’' 
INTERVAL 

(M) 

IMPACT 

TIME 

(SEC) 

0  14 

0  05 

30  54 

64  -  70 

7.14 

0.10 

28.69 

29  -  35 

0,15 

27.04 

18  -  24 

0  20 

27.53 

14  -  20 

0.20 

0.05 

31.91 

95  -  105 

5.00 

0  10 

30.00 

46  -  56 

0  15 

25  64 

26  -  35 

0.20 

28.90 

22  -  32 
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in  these  restrictions  would  require  351*  matched  filter  positions,  a  num¬ 
ber  comparable  to  the  requirements  of  a  vehicle. 

6.  Large  Capacity  Optical  Memories  for  Correlators  (il) 


(U)  In  analog  correlators  storage  is  in  the  format  of  Fourier  transform 
holograms,  and  large  numbers  of  them  are  required  so  that  the  question 
of  information  capacity  is  required  but  is  generally  undefined.  We  can 
however,  consider  the  basic  element  of  storage,  the  matched  filter.  It 
is  usually  optimized  at  some  spatial  frequency,  and  to  obtain  a  prescrib¬ 
ed  level  of  correlation,  say  95$,  one  needs  additional  spatial  frequen¬ 
cies.  This  can  be  described  by  some  factor  Nt/0  where  N  is  a  positive 
number  ^1  and  vQ  is  the  optimizing  spatial  frequency.  Also,  filters 
can  be  geometrically  placed  to  maximize  the  use  of  storage  space.  Fur¬ 
thermore,  such  techniques  as  spectrum  bifurcation,  filter  addition  and 
various  angles  of  the  reference  beam  all  help  to  increase  the  overall  ca¬ 
pacity.  Considering  these  factors,  a  capacity  can  be  defined  by  the 
equation 


25  S  P  0 

(Ni»  XF)2 
o 


matched  filters/cm 


where  the  factors  3,  G,  P  represent  spectrum  splitting,  geometrical 
placement  and  overlaying,  and  F  is  the  transform  lens  focal  length.  A 
practical  set  of  values  is  S  =  1  G  =  1.3,  NM  =  8.3  P  =  *+,  X  =  0.^88  ym  and 
F=  25  mm  to  yield  2590  MF/cm2.  Thus,  one  can  see  that  large  capacities 
are  not  out  of  the  question,  and  that  the  250-350  required  for  the 
scenarios  presented  in  this  analysis  are  a  small  fraction  of  attainable 
den  sit  i’-s  . 


Equation  2  does  incidentally  also  apply  to  the  case  of  bit-by-bit 
serial  recording  for  a  computer  usage.  In  that  case  for  G  =  1,  Nw0  = 
9.0n3  and  P  =  1,  a  system  would  yield  10^  bits/cm2  a  value  often  cited 
in  the  literature. 


(U)  There  is  a  practical  aspect  to  this  and  it  implies  that  while  the 
memories  required  for  a  given  scenario  are  not  approaching  theoretical 
limits,  they  are  exceeding  our  ability  to  normally  make  them.  Obviously, 
resort  to  computer  generation  may  be  necessary,  but  for  missiles  a 
ied i rated  system  for  such  a  "production  line"  is  feasible. 
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(U)  Table  4  —  Probability  Of  Detection  of  M-60 
Tank  in  Terrain  (U) 


THRESHOLD 
LEVEL  NO.  1 

THRESHOLD 
LEVEL  NO.  2 

TERRAIN 

%  FA 

%PD 

Bfl 

WOODLAND 

1.00 

m 

10.0 

99.9 

VILLAGE 

1,70 

Ea 

14.0 

99.9 

POWER  PLANT  IMAGE  FROM  TV  MONITOR 

PARAMETRIC  STABILIZATION  RESTRICTIONS 
AFTER  LOCK  ON 


PARAMETER 

INTERVAL 

RANGE 

0  9  -  2.000  km 

ASPECT 

•2  38 

ROLL 

•1 

PITCH 

2 

YAW 

•2 

REQUIRES  354  MF  POSITIONS 


MATCHEO  FILTER 


AUTO  CORRELATION 


Figure  7  Terminal  Homing  Scenario  With  Complex  Target 
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Matched  filter  density  can  be  increased  in  a  number  of  ways; 
these  are: 

•  Coherent  addition 

•  Noncoherent  addition 

•  Interlace 

•  Off-set  reference  beams 

•  Off-set  (signal  and  reference)  beams 

or  some  combinations  of  these,  as  well  as  a  simple  parallelling  of  com¬ 
plete  channels.  Several  of  these  are  illustrated  in  Fig,  8  and  are  self- 
evident  . 

(U)  Each  has  its  particular  advantages  and  disadvantages  both  in  fabri¬ 
cation  and  complexity  of  reading  out.  A  current  study  is  in  progress  di¬ 
rected  toward  optimizing  the  memories  for  Air  Force  applications. 

7.  Missile  Flight  Simulation  and  Guidance  Update  Determination  (u) 

(U)  The  U.S.  Army's  T6  test  missile  was  used  as  a  model  for  our  simula¬ 
tion,  with  data  and  equations  taken  from  Refs.  9-11 •  The  T6  missile  is 
6 5  inches  long,  6  inches  in  diameter  and  weighs  71  pounds  at  launch. 
Ground  launched,  it  has  a  maximum  velocity  of  1300  ft/sec  and  a  range  of 
about  5  miles.  It  can  also  be  launched  from  helicopters,  the  scenario 
used  in  the  OMFIC  analysis. 

(U)  The  simulation  is  two  dimensional,  includes  rigid  body  and  control 
actuator  dynamics,  but  neglects  launch  rail  dynamics,  thrust  accelera¬ 
tion,  and  seeker  head  dynamics.  The  T6  autopilot  is  replaced  by  a  sim¬ 
ple  pitch  damper,  a  gravity  bias  which  is  a  function  of  velocity  and  a 
guidance  law  based  on  proportional  navigation.  The  target  (tank)  is 
modeled  as  a  point  source  moving  at  a  constant  30  ft/sec  velocity. 

(U)  Preliminary  studies  were  made  in  which  the  initial  velocity  (600  to 
1500  ft/sec)  and  initial  range  (3000  to  25,000  ft)  were  varied.  In  all 
cases  the  initial  altitude  was  600  ft  and  the  flight  path  horizontal. 
Using  proportional  navigation  with  exact  line  of  sight  rate  information 
the  missile  is  well  behaved  with  miss  distances  usually  well  below  one 
f  oot . 

(u)  The  seeker  head  for  the  OMFIC  guidance  system,  however,  does  not 
have  exact,  nor  continuous,  line  of  sight  rate  information.  Instead, 
the  processing  of  visual  information  is  completed  every  A  tg  seconds, 
and  because  of  image  blowup,  this  processing  stops  at  a  range  of  R 
feet,  and  the  missile  is  guided  only  by  previously  stored  information. 
This  is  shown  in  Fig.  9.  For  this  and  all  other  scenarios  only  one  ini¬ 
tial  velocity,  800  ft/sec,  and  one  initial  range,  1500  meters,  was 
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studied.  The  seeker  head  measures  the  angle  9  +  9  (Fig.  10),  subtracts 
the  gyro-measured  angle  9  ,  and  stores  the  line  of  sight  angle  .  This 
is  repeated  every  Atg  seconds.  A  simple  parabolic  curve  fit  is  made 
from  three  consecutive  points, 

IP  =  a  +  bt  +  ct^ 

IP  =  b  +  2ct 

O 

and  V  is  extrapolated  and  used  during  the  next  interval.  This  procedure 
is  repeated  at  every  guidance  update.  The  parameters  A t  and  R_  were 
varied  and  the  results  shown  in  Fig.  11  (no  turbulence).  Only  combina¬ 
tions  of  A  tg  and  Rp  below  this  curve  satisfy  the  one  foot  miss  distance 
requirement . 

(U)  To  evaluate  the  missile  in  a  noisy  environment,  turbulence  was 
added  using  a  Dryden  gust  spectrum  (Ref.  12)  with  a  6  ft/sec  RMS  gust 
velocity.  Smoothing  was  added  to  the  autopilot  by  using  an  N-point 
least-squares  parabolic  curve  fit  for  6  .  N  was  varied  from  3  to  11 
(N=3  is  unique,  not  least-squares).  For  each  combination  of  At  ,  R  , 
and  N,  five  runs  were  made,  each  with  a  different  6  ft/sec  RMS  fust 
time  history.  RMS  miss  distances  were  calculated  and  the  same  one  foot 
miss  criteria  was  used.  These  results  are  shown  in  Fig.  11  (turbulent 
atmosphere).  The  number  of  points  N  is  not  very  significant,  with  slight 
deterioration  for  N  greater  than  six.  Even  N  =  3  (no  smoothing)  was 
noted  to  be  as  good  or  better  than  the  larger  values. 

(U)  As  an  alternative  to  proportional  navigation,  pursuit  guidance  was 
also  evaluated.  Pursuit  guidance  is  not  as  good,  however,  it  is  based 
on  line  of  sight  angle  (not  rate)  information.  Numerical  differentiation 
of  V  to  get  f  is  always  risky.  Assuming,  at  first,  exact  and  continuous 
information  f or  ip  and  no  turbulence,  it  was  impossible  to  obtain  a  pursuit 
guidance  system  that  would  satisfy  the  one  foot  miss  distance. 

(U)  A  least-squares  parabolic  curve  fit  is  a  crude  way  to  smooth  dis¬ 
crete  information;  however,  feasibility  has  been  proven.  Should  the 
OMFIC  guidance  system  some  day  reach  the  design  stage,  a  more  sophisti¬ 
cated  method  will  surely  be  developed.  Within  the  limits  of  this  anal¬ 
ysis  then,  it  can  be  seen  that  the  update  intervals  for  small  CEP  (equiv¬ 
alents)  exceeds  those  required  by  the  transducer.  Thus,  the  guidance 
update  is  the  determining  time  constant. 


8.  Preliminary  Configurations  (U) 

(U)  A  preliminary  configuration  of  the  OMFIC  -  guided  missile  system 
concept  is  shown  in  Fig.  12. 


Figure  1 1  RMS  Miss  Distances  With  and  Without  Atmospheric  Turbulence 
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(U)  At  the  acquisition  range  turn-on,  the  master  control  unit  initiates 
a  clock  generator  from  which  all  functions  are  timed.  This  unit  turns 
power  on  and  causes  the  various  sequences  to  start.  These  include  the 
transducer  cycling,  detector  array  scanning,  laser  cycling,  detection, 
and  thesholding  controls. 

(u)  With  all  units  powered  and  cycling,  the  shutter  opens  and  the  trans¬ 
ducer  is  powered  so  that  an  image  is  stored.  The  transducer  is  postu¬ 
lated  as  cycling  at  twice  per  second  based  upon  the  preceding  analyses. 
After  the  read  cycle  is  complete,  the  shutter  is  closed  and  the  laser 
activated.  This  enables  the  stored  image  to  spatially  modulate  the  col¬ 
limated  laser  beam.  The  OMFIC  operation  procedes  as  described  earlier. 
The  image  is  Fourier  transformed  by  the  holographic  transform  lens  which 
presents  the  spectrum  simultaneously  to  all  positions  of  the  matched  fil¬ 
ter  memory  bank.  At  the  latter,  correlation  of  the  image  takes  place  if 
a  tactical  target  of  interest  is  present.  The  inverse  transform  of  the 
product  of  the  matched  filter  OTF  and  image  OTF  is  taken  and  presented 
to  the  detection  array  which  has  been  activated  and  scanned.  Activation 
includes  battery  power  as  well  as  threshold  setting.  Signals  exceeding 
threshold  activate  the  target  decision  logic  where  it  is  determined 
whether  the  target  is  a  primary,  secondary,  or  tertiary  one.  This  deci¬ 
sion  is  based  upon  information  which  has  been  stored  a  priori  in  the 
matched  filter  during  fabrication.  At  the  same  time  that  the  target  de¬ 
cision  process  is  taking  place,  the  positional  coordinates  are  being 
generated. 

(U)  Upon  the  target  position  decision,  the  offset,  or  guidance  require¬ 
ment  to  reduce  offset,  is  determined.  This  information  is  then  used  to 
generate  guidance  power  signals  which  effect  an  appropriate  aerodynamic 
surface  correction. 

(U)  The  OMIFC  system  was  laid  out  in  the  confines  of  the  T6  missile  pre¬ 
viously  considered  in  order  to  establish  the  reasonableness  of  construc¬ 
tion.  This  layout  is  shown  in  Fig.  13  and  is  based  mostly  upon  labora¬ 
tory  parameters  rather  than  trade  off  studies.  For  example,  considerable 
progress  has  recently  been  made  in  extending  the  coherence  length  of 
semiconductor  lasers.  When  suitable,  the  use  of  this  laser  will  have 
significant  impact  in  reducing  the  size  and  cost  of  the  entire  package. 

9 .  Summary 

(U)  The  OMFIC  system  has  reached  a  state  of  achievement  where  its  appli¬ 
cation  to  tactical  homing  applications  can  be  considered  and  their  per¬ 
formance  analyzed.  The  two  most  critical  elements  of  the  OMFIC  technolo¬ 
gy  are  the  memory  and  the  transducer.  It  was  shown  that  memories  are 
manageable  and  the  required  transducer  cycle  times  are  achievable.  Pre¬ 
liminary  sizing  for  the  Terminal  Homing  Accuracy  Demonstration  (THAD) 
missiles  shows  the  OMFIC  system  to  be  a  reasonable  fit. 


Figure  12  OMFIC-Guided  Missile  System  Concept 
SHAPED  CHARGE  ALLOWANCE 


MATCHED  FILTER 


SHAPED  CHARGE 
ALLOWANCE 


Figure  13  Preliminary  Configuration  of  T6  (THAD)  Missile  With  OMFIC  Mounting 
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